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Characterization of Some Commercial Soaps By X-Ray 


Diffraction 


By Gopal S. Hattiangdi? 


The diffraction of X-rays by 29 commercial soaps was investigated. The patterns 


indicate that the soaps are crystalline, and the data have been used to identify the various 


phases in the solid soaps. The single long-spacing value obtained for the toilet, medicated, 


glycerin, coco, and washing soaps is indicative of the existence of a single crystal type in 


these soaps. The exhibition of two or three sets of long-spacing values by the shaving 


soaps is indicative of the presence in these soaps of two or three individual phases. 


The moisture content of these soaps indicates that the toilet and medicated soaps are 


either hemihydrates or monohydrates; the coco soaps, dihydrates; and the shaving soaps, 


monohydrates. 


I. Introduction 


The purpose of this paper is to present X-ray 
diffraction data for numerous types of commercial 
soaps manufactured in India and in the United 
States, and to correlate these with data for pure 
single soaps with a view to obtaining information, 
insofar as possible, regarding the degree of crystal- 
linitv and phase nature of these systems. It is 
hoped ultimately to correlate these X-ray data 
with the physical-chemical properties of aqueous 
soap solutions and their washing (cleansing) action. 
Results pertaining to these latter properties and 
also some information relating to phase nature, 
as revealed by electron microscopy, are presented 
in RP1973 and RP1974 [23]. 

Soaps are the metal salts of the higher saturated 
and unsaturated fatty acids and of rosin acids; 
they are sometimes understood to include salts of 
naphthenic acids and of synthetic “fatty’’ acids. 
Only the alkali and ammonium soaps, and some 


Part of this work was conducted by the author in the Chem 
stry Department, University of Southern California, Los Angeles 
Calif 


Guest worker, Surface Chemistry Section 
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organic-base salts like triethanolamine soaps, are 
freely soluble in water and therefore suitable as 
detergents in aqueous solution. In the case of 
the saturated fatty acids, the property of deter- 
gency of the corresponding soaps is limited to those 
containing at least 8 and not more than’22 carbon 
atoms, 

Several workers [2, 4, 6, 12, 13] ° have empha- 
sized the importance of thermal and mechanical 
treatment given to pure sodium soaps, inasmuch as 
their transformation from one phase to another is 
readily achieved. Although the composition of 
the fatty stock used in the manufacture of soaps 
may vary from batch to batch, and certainly from 
manufacturer to manufacturer, it may be well to 
consider briefly the soap-boiling process so as to 
understand the extent and nature of mechanical 
and thermal treatment that a commercial soap 
receives, 

Soap is obtained by saponification of fats or 
oils, the nature of the fatty stock largely determin- 
ing the characteristics of the resultant soap. The 


Figures in brackets indicate the literature references at the end of this 
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soap-boiling process involves the fact that soaps 
are soluble in dilute aqueous solutions of electro- 
lvtes but become quite insoluble when the con- 
centration of the electrolyte is increased. Conse- 
quently, after saponification of the fat or oil is 
completed, salt is added, which separates the 
system into the hydrous soap and “lye.” The 
actual boiling process requires 2 to 8 days, depend- 
ing upon the method used and the size of the charge, 
and involves at least five important stages or 
changes: (a) killing, (b) graining, (c) washing, 
(d) strengthening, and (3) fitting or pitching. 
Stages (a) and (b) of the soap-boiling process 
involve heating or “boiling” the charge with the 
aid of “open” or “closed” steam; stages (ce) 
and (d) involve successive addition of water and 
vigorous boiling with open steam, whereas in 
stage (e) heating is discontinued, and slow sep- 
aration of neat soap and “nigre’’ takes place. 
Nigre is a dark-colored layer intermediate between 
neat soap and lye and contains a high concen- 
tration of soap and salts. The soap thus obtained 
is vigorously boiled again with steam 
before it is permitted to settle, the latter process 
being carried out at approximately 65° C for 1 or 2 
days, or even longer. The following four layers 
form in the final ‘settle’: (1) a foamy top layer. 
called fob; (2) the main layer of soap, which, in 
the hot fluid condition, is known as neat soap; 
(3) nigre; and (4) a small amount of lye, called 


closed 


piteh-w ater. 

The important types of commercial soaps manu- 
the full-boiling process are milled 
form-made 


factured by 
soaps (toilet and 
soaps (household bar soaps and floating soaps), 
flake soaps (chipped and powdered soaps), and 
spray-crystal- 


shaving soaps), 


(spray-dried and 
Paste, cold-made coco, and liquid 


sprayed soaps 
lized soaps). 
soaps are manufactured by the semiboiling process, 
the cold process, ete., in which the thermal and 
mechanical treatment given to the soaps is differ- 
ent from that in the full-boiling process. 

Several years ago, MeBain [14] presented an 
interpretation of the manufacture of soap by the 
boiling process on the basis of phase-rule diagrams 
derived by him and his collaborators. Since then, 
several workers have investigated systems of 
single and mixed soaps by employing different 





Investigations of binary and con 


techniques. 
merecial-soap systems by X-ray diffraction re- 
vealed that commercial soaps crystallize as solic 
solutions with no evidence of fractional crystalli 
zation into the constituent single soaps [10) 
Ferguson, et al. [12] have shown from X-ray dif 
fraction studies that sodium soaps exist in fou: 
crystalline phases (alpha, beta, delta, and omeg: 

which can be identified by the positions and in- 


When 


two or more phases are present simultaneously, 


tensities of the !ong and short spacings. 


however, the long spacing changes, but the short 
spacings retain their values and seem to be inde- 
pendent also of the fatty stock and moisture con- 
tent. These characteristic short spacings are at 
3.65 and 2.45 A for the alpha form, 2.75 A for beta, 
3.55 and 2.85 A for delta, and 2.95 A for omega 
Ferguson, et al., found that commercial milled 
soap containing SO percent of tallow and 20 per- 
cent of coconut oil can be transformed into (a 
the beta form by milling and plodding operations, 
(b) the omega form by heating in a sealed tube 
to 190° F and slow cooling, and (c) the delta 
form by reworking the original soap at a tempera- 
ture around 50° to 60° F. 


II. Experimental Procedure 


l. Materials Used 


Of the 29 soaps investigated, 9 were manufac- 
tured in the United States and 20 in India. All 
the soaps were purchased in the retail market in 
the form supplied to customers, the soaps manu- 
factured in India being mailed to the author about 
a month before they were X-rayed. 

Of the 9 soaps manufactured in the United 
States, 6 were popular brands of toilet soaps, 2 
were coco soaps, and 1 was a “medicinal” soap. 
Of the 20 soaps manufactured in India, 6 were 
toilet soaps, 2 were glycerin soaps, 2 were “medi- 
cated” soaps, 4 were washing, or laundry, soaps, 
and 6 were shaving soaps. 

The soaps used for X-ray investigation were in 
the “air-dry” state. Their moisture content was 
determined by drying them to constant weight in 
an air oven at 107°C. The results are presented 
on a percentage basis in tables 1 to 4. 
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TABLE 2 Y-ray d ffraction lata for same miscellaneous comme ul towlel soaps 
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TABLE 3 X-ray diffraction data for some commercial washing soaps 


Soaps 


W4 


(see page 339 


Yellow ochre Pale yellow Pale yellow Pale yellow 


Very soft Dry Very dry Soft 


re content 11.46 l 4" 140 17 Ww 


Long spacings 
Number of peaks t ; i 7 ; 2. 
Intensity of third-order peak 15 , —_ Py 6 
Width at half intensity of third-order peak 0.6 0.70 O55 ate 0.70 0 
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2. X-ray Diffraction Technique 


X-ray diffraction patterns were obtained by 
using a North American Philips Co. Geiger-Miller 
X-ray spectrometer with Cu Kae radiations, the 
intensity of the diffracted rays and the diffraction 
angle being recorded automatically on an elec- 
tronic strip-chart recording potentiometer. Sam- 
ples for X-ray examination were packed into a 
circular hole % in. in diamater cut in a 1- by 3-in 
microscope slide, and smoothed off with a glass 
plate. 


scribed elsewhere (21, 22]. 


More elaborate procedural details are de- 


Bragg spacings (d/n), corresponding to inter- 
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Y-ray diffraction data for some commer tal shari ng soaps 
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planar distances of the crystal lattice, were caleu- 
lated from the positions of the peaks in the curves 
according to the relation n»A=2d sin @, using the 
weighted average wavelength of Cu Ka, and Cu Ka, 
radiations from the copper X-ray target tube 
namely, A=1.5418 A. The weighted average was 
used because the positions of the centers rathe1 
than of the maxima of the peaks were determined. 
Measurements were made to 0.01° in the value of 
the diffractron angle 2@. Values of d/n calculated 
from the measurements were reproducible to 0.02 


A at large angles and 0.2 A at angles as small as 
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The long spacings given in tables 1 to 4 represent 
the length c. sin 8 and were determined by aver- 
aging values calculated from the positions of peaks 
corresponding to the third and successive orders. 
Values were not calculated from the positions for 
the first and second order long-spacing peaks, 
because their positions cannot be determined 
accurately by the present experimental set-up. The 
intensities of the short spacing peaks are expressed 
as their amplitudes above the baseline divided by 
that of the most intense peak in the pattern, which 
is frequently the third order of the long spacing, 
or the short spacing around 4.1 A. Relative 
intensities from line to line are directly comparable 
for a given soap pattern but are not comparable 
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Figure 1. X-ray diffraction patterns of some commercial 
soaps 


from soap to soap because of different reference 
intensities in different cases, as well as because 
of the difficulty of reproducing sample surface 
sufficiently closely to obtain reproducible values 
of absolute intensities. 

Data for all the soaps, including the positions 
and relative intensities of the lines, bands, etc., are 
given in tables 1 to 4. However, since the degree 
of crystallinity can best be seen in terms of the 
sharpness of the pattern, representative one-fifth 
reductions of the original curves of intensity versus 
diffraction angle are given for six soaps in figure 1. 
The choice of the soaps and the order of their 
presentation is purely arbitrary. Values of the 
width at half intensity of the third order long- 
spacing peaks are also presented in tables 1 to 4. 
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III. Discussion 


1. Appearance and Origin of the Diffraction 
Pattern 


Examination of a typical X-ray diffraction 
pattern, such as that for soap TS in figure 1, shows 
that it consists of a series of sharp, fairly intense 
peaks between 263° and 10°, a second group of 
sharp, intense peaks between 20° and 30°, and a 
number of fainter peaks (not very evident on the 
reduced drawings) spread over the whole range. 

The peaks at small angles arise from diffraction 
of X-rays by planes of atoms whose separation is 
proportional to the length of the soap molecule, 
e.sin 8. In all the soaps investigated, successive 
orders of this distance, d/3, d/4, ete., appear as 
peaks in the diffraction pattern; d/7, and d/2 gen- 
erally occur at such small angles that values ob- 
tained with the present instrumental set-up are 
not accurate. The relative intensities of the sev- 
eral orders of the long spacing vary from soap to 
soap, being a function both of the erystal type and 
the scattering power of the various atoms for 
X-rays. 

The peaks in the intermediate angle range arise 
from the diffraction of X-rays from planes of atoms 
separated by much smaller distances than those 
described above and are referred to as_ short 
spacings. The numerical values give distances 
between planes in the crystal lattice and can be 
related directly to the distances between the mole 
cules by making assumptions about the crystal 
type. 

Individual peaks vary considerably in’ their 
angular width at half maximum intensity. <A 
good example is to be seen in comparing the peak 
at around 4.1 A for soaps TS and G2 in figure 1. 
A wide peak, such as that for the latter soap, may 
be referred to as a “halo”. The occurrence of 
halos in these patterns is indicative either of a 
molecular arrangement less orderly than that of a 
three-dimensional crystal, very small size of 
crystallites, or small spatial extent of crystallized 
regions within a different phase. 


2. Toilet Soaps 


The toilet soaps are all characterized by rela- 
tively low moisture content, varying from 3.42 to 
6.87 percent. The theoretical hemihydrate com- 
position for sodium stearate, sodium palmitate, 
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and sodium oleate is 2.86, 3.12, and 2.87 percent 
of H,O, respectively. 
lize a stock consisting of two-thirds tallow or high- 
grade fat and one-third coconut oil or other kernel 


Toilet soaps commonly uti- 


oil, and thus have a high palmitate content. It 
seems reasonable therefore to assume that all the 
toilet soaps investigated here are hydrated to an 
extent of about 0.5 to 1 mole of H,O per mole of 
soap. 

For pure single soaps, each crystal modification 
has its own field of stability with respect to hydra- 
tion or moisture. Its composition and tempera- 
ture range, and the possible equilibria with other 
soap phases, constitute the phase diagram for the 
On the contrary, 
commercial soaps are convenient mixtures made 


particular soap-water system. 


from fats and oils containing both saturated and 
unsaturated components of varying chain lengths, 
and thus have complex phase diagrams. It is in- 
teresting, however, to test the applicability of Fer- 
guson’s hypothesis [12] regarding the absence of 
the alpha phase and the presence of only three 
phases (beta, omega, and delta) in commercial 
soaps to the present data. 

The long-spacing value for the toilet soaps 
varies from 40.6 A for T10 to 42.5 A for T5 and 
T9; the mean average for the soaps manufactured 
in India is 41.7 A, and for those manufactured in 
the United States is 42 A. 
values for pure single sodium soaps reported in the 


The long-spacing 


literature are presented in table 5, in which the 
nomenclature used to designate the various phases 
is the same as that adopted by Ferguson, et al. [12]. 
Accordingly, Buerger’s [6] kappa and eta phases 
and MeBain’s [7, 8] gamma correspond to Fergu- 
son's omega, and Buerger’s epsilon and zeta cor- 
respond to Ferguson's beta. 

It is seen that the closest agreement of the long- 
spacing values for the toilet soaps seems to be with 
that for beta sodium palmitate, the values for 
sodium stearate and sodium oleate being rather 
high, and those for sodium laurate and sodium 
myristate being low. It is interesting to confirm 
the presence of beta sodium palmitate in these 
toilet soaps by comparing their short-spacing 
values with those reported by Thiessen and 
Ehrlich [18], Ferguson, et al. [12], and MeBain, 
etal. [15]. It will be evident from table 5 that the 
agreement among the four sets of data is fair: the 
strong lines around 4.3, 3.9, 3.2, 3.1, and 2.7 A 
seem to match pretty well, and so do a number of 
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other weaker lines. (The italicized values in this 
table are for the very strong, strong, and medium 
lines reported by the authors; the other lines are 
weak or very weak.) 
agreement between these values for beta sodium 
palmitate and the short spacings obtained for th: 
toilet soaps. It should be noted that the rathe: 
prominent ‘“‘halo” around 2.95 A for the toilet 
soaps is characteristic of sodium oleate [20]. Th 
halo from 36° to 42° in 2@ (centering around 2.4 A 
seems to be more characteristic of sodium oleat 
than of sodium palmitate {20}. 


There also seems to be good 


TaBLe 5. Short spacings for beta sodium palmitate 
Phiessen Thiessen 
and Ferguson, McBain and Ferguson, McBain 
Ehriich etal.{12] | et al. [15 Ehrlich etal.[12] | et al. [15 
[1s] Is} 
i i i i i i 
‘ Ww 
s yaw 2 sil 
as 5. 28 n” »* > 
41s 2. #2 2 
we 2. 48 20) 
NA s w” 2.3 2.36 
4 2. 1 
3. 48 44 2. 11 
4 4 L@ 
Is 1.78 


In commercial soaps of low moisture content, 
simple cooling of neat soap produces the omeg: 
phase, which, in many cases, can be transformed 
into the beta phase by agitation at temperatures 
where the beta form is stable. The beta crystals 
are rapidly soluble and form lather easily, and 
thus have a direct and significant involvement in 
commercial soaps. 


3. Coco, Glycerin, and ‘‘Medicated’’ Soaps 


Both of the coco soaps investigated have an 
almost identical moisture content, namely, 13.96 
and 13.99 percent, corresponding to about 2 moles 
of water per mole of soap. Their X-ray diffraction 
patterns are also identical in respect to the posi- 
tion, number, intensity, and sharpness of both the 
long- and short-spacing peaks. Both these soaps 
are reported by their manufacturers to have been 
made from pure coconut oil: the small long-spacing 
value of 34 A confirms the high laurate-myristate 
content to be expected from this oil stock; the 
main short spacings also agree reasonably well 
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sodium laurate. 

The glycerin soaps have the highest moisture 
content among these soaps. Soap G2 has a mois- 
ture content of 17.03 pereent and a long-spacing 
value.of 38.9 A, which is about 1 A longer than 
that for GI (37.8 A), whose moisture content is 
15.2 pereent. These values of the long spacing 
are substantially lower than those for Pear’s soap 

17). The long-spacing value for both G1 and G2 
does not correspond to that for any single pure 
sodium soap, and it is therefore presumed that 
these are tallow-rosin soaps of rather indefinite 
fatty acid composition. The agreement between 
the short-spacing values for these two soaps and 
that for Pear’s soap is relatively good. The glye- 
erol halo, which should normally occur around 4.3 
A [17], is not reported for both G1 and G2, because 
it is not clearly resolved; it probably is present, 
but the relatively high intensity of the closely 
situated peaks at 4.6 to 4.7 A and at 3.96 to 4.06 A 
evidently smears the glycerol halo. 

The medicated soaps have a relatively low 
moisture content, corresponding to 0.5 mole of 
H.O per mole of soap. The values of the long 
spacing for all the three soaps lie close to that for 
beta sodium palmitate and, like the toilet soaps, 
many of the strong short spacings correspond to 
those for beta sodium palmitate and a few to 
those for sodium oleate. The general nature and 
characteristics of the X-ray diffraction patterns 
for these soaps are similar to those for the toilet 
soaps. It seems therefore that the incorporation 
of medicative materials does not alter the crystal 
habit of the soap molecule, although their presence 
seems to soften the “feel” of the products. 

The single long-spacing value obtained for both 
the toilet and the medicated soaps is indicative of 
a single-crystal type rather than a mixture of 
individual components, which would exhibit sev- 


eral sets of long spacings. 
4. Washing Soaps 


Washing soap bars are characterized by high 
moisture content of the order of 11 to 17 percent, 
and a value of ¢.sin 8 which ranges from 33.1 to 
37.8 A. One of these soaps (W4 special washing 
soap) exhibits two distinet sets of long-spacing 
peaks corresponding to values of 34.3 and 41.3 A. 
These soaps are evidently derived from tallow 
acids of rather indefinite and variable fatty acid 
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with the values reported by McBain, et al. [1] for 


not correspond with those for any single pure soap 
of the acetic series. 

There is good general agreement between the 
positions, relative intensities, and half-widths of 
peaks in the region of the short spacings for the 
washing soaps considered together. Since lines 
corresponding to 2.75, 2.85, and 3.55 A are missing 
in all the patterns, one is inclined to predict from 
Ferguson's phase-identification scheme [12] that 
the beta and delta phases are not present. On the 
contrary, lines corresponding to 2.45, 3.65, and 
2.95 A are present; although it seems most un- 
likely that the alpha phase is present in’ these 
form-made washing soaps, the chances are that the 
omega phase is present to a predominating extent 
and contributes these above-mentioned lines. 

The omega phase is believed to occur in com- 
mercial-framed and in most milled soaps, the 
proportion present depending on temperature, 
moisture, stock composition, and other processing 
conditions. This phase is obtained ordinarily by 
rapid chilling of neat soap without agitation, and 
always tends to form when low-moisture soaps of 
usual formulas are slow-cooled from the molten 
condition. In general, formation of the omega 
phase is favored by higher temperatures, lower 
moisture content, and low molecular weights. In 
commercial sodium oleate, it is perhaps the only 
stable form under most conditions [12]. 


5. Shaving Soaps 


All of the six soaps investigated were in stick 
form and had a moisture content ranging from 6.7 
to 8.6 percent, which is equivalent to about | mole 
of water per mole of soap. Four of these soaps 
had two separate sets of long-spacing peaks each, 
and two soaps, S2 and Sé, exhibited three different 
sets. The lowest value of ¢.sin 8 among these was 
37.7 A, and the longest was 46.9 A. The most 
intense short spacing is around 4.0 A for each soap, 
and relatively weaker lines are present up = to 
about 2.4 A. 

Shaving soaps are manufactured usually by the 
soap boiler’s full-boiling process and are milled 
soaps. On the basis of the long-spacing values for 
sodium palmitate given in table 6, it seems that the 
shaving soaps are mixtures of two or more phases 


of this soap. The probable phases in each soap 
are given in table 4. The persistent exhibition of 
a long-spacing value corresponding approximately 
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composition, and the values of ¢.sin 8 therefore do 

















to that of alpha sodium palmitate is rather sur- 
prising, since this phase, unlike beta, omega, and 
delta, is rarely encountered in the usual commercial 
soaps. The alpha phase is usually the starting 
point for the formation of the other phases, how- 
ever, and its presence in the shaving soaps is 
indicative of incomplete transformation to the 
other three phases. Transformation of the alpha 
form to the beta form is achieved by removal of 
water either by heating to about 52° C in an open 
or closed tube, by evacuation at room temperature, 
or by grinding at room temperature [4]. Also, 
water can be incorporated with the omega form to 
produce beta, and water can be driven off from 
beta to reduce it to omega. The formation of the 
delta phase is favored by high molecular weights, 
low soap content, and low temperatures. It is 
obtained in sodium palmitate by quenching nigre 
and in sodium stearate by ordinary cooling. 
Delta is also formed spontaneously at room 
temperature in 50 percent aqueous mixtures of 
alpha, beta, and omega palmitates. In addition, 
extrusion of sodium palmitate of high moisture 
content at and somewhat above room temperature 
converts the beta to the delta phase [12]. 


TABLE 6 Literature values of long spacing c.sin B, for 


sodi im 80a ps 


Phase 
Soap = 
Alpha Beta Omega Delta 
| 1. 5 [12] | 
Na-Laurate 1.8 [1] 
l31 » (9, 10] | 
136.0 [10] } 
Na-My t t 12 
136.2 112 
| 7 11.9 [15 9.6 [12] | 
Na-Palmita ‘he 6 [12] s2.1 [12 0.2110 pes 12 
23111 } 
1.4 [12 45.8 [12 44.1 [16] 
51.8 [3 5.2 [16 44.4 [12]. 
51.9 (16) B.3 (a 17 
Na-Stear ate - “ - ~ a ‘I.- 45.0 [12 
52.7 [16 wy 8 44.9 (9, 10] . 
45.0 [16] 
$5.6 [15] 
Na-Oleate 44.9 [20] 


6. Particle Size, Degree of Crystallinity, and 
Crystal Form 


A crystalline soap, well developed in three 
dimensions, exhibits a sharp, rather intense long 
spacing and several short spacings that are well 
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defined and that lead to sharp peaks in the X-r.y 
diffraction pattern. Symmetry of the = cryst 


lattice in certain cases may reduce the number 


short spacings, but the spacings remain sharp 
and fairly intense. Several sharp short-spaci 
peaks combine to give a broad peak of a sort re- 
ferred to as a halo, but when this occurs, it is 
usually found that this peak is more asymmetr\: 
than the halos obtained for liquids and for less- 
crystalline arrangements. 

The decrease in the degree of regularity may 
take place along one dimension at a time. If the 
regularity along the c-axis falls off, the long spacing 
decreases in intensity to the point of disappearanc 
If the decrease in regularity occurs along the othe: 
axes of the crystal arrangement, the long spacing 
remains relatively intense, but the number of the 
sharp short spacings decreases. In the extreny 
case, random orientation of the chains within 
sharply defined layers leads to sharp long spacings 
plus a liquid like halo (broad and nearly symmetri- 
cal) centered at an angle corresponding to thy 
average separation of the hydrocarbon chains 
within the layers. 

Different crystal forms are characterized by 
different unit cells, the most easily observed char- 
acteristic being perhaps the value of the long 
spacing, c.sin 8, from which the angle of tilt, 8 
can be evaluated. The short spacings obtained 
for any arrangement of ions are not so easily cor- 
related with structure. However, it might | 
that the translations along the a- and b-axes of 
the unit cell should essentially be the spacings 
obtained, except for their absence due to extine- 
tions corresponding to symmetry properties of 
the given crystal. 

A decrease in the intensity of the different 
maxima, but with no change in the sharpness or 
number of the peaks, may be interpreted as being 
caused by dilution or by a decrease in the size of 
the crystalline particles present. An increase in 
the half-widths of ordinarily sharp peaks, but 
with no change in their intensity or number, is 
usually caused by the transformation of an orderly 
crystal form into a relatively disorganized liquid- 
crystalline form. .A decrease in the number of 
lines with no change in the intensity, or sharpness. 
of the peaks may be interpreted as being due to the 
transformation of the material from one phas 
state to another. 
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IV. Summary 





1. The diffraction of X-rays by 29 commercial 
soaps, including toilet, glycerin, medicated, coco, 
washing (laundry), and shaving soaps, was in- 
vestigated. The degree of sharpness of the long- 
and short-spacing peaks for the various soap 
patterns, including those for the glycerin soaps, 
indicates that the soaps are crystalline. 

2. On the basis of their moisture content, the 
toilet and medicated soaps were found to be either 
hemihydrates or monohydrates, the coco soaps 
dihydrates, and the shaving soaps monohydrates. 

3. The X-ray diffraction data have been used 
to identify the various phases present in the solid 
soaps. The values of the long and short spacings 
for the toilet and medicated soaps seem to indicate 
that these soaps consist essentially of beta sedium 
palmitate and to a smaller extent of omega sodium 
oleate. The glycerin soaps are evidently tallow- 
rosin soaps of rather indefinite composition; they 
are not amorphous or undercooled liquids but 
contain very fine crystalline material. The coco 
soaps have a high laurate-myristate content of 
the omega phase. The washing (laundry) soaps 
exhibit long spacings that do not correspond to 
values for any single pure sodium soap, but their 
short spacings indicate the presence of theomega 
phase. The shaving soaps are characterized by 
the existence of two or more individual phases, 
probably of sodium palmitate, as evidenced by the 
presence in the X-ray diffraction patterns of two 


or more distinct sets of long-spacing lines. 


The author is very grateful to Robert D. Vold, 
University of Southern California, Los Angeles, 
Calif., for use of the X-ray spectrometer and for 
general discussions on the interpretation of soap 


phases. 
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Characterization of Alkali Soaps by Electron 
Microscopy 


By Gopal S. Hattiangdi' and Max Swerdlow 


The morphological differences between several closely related alkali soaps have been 


determined by the techniques of electron microscopy. The erystalline alkali soaps consist 


of an interlocked mesh of bundled fibers, whereas the liquid-crystalline soaps exhibit a flagellar 


and/or featureless phase. The micellar groupings for the individual soaps are unique, and 


a scheme based upon the electron micrographs is suggested for their characterization. The 


structures of lithium, sodium, and potassium palmitates are closely related. The laurate, 


palmitate, and stearate of sodium also exhibit homologous structures. 


The results of these morphological investigations offer a rapid and accurate means of 


identification and have been successfully applied to the problem of characterizing commercial 


soaps of unknown composition, 


I. Introduction 


When alkali soaps are heated, desiccated, 
hydrated, or ground mechanically, they pass 
through a remarkable series of forms, each con- 
stituting a definite stable phase existing over a 
definite range of temperature and or composition. 
Attempts have been made to characterize these 
phases by means of X-ray diffraction patterns, but 
the data reported in the literature are rather 
anomalous. Buerger, Smith, Ryer, and Spike 
1]}* believe that sodium soaps eXist in seven 
distinet crystalline phases S £h & 6, a, and n); 


‘? 


in addition, they suggest several others (a’, a 


n. 9’, ete), which are closely related to these 
seven. Ferguson, Rosevear, and Stillman [3] 


report only four solid soap phases (a, 8, 6, and @), 
their omega phase corresponding to Buerger’s 
kappa and eta phases and their beta, to Buerger’s 
epsilon and zeta phases. Also, Ferguson’s omeg: 
phase is believed [3] to be identical to the gamma 
phase reported several years earlier by deBrette- 
ville and McBain [2]. Added to this rather con- 
fusing use of different designations for the same 





' Guest worker, Section of Surface Chemistry. 
? Figures in brackets indicate the literature references at the end of this 
paper 
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phase is the large amount of somewhat conflicting 
data reported by these investigators, a general 
idea of which is presented elsewhere {4}. 

The electron microscope, with its high resolv- 
ing power, has been used with considerable suc- 
cess in recent vears to study a variety of col- 
loidal systems. The curd fibers of sodium laurate 
have been investigated by Marton, MeBain, and 
Vold [5] with the aid of an electron microscope, 
and have been shown to consist of a mass of fibers 
that are thin ribbons whose widths are reported 
to be integral multiples of approximately twice 
the length of the soap molecule 

The present investigation is of an exploratory 
nature. It is designed to reveal morphological 
differences between eight pure alkali soaps repre- 
senting different homologous cations and anions, 
and to characterize some commercial soaps of 


unknown composition. 


II. Materials Used 


The soaps investigated were the palmitates of 
lithium, sodium, potassium, rubidium, and ces- 
ium, and sodium laurate, sodium stearate, and 
sodium oleate. All the palmitates were prepared 
by neutralizing a solution of palmitic acid in 
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95-percent ethyl alcohol with an alcoholic solution 
of the corresponding alkali hydroxide and drying 
the product to constant weight at 105° C. The 
sodium soaps were prepared by neutralizing a 
solution of the fatty acid in 95-percent ethyl 
alcohol with sodium ethylate; sodium laurate and 
sodium stearate were dried to constant weight at 
105° © 
duced pressure at 40° C. 
fatty 
tabulation: 


and sodium oleate was dried under re- 
The constants of the 
following 


acids used are given in. the 


E-quiv- lodine Melt- 
Acid alent num- ing 
weight ber point 
( 
Lauric Kastman Whodak 
Company's best quality 202.6 0.05 15. 0 


Palmitie (FRastman Kodak 
Company's best quality 257. 0 0 62. 2 
Stearic (Kahlbaum, Scher- 


ing’s pure quality 284. 3 0 69.0 
Oleie (British Drug Houses 
pure quality 288.3 85.8 


Two toilet soaps (T3 and TS), one manufac- 
tured in the United States and the other in India, 
and samples of washing (W5) and shaving (S5) 
soaps (both manufactured in India) were the same 
as those used for the X-ray diffraction studies [4]. 
They were obtained from the retail market and 
were examined without any purification. Chem- 
ical analysis disclosed that the cation was sodium 


in all of the four soaps. 


Ill. Experimental Procedure 


The soaps were dispersed in anhydrous n-butyl 
alcohol by milling with a mullite mortar and 
pestle for 5 minutes in an air conditioned room at 
24°C and humidity. The 
suspension Was cast on Formvar (polyvinyl for- 


60-percent relative 
mal) plastic films supported by 200-mesh stainless 
steel specimen screens. Upon drying, the speci- 
mens consist primarily of the undissolved soaps, 
since their solubility in n-butyl alcohol is very low. 

Since metallic shadowing adds contrast as well 
three-dimensional surface 
details, a very thin coating of evaporated chro- 


as a aspect to the 
mium was deposited obliquely (at an angle whose 
tangent equals 5 em/20 cm) upon the soap speci- 


ons : hl 
mens. This was achieved in a manner similar 
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to that of Williams and Wyckoff [8] by slowsy 
evaporating 100 mg of pure (electrolytically de 
posited) chromium metal with a current of 30 
amp at 15 v in about 10 min at a pressure o/ 
10-° mm of Hg. The chromium was outgassed 
for about 20 min with 10 to 20 amp of current 
During th 
deposition of the metal, elevations and depression 


through the tungsten-wire basket. 


on the surface of the soap specimen cast shadows 
characteristic of its contour; where no metal is 
condensed, transmission of the incident beam is 
maximum, and thick 
electron scattering is maximum and little or no 


where elevations occur 


electrons are transmitted. More details concern- 
ing the technique of metallic shadowing used here 
are reported elsewhere [7]. 

The shadowed specimens of the soaps were ex- 
amined and photographed with an RCA--EMU 
(50 kv) electron microscope. The sampling plan 
included at least eight separate specimens for each 
soap. Representative micrographs of character- 
istic areas for the individual soaps are presented 
in figures 2 to 13, the linear magnification in each 
cause being denoted below each reproduction 
Negative prints are shown in the figures. For 
proper interpretation of the surface features and 
structures, they may be viewed as though the 
soap formations were obliquely illuminated. 


IV. Results and Discussion 


One of the most characteristic features of the 
soups when viewed through the electron micro- 
scope is that many of them exhibit an interlocked 
mesh of fiber bundles of varying diameters and 
different degrees of twist but with a general tend- 
ency toward retaining both a criss-cross and a 
parallel structure (fig. 7, A). The diameters of 
the fibers are dependent not only upon the nature 
of the soap but also upon the concentration and 
rate of crystallization, inasmuch as denser areas 
of the same soap specimen exhibit stouter fibers 
than those found in regions of lower concentra- 
tion. Although any measurement and 
quent interpretation of the ultimate fiber widths 


subse- 


in terms of absolute units seems questionable, a 
qualitative comparison of dimensions and struc- 
tures between soaps is reasonable. The growth 
of the fibers in parallel rows at right angles to 
each branch, either in the same or in different 


planes, is interesting and is probably related to 
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the monoclinic angle (the angle of tilt of the 
hydrocarbon chains, 8). When 8=90°, preferred 
orientation of the lateral branches is a maximum. 
The packing of the hydrocarbon parts of the soap 
molecules can be achieved in several ways. One 
possible manner, constructed from known values 
for the bond distances, bond angles, and atomic 
radii [6], is shown in figure 1, in which the packing 
of the soap molecules in the fiber is restricted 
to a lateral association of the hydrocarbon chains. 

Another outstanding fact revealed by the elec- 
tron microscope is that each soap exhibits unique 
and distinet features, such as a curdy mass, an 
octopean mass, filamentous, hairy, frond-like, 
and sheaf-like formations. The presence of these 
formations as well as the increase or decrease in 
their degree of crystallinity, internal and external 
structure, ete., are dependent upon the position 
of the soap in the cationic or anionic series. That 
these formations are characteristic and reproduci- 
ble is shown consistently by the presence of 
approximately the same number, size, and shape 
of such formations in different specimens of the 
same soap prepared under identical conditions. 
These are probably due to a type of structural 
unit, such as a micellar grouping, within the soap 
fibers and are probably related either to the mosaic 
structure of the surface of crystals or to the crystal 
of the individual phases. Before any correlation 
of statistical significance can be made with regard 
to external morphology and molecular arrange- 
ments, further studies are necessary. Whatever 
the interpretation, these patterns serve as excellent 
guides for a quick characterization of the alkali 
soaps notwithstanding possible artifacts. 


l. The Alkali Palmitates 


Lithium palmitate is characterized by thick 
curdy clumps. Whether present singly or joined 
together with one or more of such clumps, they 
always exhibit a tendency for the curd to grow 
out ina frond-like formation that is usually opaque 
to the electron beam and lacks fine structure (fig. 
2,A,B, and C). The primary pattern for sodium 
palmitate is revealed in the less dense areas as a 
very fine network of interlaced fibers that show a 
perceptible periodic structure (fig. 3,C). On this 
background are developed mounds similar to those 
for lithium palmitate but with fine spicular fibers 
fanning out from them (fig. 3,B). In the less con- 
centrated regions, these mounds have the con- 
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figuration of the fan palm (fig. 3,A). The re- 
semblance of the sheaf-like appearance found in 
figures 2,.C and 3,A is noteworthy. Potassium 
palmitate, like sodium palmitate, appears basically 
crystalline (fig. 4,B and C), and the fibers show a 
spiral periodicity that is not exhibited by lithium 
palmitate. Furthermore, the fibers tend to pack 
laterally with their periodic structure in register 
(fig. 4,B). In addition to this highly crystalline 
state, this soap exhibits a featureless liquid-like 
phase (fig. 4.4). In contrast to the somewhat 
analogous behavior of lithium, sodium, and potas- 
sium palmitates, rubidium and cesium palmitates 
show almost entirely different habits. Rubidium 
palmitate is characterized by a liquid-like phase 
from which small spheroidal particles seem to sep- 
arate (fig. 5,B), the larger aggregates being present 
as terraced mounds of irregular shape (fig. 5,A 
and C). Cesium palmitate also exhibits a liquid- 
like phase that is present in the form of amorphous 
ribbon-like films (fig. 6,A). Figure 6,B and C 
show thin long fibers, which probably have a 
greater degree of crystallinity, developing from an 
agglomerated mass of this liquid phase. 

In summation, it might be said that of the five 
homologous alkali palmitates, lithium, sodium, 
and potassium exhibit a somewhat similar appear- 
ance, Whereas rubidium and cesium are dis- 
tinetly different. The first three soaps are differ- 
entiated by the mounds of curd-like material 
from which fibers grow. The degree of erystal- 
linity, orientation, and periodicity of these fibers 
are least for lithium palmitate and most for potas- 
sium palmitate. Rubidium and cesium  palmi- 
tates have a predominant liquid-like phase unlike 
the other three soaps, which exhibit a linear net- 
work of interlaced fibers. Rubidium palmitate is 
characterized by very small terraced mounds of 
eccentric shapes, whereas cesium palmitate ex- 
hibits agglomerated masses, from which thin long 
fibers grow out in all directions. 


2. The Sodium Soaps 


Sodium laurate has a basic structure of curd- 
like fibers intermeshed in a more or less regular 
parallel formation (fig. 7,A). The resemblance 
between the unshadowed specimen (fig. 7,B) and 
that published by Marton, McBain, and Vold [5] 
is very interesting because the soap structures seem 
to be independent of the two different methods of 
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specimen preparation and of the ionic or monionic 
nature of the solvents used. In the less concen- 
trated areas, the bundled fibers appear to be 
flattened out into a sheaf-like pattern (fig. 7,C). 
There also appears a thorny protrusion in these 
fibers. Sodium palmitate exhibits finer fibers, 
and sodium stearate has the finest (fig. 8,A), but 
both these soaps have a parallel bundling similar 
to that for the laurate. Sodium stearate has 
extremely fine fibers developed in dendritic forms 
(fig. 8,C). By contrast, the basic structure for 
sodium oleate consists of a very stout fibrous 
network, the less concentrated areas showing a 
fine swirled flagellar formation (fig. 9,A). From 
this predominant phase large micellar groupings 
develop with an octopean structure—a semicrystal- 
line mound in the center and curling fiber bundles 
of varying diameters emanating in all directions 
(fig. 9,B). 

In summation, it might be said that the sodium 
soaps of the saturated fatty acids exhibit a homol- 
ogous behavior, characterized by a meshwork of 
interwoven fibers laid out in a regular linear for- 
mation and a sheaf-like bundling of the finer 
fibers (fig. 8,B, 3,A, and 7,C). The fineness of 
the sheaf-like formations varies for the sodium 
soaps in the order: laurate < palmitate < stearate. 
Sodium oleate exhibits a very different behavior, 
characterized by a lack of linearity and is probably 
due to the essentially liquid-crystalline nature of 
this soap. 

The increase in the periodic spacings of the soup 
fibers with increasing molecular weight of the soap 
(LiP< NaP<KP and NaL< NaP< NaStr) sug- 
gests their possible correlation. This is brought 
out in the tabulation below, in which data are 
presented for the periodicity for the alkali soap 
fibers as found in the original electron micrographs. 


\Molee- Approxi 
Soap ular mate 
weight periodicits 
A 
Sodium palmitate 278. 3 500 
Potassium palmitate 294. 4 675 
Sodium stearate 306. 3 750 


The periodic structure that results from a type 
of stress imposed upon the growth of the soap 
fibers is exhibited by those soaps that are crystal- 
line in nature and are salts of the saturated fatty 
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acids. Other factors that have an importan 
bearing on the amorphous and/or crystalline nature 
of the soaps are the degree of hydration, composi 
tion, temperature, pressure, extent of mechanica| 
treatment, and rate of precipitation. Assuming 
that all these factors are similar for the various 
soaps, it seems reasonable to conclude that for the 
soaps with the same anion, the strain is probably 
caused by the heavy cation attached to the fatty 
acid radical; and for the soaps with the samy 
cation, the strain is brought about by the inertia 
of the heavy hydrocarbon chains (fig. 1). 


3. Application to Commercial Soaps 


Electron micrographs of the samples of two 
toilet soaps exhibit a network of fine interlaced 
fibers (figs. 10, B and 11, C); numerous mounds 
of curd-like material from which fine fibers show- 
ing traces of periodic structure grow out (figs. 
10, A and 11, A); and a swirled flagellar formation 
(figs. 10, C and 11, B). The first two character- 
istics are typical of sodium palmitate and the third 
of sodium oleate. Toilet soaps are usually made 
from oils of high palmitate-oleate content. Anal- 
ysis by X-ray diffraction cannot definitely identify 
the components of the binary system, since only 
an average value for the two long spacings and 
overlapping peaks for the short spacings are ob- 
tained [4]. 
electron microscopy are in excellent 


On the other hand, observations by 
agreement 
with chemical analysis, and in such cases prove 
to be not only more rapid but more accurate than 
X-ray diffraction or other methods of analysis. 
The washing (laundry) soap (fig. 12) is charae- 
terized by an essentially liquid-like phase from 
which aggregates of indefinite shape and little fin 
structure separate. The extremely primitive 
state of these aggregates and the absence of any 
appreciable number of fiber or flagellar formations 
makes identification of this soap uncertain, al- 
though the possible existence of sodium palmitat: 
X-ray diffraction data [4] show that 
the long-spacing value for this soap does not 


is indicated. 


correspond to that for any single pure sodium soap 

The shaving soap consists of a very fine network 
of interlaced fibers, which show a periodicity o! 
about 600 A (fig. 13, C), and curd-like mounds 
from which fine short fibers emanate (fig. 13, A 
and B). The first may be expected from a mix 
ture of sodium palmitate and sodium stearate, anc 
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he second from sodium palmitate alone. X-ray 
diffraction data [4], however, indicate the presence 
of several phases of sodium plamitate and do not 
reveal the presence of any other type of soap. 


The authors are very grateful to Robert D. Vold 
for some of the soap samples, to James W. McBain 
and Howard F. MeMurdie for their critical com- 
ments on the interpretation of this study, and to 
Esther Golovato for her very valuable assistance in 
processing the many electron micrographs. 
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Some Physical Chemical Properties of Aqueous Solutions 
of Soaps and Soapless Detergents 


By Gopal S. Hattiangdi,'! William W. Walton, and James I. Hoffman 


A study has been made of the effect of concentration on the surface tension, electrical 


conductance, pH, opacity, and foam formation of aqueous solutions of 30 ty pical commercial 


soaps and 6 popular soapless detergents. 


The soaps included the following types: Toilet, 


“medicated,” glycerin, coco, washing (laundry), and shaving soaps. The results are explained 


on the basis of the colloidal nature of the solutions and the soap phases that have been found 


in the solid soaps by studies of X-ray diffraction and electron microscopy. The significance 


of the results in terms of micelle formation and the relationship between micelles and deter- 


gent action are discussed, 


I. Introduction 


Running water was the earliest “detergent” 
used, the cleansing being caused mainly by 
mechanical action and partly by solvent action. 
The relatively superior washing action of aqueous 
solutions, suspensions, and pastes of various herbs, 
plants, and nuts was quickly realized, to be fol- 
lowed a few centuries thereafter by the manufac- 
ture and use of saponified oils and fats (soap) 
In more recent times, synthetic soapless detergents 
are being commonly used by the housewife for 
laundering, dish-washing, and cleansing in general. 
Whether it is the ancient herb, the more recent 
soap, or the modern soapless synthetic compound, 
the cleansing action in all cases is based upon the 
same principle —the replacement of the undesired 
dirt by a detergent, leaving the dirt in a suspended 
and protected state so that it is easily removed by 
running water, 

The criteria used in purchasing soaps and soap- 
less washing compounds are usually their appear- 
ance and texture, the quantity of suds they 
produce, and, with some critical housewives, the 
action of the compound on the skin and hands. 
There are, however, no universally accepted 
physical-chemical standards to which soaps are 
referred with regard to their washing power. Some 
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pertinent experimental methods suggest them- 
selves if the prerequisites for efficient washing are 
critically reviewed and analyzed. 

In trying to wash a greasy dish or a dirty fabric, 
a superior wetting or surface-active agent is needed. 
After the greasy dirt is removed from the surface 
of the dish or fabric, the dirt particles must be 
suspended in the aqueous medium. This can be 
achieved by compensating the charge around the 
particles. When the charge is neutralized, both 
the uncharged dirt particles and the clean surface 
of the dish or fabric must be protected with an 
adsorbed film of the detergent so that the removal 
of the dirt with running water may be facilitated. 
The use of an efficient detergent also enhances the 
degree of solubilization of the insoluble dirt, owing 
to its sorption by, or inclusion inside, the colloidal 
particles of the detergent. The factual data on 
these factors may be obtained by making direct 
measurements of surface tension, interfacial tension, 
carbon number, gold number, and dye number. 

Other factors, which may be considered signifi- 
cant by the housewife, are the quantity of suds 
produced (properties of foams), the irritation to 
the skin and hands (control of pH), and the 
clarity of the solution (opacity of turbidity). 

The literature is profuse with physical-chemical 
data for aqueous solutions of pure alkali soaps, 
stress being laid in most of these investigations on 
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the colloid chemical viewpoint. Krafft and Stern 
|4}* suggest that soaps must be in solution in order 
to exhibit appreciable detergent action. Me- 
Bain [6] has shown conclusively that the detergent 
action of aqueous soap solutions is due to their 
colloidal nature, and that this action is influenced 
by the chemical constitution of the detergent only 
insofar as this affects their behavior as colloids. 
Investigations conducted during the past 10 vears 
have shown, however, that sodium soaps of the 
saturated fatty acids display probably three al- 
ternative crystal forms, as well as a sequence of 
about 5 or 10 stable mesomorphic forms, each 
having different and unique physical-chemical 
properties [9]. 

The purpose of this paper is to present data for 
several physical-chemical properties of aqueous 
solutions of numerous commercial soaps and soap- 
less detergent materials, and to correlate these 
with the structure and phase nature of the solid 
soaps as revealed by investigations with X-ray 
diffraction [2] and electron microscopy [3], and to 
interpret the significance of the results in terms 
of fundamental colloid chemical hypotheses. 


II. Materials Used 


Thirty typical commercial soaps (toilet, “ medi- 
cated,” glycerin, coco, washing (laundry), and 
shaving soaps) and six popular synthetic soapless 
detergents were investigated. The commercial 
soaps were the same as those used in previous 
investigations [2, 3]. Of the soapless detergents, 
two (D1 and D5) are reported by their manufac- 
turers to be nonalkaline compounds, one (D6) as 
containing “a special ingredient that dissolves 
grease on contact,” and one (D2) as a new kind of 
soap. The active principles of the soapless deter- 
gents D1, D8, D4, and D5 are probably a fatty 
alcohol sulfate, an alkyl-aryl sulfonate plus an 
aleohol sulfate, a sulfated monoglyceride of 
coconut oil, and an alkyl-aryl sulfonate, respec- 


tiv ely. 
Ill. Experimental Techniques 


Aqueous solutions corresponding to 10, 5, 2.5, 
5 


hed 0.625, 0.313, 0.156, and 0.078 percent, by 
weight, of the soap or soapless detergent were 
prepared by dissolving the appropriate amount 


3 Figures in brackets indicate the literature references at the end of this 


paper 
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of the soap (or soapless detergent) in hot (80° to 


100° C) distilled water. The solution was then 
cooled, transferred to a 250-ml volumetric flask, 
and fresh cold distilled water added until the 
solution reached the standard mark. 


1. Surface Tension 


A du Nouy tensiometer was used for making 


surface-tension measurements. The platinum 
ring had a circumference of 6.014 em and a ratio 
of R/r equal to 53.5. The instrument was cali- 
brated against known standard weights and also 


against conductivity water. 
2. Electrical Conductance 


The electrical conductance was measured with 
a dip-type conductivity cell, with a cell constant 
of 0.8307, and a conductivity bridge manufac- 
tured by the Solar Manufacturing Corporation, 
New York (model-RC1; 20 watts; 100 to 125 
volts; 50 to 60 cveles). The null point was 
detected with the help of a magie-eye, which 
flickered when the switch multiplier was turned 
from one position to another, indicating that the 
null point occurred in that region. The multi- 
plier was then kept steady and the resistance dial 
moved until the green portion of the magic-eve 
receded to a maximum. This was the null point. 
In the absence of accurate knowledge of the exact 
equivalent weights of the commercial soaps and 
soapless detergents, values for the specific con- 
ductance rather than the equivalent conductivity 
have been calculated and are presented in this 
paper. However, the data in figure 2 were cal- 
culated by the equation 
: 1,000“ K 
weight conductance 7 
where A equals specific conductance, and C equals 
weight of soap or soapless detergent in grams per 
liter. This was done to obtain curves similar to 
those that would be obtained if equivalent con- 
ductance were plotted. 


3. pH 


The pH of the various soap and soapless deter- 
gent solutions was measured with a glass-electrode 
pH meter. The instrument was standardized 
against standard solutions of Borax (pH equals 
9.18 at 25° C for a 0.01-molar solution). 
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4. Opacity 


The opacity of the various solutions was meas- 
ured with a Klett-Summerson photoelectric color- 
imeter. A projection type of electric lamp (100 
watts; 110 volts) fitted with a green filter (540 
millimicrons) was used as the source of light, and 
the current produced by the photocells was 
measured by means of a suspension wire type of 
galvanometer. The opacity of the system, that is 
the reciprocal of the intensity of transmitted light, 
was measured in terms of the divisions on the 
potentiometer slide-wire that are necessary to 
nullify the potential applied to the galvanometer. 


5. Foaming 


The rate of growth of foam was determined by 
measuring the height to which the foam rises under 
a standard bubbling rate. The bubbler was pre- 
pared as follows. A glass tube, }\ in. in diameter 
and 18 in. long, was sealed at one end with small 
pieces of copper wire (0.01 mm diameter) fused 
into it. When the end of the tube had cooled it 
was dipped into nitric acid, which dissolved the 
copper wire and left 20 fine holes of uniform 
diameter. The open end of the tube was then con- 
nected to an air-pressure valve via a flow regulator, 
and the bubbler was dipped into a 1,000-ml meas- 
uring cylinder containing 200 ml of the solution to 
be tested. The time was noted on a stop watch. 


IV. Results and Discussion 


Data for aqueous solutions of some typical soaps 
and soapless detergents are presented in figures l 
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Figure 2. Variation of specific conductance with con- 


centration. 


T3 toilet soap; @ C1 coco soap; A G2 glycerin soap; 4 W5 washing soap; 
85 shaving soap; §§ D5 soapless detergent. 


to 4, in which the surface tension, weight con- 
ductance, pH, and opacity of the various solutions 
are plotted with respect to concentration. Data 
for all the soaps and soapless detergent solutions 
at a concentration of 2.5 percent are presented in 
tables 1 and 2, since such a tabulation permits 
quick comparison and evaluation. 
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Figure 3. Variation of pH with concentration 


T3 toilet soap; @C1 coco soap; A G2 glycerin soap; & W5 washing soap, 
85 shaving soap; §§ D5 soapless detergent 


1. Surface Tension 


The curves obtained by plotting surface tension 
with respect to concentration are very similar for 
the various toilet soaps but differ somewhat from 
those for the coco, glycerin, washing, and shaving 
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TARBLe 1 Some physical-chemical properties of 2.5 percent 


aqueous soap solutions 
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soaps. Solutions of the toilet soaps are charac- 


terized by a very low value of surface tension at 
small concentrations (about 0.1 to 0.2 %), which 
increases to a Maximum value around 1.25 percent 
of soap content. The surface tension then de- 
creases slowly to a minimum value, around 2.5 
percent of soap, which either remains constant or 
increases very slowly as the soap content is in- 
creased. The curves thus show two minima, 
around 0.1 to 0.2 percent and 2.5 percent of soap. 

The solutions of coco, glycerin, washing, and 
shaving soaps are characterized by the surface 
tension passing through a minimum value in very 
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dilute solutions (about 0.2 to 0.3 % of soap), 
and then following a somewhat irregular and in- 
definite path up to about 1.25 percent. The 
value of surface tension then increases somewhat 
to reach a maximum, which is either retained or 
decreases slightly as the concentration of the soap 
is increased. In contrast to the data for solutions 
of the toilet soaps, solutions of this group of soaps 
do not show a second minimum around 2.5-percent 
concentration, 

In either case, the minima in the region of lowe 
concentration are due to the surface film being 
saturated with a single layer of soap molecules, 
whereas the minima in the region of higher con- 
centration indicate the formation of soap micelles 
with single or multiple charges. The steady value 
of surface tension obtained beyond a concentra- 
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tion of about 2.5 to 4 percent of soap indicates 
that both the surface and the interior of the system 
are saturated with respect to the charged micelles. 
The exhibition of two regions of minimum surface 


tension for solutions of the toilet soaps and of only 
one such region for solutions of other types of 
soaps is probably to be explained as being caused 
hy the presence of many sizes and kinds of colloidal 
particles in the former and of a constant particle 
size in the latter. 

The solutions of the soapless detergents are 
characterized by a continuously decreasing value 
of surface tension with increasing concentration, 
an almost constant value being reached after 
about 1.25 percent. The shape and nature of the 
curves for the solutions of the soapless detergents 
are thus very different from those for the soap 
solutions. It seems reasonable to conclude from 
these data that the concentrations at which mono- 
layers and ionic micelles are formed are much 
lower for the soapless detergent solutions than 
those for the soaps. In terms of washing action, 
this may be interpreted as indicating that at lower 
concentrations the soapless detergents are better 
surface-active cleansing agents than the soaps 


2. Electrical Conductance 


In the absence of data regarding exact equiva- 
lent weights of the soaps and soapless detergents, 
values of the specific conductance rather than of 
equivalent or molar conducitivity have been 
presented in tables 1 and 2, and values of the 
weight conductance were plotted in figure 2 (see 
p. 363). The curves in figure 2 show that the 
conductance for the various soap solutions rises 
rather sharply in dilute solutions to a maximum 
around 0.25 pereent and then decreases to about 
2.5 percent, bevond which an almost constant 
value is obtained. The decrease in conductivity 
is due to the replacement of ions by colloidal 
particles which, although conducting, have a lower 
equivalent conductivity than the ions from which 
they are formed. The variations in the regions of 
maxima for the different soap solutions is probably 
caused by the presence of many sizes and kinds of 
colloidal particles in equilibrium with one another 
and to the shifting of the equilibrium position 
with concentration and soap type. 

Whereas soap solutions exhibit a high electrical 
conductance in dilute solutions and low and de- 
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creasing values with increasing concentration, 
solutions of the soapless detergents are character- 
ized by very high values of conductance, which 
decrease rapidly up to 0.5 percent, then increase up 
to about 2 percent and finally decrease slowly with 
increasing concentration. The formation of col- 
loidal particles thus appears to take place at very 
low concentrations in solutions of the soapless 
detergents, and the relatively high charge carried 
by them tends to make them better cleansing 
agents than commercial soaps in solutions of 
similar low concentrations. 


3. pH 

The pH-concentration curves in figure 3 show 
that the values of pH for the soap solutions in- 
crease very sharply at low concentrations up to a 
maximum around 0.2 percent, followed by an 
equally sharp decrease and a minimum between 
1 to 2.5 percent. A further increase in the con- 
centration of the soap does not appreciably change 
the pH value. 

Table 1 shows that the values of pH average 
about 10 for solutions of different soaps at a con- 
centration of 2.5 percent. The presence of excess 
alkali in aqueous soap solutions favors the forma- 
tion of ionic micelles [6] and helps to clarify the 
solutions by driving back hydrolysis and simul- 
taneously increases the viscosity of the system. 
The relatively high alkalinity of the various soap 
solutions is caused by their readiness to hydrolyze 
in water and also by the free (unreacted) alkali 
in the commercial products. In either case, the 
alkali has an irritating effect upon the skin. 

The solutions of the soapless detergents may be 
either acidic or alkaline, depending upon the nature 
of the major constituent of the products. As with 
the alkaline soap solutions, irritation to the skin 
may be caused by these alkaline or acidic solutions, 
and likewise solutions of these soapless detergents 
would seem to be best used in greater dilution. 


4. Opacity 


The opacity-concentration curves for the various 
soap-water systems are shown in figure 4 and 
indicate that the solutions at lower concentra- 
tions are clear and transparent. With increasing 
concentration the solutions become at first cloudy 
and translucent, and then dense and opaque at 
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higher concentrations. It is presumptuous to 
attribute the turbidity of the soap solutions to 
their colloidal nature alone, because perfectly 
transparent solutions of various colloidal materials 
and translucent to opaque systems of crystalloidal 
substances are obtainable. Also, the insoluble 
acid soap may sediment out almost entirely on 
standing, leaving the bulk of the soap in a clear 
transparent solution, or it may remain partly 
or wholly suspended in particles whose dimensions 
range all the way from coarse suspensions down 
to the smallest resolvable particles. 

In concentrated as well as dilute solutions, the 
opacity of the coco and glycerin soaps is substanti- 
ally low. The former made from coconut oil, 
have a high laurate-myristate content; the opacity 
of solutions of soaps of the lower fatty acids is 
lower than that for corresponding solutions of 
The high trans- 
pareney of the glycerin soap solutions probably 
results from the presence of glycerin and sugars 
that are adsorbed on the surface of the soap 
suitably index of 


soaps of the higher fatty acids. 


particles and change the 
refraction. 

Solutions of the soapless detergents exhibit a 
very high transparency at both low and high 
concentrations. The initial low values of opacity 
occur in regions of concentration in which there 
is low surface tension and relatively high electrical 
conductance, and is probably due to the forma- 
tion of monolayers and primary aggregates of the 
detergent) molecules. The © slight 
opacity that starts around 1.25-pereent concen- 
tration corresponds to the region in which the 


increase in 


surface tension reaches a constant value, and the 
weight conductance levels off and then decreases. 
This is probably to be ascribed to the formation 
of secondary aggregates and ionic micelles. 


5. Foaming 


Measurements of the height to which the foam 
rises when a standard bubbling rate is employed 
indicate that there is no appreciable change in 
the rate of growth of foam with variations in the 
concentration of a given soap or soapless deter- 
gent. However, the rate is’ probably — slightly 
slower in dilute solutions and somewhat faster in 
concentrated solutions. 

The first signs of foam formation are shown in 
solutions containing approximately 0.0035 percent 
of soap, bui at this extremely low concentration, 
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the foam is very unstable and breaks down as 
soon as it is formed. Relatively stable foam is 
produced beyond about 0.007 percent. The for- 
mation of foam in these extremely dilute soap 
solutions shows that the stabilizing action of ad- 
sorption layers attains its maximum at surface 
concentrations much below the saturation of the 
surface layer (about 0.1 to 0.2%). This stabiliza- 
tion arises not only by a mere enhancement of the 
mechanical properties of the surface layer, but 
also by a certain degree of mobility of the stabiliz- 
ing molecules necessary for the rapid bridging of 
ruptures in the surface layers. 

The foam produced from solutions containing 
0.078 to 10 percent of soap or of soapless detergent 
is very stable and remains almost unchanged for 
an hour or more if the system is not disturbed by 
mechanical or thermal agitation. The require- 
ments for a stable foam, enabling the bubbles to 
withstand mechanical shocks, are a lowering of 
the boundary tension by adsorption and suitable 
mechanical properties of the adsorbed stabilizing 
film. It is worth while remembering, however, 
that the power of foaming is not always or entirely 
For in- 
stance, saponin solutions, whose surface tensions 
are only slightly less than that for water, bave a 
very great power of foaming on account of the 
formation of a solid or highly viscous film of an 
adsorbed material on the surface. It should be 
very worth while to examine the decay of foam 
under standard conditions of agitation and the 
stability of single bubbles liberated beneath the 
liquid surface, in the ease of both the commercial 


to be ascribed to low surface tension. 


soap and soapless detergent solutions. 


6. Correlation with Data on the Phase Nature and 
Structure of the Solid Soaps 


X-ray diffraction studies of commercial soaps 
have shown [2] that the toilet and “‘medicated” 
soaps consist essentially of beta sodium palmitate 
and small amounts of omega sodium oleate, the 
coco soaps consist primarily of cmega sodium 
laurate and sodium myristate, the shaving soaps 
consist essentially of two or more distinct phases 
(alpha, and/or beta, omega, and delta) of sodium 
palmitate, and the constituents of the glycerin and 
washing soaps are rather indefinite, the latter being 
probably in the omega phase. Electron micro- 
scopic investigations [3] also indicate that sodium 
palmitate may be in the washing soaps. 
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With a view to obtaining some information 
regarding the general physical-chemical properties 
of the solid soaps and their phase nature, some 
qualitative experiments were conducted on the 
hardness of the soaps, their solubility, and their 
reaction to water. In general, it was found that: 

The hardness of the soaps, based on the resist- 
ance offered to the fine edge of a steel spatula on 
application of a 100-g load, decreases in the 
following order: Toilet, ‘“medicated’’, washing, 
shaving. 

The solubility of the coco, glycerin, and ‘“medi- 
cated”? soaps in warm water is the greatest of the 
groups and decreases in the order: coco, glycerin, 
medicated. The shaving soaps are more soluble 
than the toilet soaps and give transparent to 
translucent solutions on cooling, whereas the solu- 
tions of the toilet soaps set to opaque gels or curds. 
The solubility of all the five types of soaps is, 
however, relatively low in cold water and in hard 
water, 

The percentage of soap rubbed off the cake 
(bar), when used in water against a laboratory 
towel, decreases in the following order: Shaving, 
toilet, washing. 

On standing with water, the toilet and shaving 
soaps swell and disintegrate, whereas the washing 
soaps crack somewhat with little or no swelling. 

Data in tables 1 and 2 show that the values of 
most physical-chemical properties are distinctly 
different for solutions of soaps as against those of 
the soapless detergents. Aside from this, the most 
striking behavier is perhaps the ability of the 
soapless detergents to dissolve rapidly and work 
in cold as well as in hard water. A considerable 
part of the soap used in household applications is 
probably employed to soften hard water in order 
to make the additional soap effective. On the 
contrary, the soapless detergents usually work 
equally efficiently in hard and in soft water, none 
being wasted in the task of softening the water. 


7. Micelles and Their Applications in Cleansing 
Action 


There is some difference of opinion as to the 
origin, manner of formation, structure, size, ete. 
of the colloidal aggregates termed loosely as 
“micelles”. McBain [7] postulates at least two 
types of micelles (ionic and lamellar), whereas 
Hartley [1] postulates only one, namely, the ionic 
micelle. Solutions of different colloidal electro- 
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lytes exhibit different phenomena with varying 
concentrations, and this behavior can be explained 
only if one assumes the existence of more than 
one type of micelle. It is the general belief that 
the ionic micelle in soap solutions, which is of 
significance in the present investigations, con- 
sists of a loose spherical aggregate of the paraffin 
chain ions with the polar groups oriented toward 
the aqueous phase. Recent investigations by 
Mattoon, Stearns, and Harkins [5] have shown 
that aqueous solutions of soaps and soapless de- 
tergents give a Bragg spacing that is independent 
of concentration and that is equal approximately 
to the double length of the molecule. They have 
designated this as the micelle thickness, their con- 
ception of a micelle being cylindrical, rather than 
spherical or lamellar. 

The evidence for the existence of ionie micelles 
in solutions of commercial soaps and soapless 
detergents is based primarily on data on electrical 
conductance and surface tension. The former 
show the formation of micelles by the rapid de- 
crease in conductivity with increasing concentra- 
tion—the ionic micelles have a lower conductivity 
than the ions. Surface-tension data indicate the 
formation of monolayers of soap molecules on the 
surface of the solutions at very high dilutions and 
the formation of colloidal aggregates (ionic 
micelles) at slightly higher concentrations. Elee- 
trical conductance and surface-tension data for 
the solutions of the soapless detergents indicate 
that micelle formation reaches a very high value 
at relatively high dilutions. 

The mechanism by which solutions of soaps 
and soapless detergents dissolve relatively insol- 
uble substances, such as grease and dirt, can be 
explained simply in terms of either concept of the 
ionic micelle and the action of mixed monolayers: 
(a) the hydrocarbons in the dirt are dissolved in 
the interior of the micelle; (b) the polar com- 
pounds in the dirt are removed by adsorption on 
the surface of the micelle, with the hydrocarbon 
part inside the micelle and the polar groups in the 
aqueous phase; and (c) compounds that are soluble 
in water but insoluble in hydrocarbons are taken 
up by the polar group of the micelle and thus act 
as peptizing agents. 

Since the degree of solubilization is related to 
the number and size of the micelles present [S}, it 
seems reasonable to conclude that greater solu- 
bilization and hence more efficient cleansing may 
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be achieved by using soaps in concentrations 
ranging from 2.5 to 4 percent and soapless deter- 
gents in less concentrated solutions (about 1.25% 


or less 
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Doubly Charged Ion Spectra in Mass Spectra 


of Hydrocarbons 


By Fred L. Mohler, Evelyn G. Bloom, Edmund J. Wells, Jr., Jonathan H. Lengel,’ and 


C. Edward Wise 


A compilation of doubly charged ions in the mass spectra of 148 hydrocarbons has been 
made, based on mass spectra published in the American Petroleum Institute Catalog of Mass 
Spectral Data. In hydrocarbons with three, four, or five carbon atoms, the most probable 
double ionization process involves loss of all but two or three hydrogen atoms without 
breaking carbon bonds. Unsaturated molecules give larger doubly charged ion peaks than 
saturated molecules, and ions with an even number of hydrogen atoms tend to be more 
abundant than those with an odd number. In saturated hydrocarbons with six or more 
carbon atoms, double ionization with breaking of carbon bonds is most probable, but in 


benzenes and highly unsaturated molecules, double ionization without loss of carbon atoms 














is probable. These heavier molecules show distinetive differences depending on moleec- 


ular structure. 


The ratio of ion current at 70 volts to that at 50 volts ionizing voltage is an indication 


of the appearance potential of the ion. It increases as the number of hydrogen atoms 


removed increases, and the results indicate that hydrogen atoms are removed in pairs 


to give Ho. 


I. Introduction 


The mass spectra of hydrocarbons and other 
polyatomic molecules in general show some small 
peaks that arise from doubly charged ions. Ions 
with a double charge are characterized by an ap- 
parent mass equal to half the molecular weight of 
the ion. If the molecular weight is an odd num- 
ber, the doubly charged ion will be distinguished 
by half-integer values; if it is an even number, 
then it will often happen that the apparent mass 
of the doubly charged ion is coincident with the 
mass peak of a singly charged ion. For this 
reason the data on the occurrence of doubly 
charged ions are necessarily incomplete. 

The application of the mass spectrometer to 
chemical analysis has resulted in the careful 
measurement of mass spectra of at least 200 hy- 
drocarbons; and these spectra afford a vast 
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amount of data on the existence of doubly charged 
ions, but no attempt has previously been made to 
collect and correlate these data. The Catalog 
of Mass Spectral Data of the American Petroleum 
Institute [1] * has been the source of data used in 
this compilation, and the tables contributed by 
the Mass Spectrometry Section of the National 
Bureau of Standards have been used except in 
cases noted in the text. There are distinet ad- 
vantages in using data of one laboratory insofar 
as possible, for there are minor differences between 
laboratories arising from instrumental differences 
and different methods of operation. In three 
cases, ethane [2], ethylene [3], and benzene [4], 
appearance potentials of some doubly charged 
ions have been published. For most molecules 
all that is known is the relative intensity of the 
doubly charged ions at one or two ionizing vol- 


nome 
tages. 





? Figures in brackets indicate the literature references at the end of this 


paper. 
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II. Experimental Procedure 


Mass spectra have been measured with a 180° 
Consolidated mass spectrometer following recom- 
mended [5]. The instrument is 
equipped with automatic electron current control, 
and the ionization chamber is thermostatically 
maintained at 245° C. The compounds used are 
in most cases the pure NBS Standard Samples. 

The Catalog of Mass Spectral Data [1] gives 
intensities of mass peaks relative to the maximum 
peak taken as 100. The limit of this instrument 
is about 0.01 percent of the maximum deflection, 
and as spectra are always obtained with a pres- 
sure sufficient to give a nearly full seale deflection 
at the maximum, the lower limit remains about 
the same in all spectra. Peaks of doubly charged 
ions rarely exceed 1 percent of the maximum, and 
most of the peaks are less than 0.1 percent. The 
reproducibility on these small peaks is about 
10 percent of the peak height, or 0.01 percent of the 
maximum for peaks less than 0.1 percent. 

The tables in the catalog give relative intensities 
with ionizing voltages of both 50 and 70 volts. 
This paper tabulates values at 70 volts and in- 
cludes some data on the ratio of ion current at 
70 volts to that at 50 volts. The relative intensity 
at 70 volts is greater than at 50 volts and often 
much greater. The APL Catalog gives spectra 
from several laboratories for many of the com- 
pounds, and in general there is good agreement on 
these small doubly charged ion peaks in all cases 
where they are reported. 


III. Results 
1. C, Hydrocarbons and Methane 


procedures 


No doubly charged ions in the mass range 6 to 8 
are observed in methane. Table 1 lists the doubly 
charged ions in CH;.CH;, CH,=CH, and CH= 
CH. The doubly charged ions all contain two 
carbon atoms, and those ions of even molecular 
weight are coincident with the singly charged ions 
CH*, CH,* and CH,*. However, in two cases 
there is a basis for estimating the intensity of 
doubly charged ions of even molecular weight. 
Kusch, Hustrulid, and Tate [3] report a small peak 
in ethylene at mass 14'5 from the ion C’C’H,**. 
The CYH,** peak will be about 50 times as large 
as the 144) peak. They used a mass spectrometer 
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of greater sensitivity than the Consolidated in 
strument. Similarly, in acetylene a peak has bee: 
observed at 13% from the ion C?C™H,**, and th: 
C?H,** peak is computed from this. The dis 
tinctive differences between the three C,) com 
pounds are of a type not found in the heavier 
hydrocarbons. 

TaBLe 1. Doubly charged ions of Cz hydrocarbons 


Relative intensities in percentage of maximum peak in spectrum 


Compound CoH C He CoH CoH Cals 
Ethane 0.07 0.70 
Ethylene 10. 007 42 0.3 
Acetylene 73.9 


Kusch, Hustrulid, and Tate, Phys. Rev. $2, 843 (1937 The value f 


(>)Hy** is an estimate based on a small peak at 144% from CUCUHg*+ 
* The value for CgH ** is based on a peak at 1344 from CCH )*+ 


2. C, Hydrocarbons 


Table 2 lists doubly charged ions found in the 
five hydrocarbons CH;.CH,).CH,, CH,).CH»).CH), 
CH,=CH—CH;, CH,=C=CH:, and CH, 
C=CH. All observed doubly charged ions con- 
tain three carbon atoms, and the apparent masses 
fall in a mass range 19 to 21 containing no singly 
charged ions. Delfosse and Bleakney [6] have 
previously published data on doubly charged ions 
in propane, propylene, and propadiene at 100-volts 
ionizing voltage. 

TABLE 2. Doubly charged ions of C3 hydrocarbons 


Relative intensities in percentage of maximum peak in spectrum 


Compound CaHye CyHy CyHy CyHs CiHe 
Propane 1.10 0.61 1.12 0. 26 0.02 
C yclopropane 3.33 1.44 2.40 73 02 
Propylene 3.23 1% 2.35 67 4 
Propadiene ! 2.79 3. 20 3.61 
Propyne 3.23 2.34 0.80 


From unpublished mass spectra at a %)-volt ionizing potential. 


With one exception, the relative intensity of 
the ions C,H,**, C,H,** and C,H,** are similar 


in the five molecules in contrast to some distinctive 
differences in the spectra of the singly charged ions. 
Propyne differs from the others in that the C;H,** 
peak is small, and this is of interest because other- 
wise the mass spectra of propyne and propadiene 
are nearly identical. 
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3. C, Hydrocarbons 


Table 3 lists the doubly charged ions in ten C, 
hydrocarbons. C.H,* ions of 
mask doubly charged C, ions of even molecular 
weight, but in 1,3-butadiene a small peak at 27% 
from C'C?H,** permits a rough estimate of the 
intensity of the C?7H,*t peak. 


mass 25 to 29 


Tasie 3. Doubly charged ions of Cy hydrocarbons 


Relative intensities in percentage of maximum peak in spectrum 


Compound CsHy Cally CyoHy Call Cay Cas Cie) Cay 
Butanm 0.04 ool 0. 02 0.36 0. OS 0.05 
Butane 05 Led (4 21 o2 

1-Butene On 03 02 1. 06 17 06 
s-2-Butene 02 01 01 1.55 27 17 
rna-2-Butene 1. 6 29 27 

i-Butene On 03 038 0.95 09 

1,2-Butadiene 03 02 1.07 cD) 

1,3- Butadiene ay 02 Low ‘a LO 

1-Butyne 06 os os 3. 42 a 

2-Butvne Ol $. 25 46) 


Estimated from a smal! peak of CC 3"2H¢*+ 
The relative intensities are roughly similar in 
all these spectra, and there is no evident relation 


hetween the intensity distribution in the mass 


TABLE 4, 


spectra of doubly charged and singly charged 
ions. C,H,;** is in every case the largest peak, 
and the height is least in C,H isomers, greater 
in the C,H, isomers, and largest in C,H, molecules. 
There is a small but real difference between the 
cis- and trans-2-butene (the singly charged ion 
spectra are very similar). 


4. C; Hydrocarbons 


Table 4 lists doubly charged ions of 12 C; 
hydrocarbons. It includes all C;sHy» and C;Hio 
isomers except some rare cyclics but only two 
of the many C;H, isomers. As in the case of 
C, and C; hydrocarbons, the relative intensities 
in all these spectra are similar. C;H,** is the 
largest doubly charged peak and it is smallest 
in C;Hy isomers, larger in CsHy isomers, and 
largest in the C;H, isomers. Ions with three 
or four carbon atoms give small peaks often near 
the threshold for observation. lons with an even 
number of hydrogen atoms tend to be more 
abundant than those with an odd number. This 
is probably a fairly general rule for doubly charged 
ions, partially concealed by the fact that there is 
often masking of ions with a double charge and 
an even number of H atoms. 


Doubly charged ions of Cs hydrocarbons 


Relative intensities in percent of maximum peak in spectrum 


Compound CsHe Css Cala  CyHy 
n-Pentane 0.04 
1-Pentane . 06 
Neopentane 4 
Cyclopentane 
Pentene 0.02 0.01 0.01 03 
cis-2-Pentene . 06 
trans-2-Pentene ol ol . 06 
2-me-I-Butene 06 
}-me-1-Butene . .10 
2-me-2-Butene ol Ol .4 
Cyclopentene 03 02 ol .02 
2-me-1,3-Butadiene ? 03 01 ol 15 


These also have a small peak at C4Hs*+. 


5. C, Hydrocarbons 


Table 5 lists 23 C, hydrocarbons including all 
C,H, isomers and most of the common C,H, 
isomers but only two C,H, molecules. Doubly 
charged Cy ions with an even number of hydrogen 
‘toms are masked by singly charged ions. These 
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CysH CsHe CsHs CsHy CsHs CsHe CsH;, CsHe (gH, 
0. 33 0. O08 0.13 0.05 0.07 0.02 
~27 05 .09 . 02 .08 
to) 02 04 oO 
0.01 | .19 44 .16 24 .05 0.01 
ol 1.038 24 4 -2i 22 .10 4 0.02 
08 1.23 37 53 >) .22 a OS 02 
Ol 1. 52 37 BB . 26 . 24 .10 ay 03 
01 oO .21 20 15 18 06 4 ol 
Ol iA 14 15 Os on 04 
.02 1.15 25 .39 17 29 07 12 .02 
On 2.43 4 1.00 .70 74 3 13 
10 4.31 1.07 1.07 sO . 85 uw) 41 


2 Also has a peak Cg**0.08, 


to C 
show considerable variations and do not 
the rule that double ionization without loss of a 
C, doubly charged 


spectra in contrast C,, and C, molecules 


Se 


follow 


carbon atom is most probable. 
ions are missing in C,H, isomers and are relatively 


weak in the olefins. The intensity distribution 
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TABLE 5. Doubly charged ions of Cy hydrocarbons 


Relative intensities in percentage of maximum peak in spectrum 





Compound C Hy CsHy CsH; CsHy CsHs CsHe¢ CsH- CsH¢ CsH» CH CoHs CeHy CoH, 
n-Hexane 0.02 0.07 0.02 0.06 0.03 0.05 0.02 0.01 
2-me-Pentane 02 ON 02 05 93 On 03 02 0.02 
}-me-Pentane 02 ox 03 06 03 03 
2,2-me>- Butane 02 03 03 07 02 03 03 03 ol : 
2,3-mez- Butane 01 06 01 04 . 02 02 02 02 j 
C yclohexane 03 06 O68 02 O4 0.01 0.07 0.03 
me-C yclopentane .12 02 OS OH 07 03 n 02 1 
1-Hexene ol 15 O44 Tr) O5 03 02 on | 
cis-2- Hexene 02 nw O4 12 02 02 ol O05 02 3 
trans-2-Hexene 02 14 4 10 02 02 ol O5 Ol ] 
trans-3-Hexene 02 28 09 14 1! 24 15 02 412 09 ol 4 03 
}-me-1l- Pentene 02 18 O5 10 O07 Oo 02 
4-me-1-Pentene 17 O5 Is 13 17 13 03 07 ol 
2-me-2- Pentene 02 19 O5 09 05 03 14 0s O05 05 
}-me-cis-2- Pentene ol 14 4 12 04 03 ol 02 4 02 
3-me-trans-2-Pentene ol 19 4 10 O4 03 02 03 O5 03 
#-me-cis-2-Pentene 32 oo 16 il 0 07 02 02 ol ; 
4+ me-trans-2-Pentene ny as 17 12 w O7 2 02 ol 4 
2-eth-1- Butene 4 3 O05 10 ie ie 03 02 oo O38 
},3-me>- l- Butene 02 18 On} o7 O68 02 of 
2,3-me>-2- Butene 17 O05 O7 02 02 Ol ol ON 
Cyclohexene _ . 03 11 02 05 02 25 il 53 oo 
1,5 Hexadiene 03 “4 04 O5 03 O7 o O5 


This also has peaks CaHy**0.03, CyHy**0.02, CsH4**0.02, CsHs**0.01, 


TABLE 6. Doubly charged ions of C; hydrocarbons 


Relative intensities in percentage of maximum peak in spectrum 


Compound CsH CsH CsHy CsHs Coste CsHy CsHe CsHy CsH C.H CoH, CoH Cr;H, CrHs CrHy CrH 
n-Heptame 0.03 0.02 0.02 0. 02 
2-mne-Hexane v2 02 02 0.02 01 
}-me-Hexane 03 Ts] 02 Ol 
teth-Pentane “2 Ol 
2,2-me?-Pentane a2 oo “3 0.01 2 07 
2,3-me;-Pentane 2 4 ol 
2,4-mes-Pentane Os a2 ol og 0.05 0.05 0.09 0.06 02 0 
3,3-me-Pentane ss 0 02 
2,2,3-mes- Butane v2 03 02 0 
me-C yclohexane Os 03 4 4 03 0.05 0.20 0.06 % 0.02 ON 
eth-C yclopentan ss ol (2 2 ol ol v2 16 5 19 02 2 
1,1-mes-C yclopentane . 4 “2 ol ol Is 12 4 0 19 ol 02 
cis-1,2-me2-C yclopentane 07 0 03 Ol 0 Ta] 02 16 O5 21 08 02 
trans-1,2-me>-C yelopentane i O5 v2 02 02 i4 Oo4 2 ‘] 02 
cig-1,3-me;-C yclopentane 7 re | 038 OS 07 Os 13 17 M4 21 “2 m4 
trang-1,3-mes-C yclopentane 7 4 4 Tt) il 07 12 16 0 *% 03 4 
1-Heptene Ww 07 2 Os 
4,4-me>-1-Pentene 4 02 08 Ww 1" 02 02 
2,3,3-mes-1-Butene ? OF (2 11 ty 05 (4 
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ong the C; ions is frequently much like that in 

. C; hydrocarbons, though the intensity is less. 

ins-3-hexene is an exception, for here C;H,** 

the largest peak, and this peak is small or absent 

all other Cy hydrocarbons. The benzene 
spectrum is given in table 9. 


6. C; Hydrocarbons 


Table 6 lists doubly charged ions of 19 C, 
ivdrocarbons. lons C;H,** and C;H,** would 
be masked by singly charged ions in C;Hj,, and 
(.H,, isomers. The nine C;H,, isomers give a 
variety of mass spectra of doubly charged ions 


TABLE 7. 





with Cs; ions the most prominent, except in one 
5 
case, 2,2-dimethylpentane, which has a maximum 


peak at C,H,**. The only other C; with a maxi- 
mum at this peak is 4,4-dimethyl-l-pentene. 
This is a structurally similar molecule, except that 
it contains a double bond. The cyclies all give 
quite similar doubly charged ion spectra with 
C,H,** the maximum peak. 
spectrum of methylbenzene. 


7. C, Hydrocarbons 


Table 7 lists doubly charged ion spectra of 
15 CH, isomers, 22 cyelic and olefin CH), 


Table 9 gives the 


Doubly charged ions of Cy hydrocarbons 
4 g A y 


Relative intensities in percentage of maximum peak in spectrum 


Compound CsH, CsHy CsHe CsHg C3H CoH; 
Octane 
2-me-Heptane ool ool 
me-Heptane 2 ol 
j-me-Heptane aw we 


me -Hexane 
»Hexane ol vl 
»Hexane ol 








2.4 
mep-Hexane 

mn Hexane oy “ 

1-me)-Hexane ol 0] 

2,3-mer Pentane 

2,4-mes-Pentane oO 0 ool 





}-mes-Pentane ol 
4-mey- Pentane 


j-me,y- Butane 


Octene 


s-4-Octene 


2.4,4-me-1-Pentene 2 
2,4.4-mey-2-Pentene Os 
‘ ohera 
1, 1-me ol 
8-1, 2-1 
rana-|,2-Ie 
s-1 4-1 Ti 
franas-1,4-me Ol 
s-1,.b-mMe ed 
frana-1,3-me wn 
ethyl as al oo 


( yclopentanes 


Propyl as (rR 
Propyl "n ol 
me-l-eth v2 us 
s-1-me-2-eth (2 
1.1, 2-1 nm LP 
}-I¢ mi) 2 02 O5 
cis, cis-1,2,3-m« ol mi ol 
&, frana, cia-1,2,3-me 2 "2 LD 
s, frans, cis-1,2,4-me OD ol ol ol 0.04 
cis, trans-1,2,4-me ol al " ol “4 
1-Vinyl-l-eyclohexene "4 


mm OO] O04 0.08 
a ol 
i 
O08 ws i4 4 
(h 
ool (4 ol 
1 ow ts a 1) 
re 
in 
“ ira 
4 WF Ks 4 ag oe oOo 
hy 
ir iM 
tr 1S) Led 
hy 
Wt ik 
Oo; iM 21 iM Os oo ol ww) 
0.02 a» 1! tis 17 lw ol 2 
yA} 07 rY 1] Os it “i (4 
21 0 2 4 (4 ol 
AY 0. A (rR in 
i“ 21 i | On tt 
(2 ol i7 Os 22 4 4 ol ol 
ol 4 Os 14 ns ny i 
ol 2) Ww y ry iM) as 
Ww aw ms aK as 
11 (4 1 a] 
Pall ww uM) 7 iM ol 
1 0 1 ta al 02 
l 7) lt ww i“ 
22 in a os m4 ths 
(4 Lid wv 15 % in Ww iw iM 24 
12 0. lt wm m 
Vt 0 22 (i a 
12 4 i4 (1 
iw On l i“ ol 
0 2. 20 11 Ww 


Also gives peaks C;H)++0.18, C;H)++0.02, C;H4+*0.07, CgH5*+0.04, CsH5**0.01, CgH;**0.09, CoHy**0.20, Cel) **0.01, 


Doubly Charged Ions in Mass Spectra 


373 





a: =) 7 


— 


be Bs 





TaBLe 8. Doubly charged ions of Cy hydrocarbons 


Relative intensities in percentage of maximum peak in spectrum 


Compound Crs CrHy 
n-Nonane 
2,2,5-mey- Hexane 0.04 
2,2,4-mey- Hexane as 
2,2,5-mey- Hexane 07 0.01 
2,5, 3-mes- Hexane 


2,3, 5-me)- Hexane 
2.4.4-me>- Hexane 
3, 4-mnes- Hexane ol 
-ethe-Pentane 
2,2,3,3-me,y- Pentane 
2,2,3,4-me,-Pentane i ol 
2,2,4,4-me;- Pentane Te ol 


 3,3.4-mee Pentam 


( peer ‘ 
Butyl Os 
i-Butyl! 13 0s 
Cyelohexa 
n-Propyl 0 
#-Propyl mn Os 
1,1,3-mne 4 ow 
1-Nonene iM 
Phis also has a peak of 0.02 at CsHy* 


isomers, and 1 C,H, molecule. Cy ions with an 
even number of hydrogen atoms would be masked 
by singly charged ions. Three C,H, isomers not 
included in the table give no observable doubly 
charged ions. These are the three octanes that 
have ethyl side chains; 3-ethylhexane, 2-methyl-3- 
ethylpentane and 3-methyl-3-ethylpentane. Most 
of the CH. molecules give small doubly charged 
peaks, and in 12 cases the maximum peak is 
C-H.**. The largest peaks and the greatest 
number of peaks are in 2,2,4-trimethylpentane. 
Among the C,H, cyelies, the most abundant 
doubly charged ions are in 1,1,3-trimethyleyelo- 
pentane and of the four olefins 2,4,4-trimethyl-2- 
pentene has the largest peaks. These three com- 
pounds have a structural similarity that is parti- 
ally concealed by the rules used for numbering 
the side chains. The C,H. molecule, 4-ethenyl-1- 
cyclohexene, shows many peaks not found in any 
other molecules in this table. The occurrence of 
C, ions ranging from C,H,** to C.H,,** is signifi- 
eant. The ions C;H,.** and C;H,** would be 
masked in most of the other compounds in table 7. 
The spectrum resembles somewhat the alkyl 
benzenes of table 9. 
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8. C, Hydrocarbons 


Table 8 lists doubly charged ions of 13 C,H] 
isomers and six C,H,, isomers (cyelics and on 
olefin). One nonane, 3,3-diethylpentane, shows 
no doubly charged ions. In the octanes it was 
also found that isomers with ethyl side chains giv: 
no doubly charged peaks. With one exception th 
doubly charged ions of these molecules are C; ions 
The exception is 1,1,3-trimethyleyclohexan 
with two small peaks at C,H,** and C.H,,* 
This molecule gives the most intense doubl) 
charged spectrum of the five cyelics listed. It was 
found that 1,1,3-trimethyleyclopentane gave tli 
largest peaks of any of the C,H), cyelies listed i 
table 7. A following section gives evidence that 
peaks at 96, 97, and 98 are C;H,,.**, C;H,,**, an 
C;H,,**; and not C,**+, C,H**, and C,H,**. 

It is of interest that n-butyl-cyclopentane an 
n-propyleyclohexane have similar spectra and tha! 
the iso-butyl and jso-propyl are also similar 
each other and different from the other cyelies 

In general the doubly charged spectra of ( 
hydrocarbons are much like those of Cs hydro- 
carbons. C;H,*t* is the maximum peak in mos! 
cases with only two striking exceptions. 


Journal of Research 











ee a ee 


+ eR NE 


Nt. 





on 
ows 
Was 
rive 
the 
OTs 


an 
l 

bly 
Wiis 


the 


hat 


ana 


Th 


ch 











Taste 9. Doubly charged ions of benzene and alkyl benzenes 


Relative intensities in percentage of maximum peak in spectrum 


Compound C.H | CoH CsHs CeHs | CrH CrH2 | Cr 
Benzene 0.05 1.15 0.35 73.0 
Me- Benzene * -4 0.13 1.73 0. 
Eth- Benzene 4 02 10 
1,2-me;-benzene 4 15 
1,3-me; benzene 03 ol 16 
1,4-mer benzene 06 

Also has peaks at C,Hy** 0.08 and CsHy** 0.05 


Based on a peak of 0.20 at 39.5 from C!Cs"He* 
9. Benzene and Alkyl Benzenes 


Table 9 gives doubly charged ions of some aro- 
maties. These have been listed in a separate table 
for convenience of tabulation, because they include 
many ions not found in other hydrocarbons. The 
peak heights are larger than in other C,, C;, and 
(, hydrocarbons. 

Benzene has a very large C,H,** peak and only 
one small C,; peak in contrast to most of the 
spectra in table 5. Cyclohexene is the only com- 
pound in table 5 that has large Cy peaks. 

Methylbenzene qualitatively the 
C-H,, eyelies (table 6) except for differences 
definitely related to the fact that there are 8 H 
atoms instead of 14. C;H,** and C;H,** 
be masked by singly charged peaks in the cyeclics. 


resembles 


would 


The CH, aromatics give spectra quite unlike 
any of the Cy molecules of table 7 except the 
C.H,, molecule, 4-ethenyl-l-cyclohexene. It dif- 
fers from this in having Cy, ion peaks that are 
larger than the C; peaks. The estimate of CcH,,** 
is based on very small isotope peaks and is quite 
uncertain. The fact that C;H.** 
ing here in contrast to other Cy molecules un- 


is small or miss- 


doubtedly is a consequence of the fact that dis- 
sociation by removal of CH, gives C;H;**, whereas 
multiple dissociation or rearrangement is required 


to vive CH, i 


10. Hydrocarbons with 10 or More Carbon Atoms 


There is little published data on these heavier 
hydrocarbons. n-Decane, CyH», gives only two 
small peaks, C,;H,**0.02 and C,H,**0.02. 

1-Decene, CioHx», shows only one doubly charged 
ion, C7H.**0.038. 
octene, whereas I-nonene gave this ion and C;H,** 


This is also the only ion in 1- 


Doubly Charged Ions in Mass Spectra 


12 


CrHy CrHs CrHe | CH: | CHa |) CsH CsHs Cay) CsHe Coablio 


0.97 0.49 4.87 2.35 3.45 

05 O45 05 Pi) “2 0.05 0. OS 0.32 0.24 

O08 O5 O05 2S 0s On ll 73 6 40.2 
07 O5 4 21 ol 10 12 73 -- 42 
15 O4 02 22 02 07 61 52 “5 


3 Also has peaks at CsH )** 0.17, CsHy** 0.02, and CsHy** 0.11 
* Based on a peak at 53.5 from C™C;2Hyo**. 


The four butyleyclohexanes give doubly charged 
spectra like the C,H). cyclics, except that the 
peaks are smaller. All the ions have seven carbon 
atoms. 

The Atlantic Refining Co. has published mass 
spectra [7] of the three  diisopropylbenzenes 
(C\,Hjs) at an ionizing voltage of 50 volts. These 
include the heaviest doubly charged ions recorded 
to date. The three spectra are similar, and 1,3- 
diisopropylbenzene gives the following peaks: 
C,H,;**0.37, C,H,**0.08, C,H, **0.05, CyoH,* *0.75, 
CyoHy**1.25, CroHy**2.79, CroH3**0.53, and 
C,,H,;* 70.65. 


IV. Ratio of Peaks at 7O- and 50-Volt 
Ionizing Voltage 


The API tables gives mass spectra obtained 
with ionizing voltages of 50 and 70 volts and, in 
general, the doubly charged peaks are considerably 
The value of the 
ratio of the currents at 70 and 50 volts, 7(70) 7(50), 


larger at the higher voltage. 


is of interest, as it gives at least qualitative evi- 
dence as to the appearance potential of the ion. 

A plot of ion current as a function of ionizing 
voltage rises at first almost linearly from the 
appearance potential but gradually becomes con- 
cave to the axis of abscissas approaching a flat 
maximum value at roughly 30 volts above the 
threshold. Llonizing voltages of 50 and 70 volts 
fall near the flat maximum for the more probable 
single ionization processes, and the ratio 7(70)/i(50) 
is nearly unity. However, if the appearance 
potential is high, then the ratio will be much 
greater than unity and approach infinity as the 
appearance potential approaches 50 volts. 

There are a few direct measurements of appear- 
ance potentials of doubly charged ions and these 
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are correlated with the observed values of i(70)/- 
(50) in table 10. There is a rather large experi- 
mental uncertainty in the ratio amounting to at 
least 10 percent, and one can only conclude that a 
value of 1.2 to 1.3 corresponds to a potential 
between 30 and 36 volts and 1.4 to roughly 40 
volts. Above 40 volts the ratio will increase 


rapidly with increasing appearance potential. 


Tasie 10. Correlation of appearance potential and the 


ration of ton currents at 70 and 50 volts 


Compound Reference lon Appearance i(70)/i(50 
potential 


CoH CoH s2 1.2 
CoH, } C,H ue 

CoH i CH 7 1.0 
CH CoH 4 1.2 
CoH CyHy us 14 


Observed values of the ratio 7(70)/7(50) for the 
more abundant doubly charged ions in hydro- 
carbons with two to five carbon atoms and for 
benzene are listed in table 11. It is found that 
the ratios for the same ions in different isomers 
are nearly equal. In general, the range of values 
indieated is no greater than the experimental 
uncertainty. This is true even when the isomers 
are chemically different, as in the case of eyelies 
and olefins. The table is arranged to show that 
the ratio increases as the number of hydrogen 
atoms removed increases, and that ionization 
processes involving removal of equal numbers of 
hydrogen atoms from different molecules give the 
sume ratio within the range of experimental un- 
certainty. There is not a progressive decrease in 
the ratio as the number of H atoms removed 
becomes less, but the ratio remains nearly the 
same for removing 7H and 8H and similarly for 
5H and 6H, 3H, and 4H and for 1H and 2H. 

Interpretation of ratios in terms of appearance 
potentials.—The interpretation of this is fairly 
obvious. Studies of appearance potentials of 
singly charged ions show that frequently but not 
always, the hydrogen atoms are removed in pairs 
to form H,. Thus Delfosse and Bleakney [6] find 
that in propane, C,H,* is produced by the transi- 
tion C,H,-C,H,*+H,+H requiring 14 volts. 
The process C,H.—-C,H,*?+2H, oceurs at 14.7 
volts, and C,H.--C,H,* +2H,+-H requires 15.7 
volts. As it requires about three electron volts 
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Tas_e 11.—Ratio of peaks at 70- and 50-volt ionizing volta, 


lonization precess 170) /1(50 Remarks 

CsHp eo CsHy*+10H 15tots 3 isomers 

CysHy oC sHy*+8H 258 to 3.7 7 isomers, mean 3.5 
CyHy oC lly +7H 2.2 and 3.5 2 isomers 

CyHe oC Hy *+6H 2.0 Propane 

CsHy oC 3H +0H 2.2 and 2.45 2 isomers 

CsHyw oCsHyt+6H Litol9s 7 isomers, mean 1.7 
CyHe oC ly *+5H iY Propane 

C\He oC Hy? +5 LAS to LSS 4 isomers 

CyHy eC Hy *+5H 20 2 isomers 

CsHe oCsHy*+5H 17 and 1.9 2 isomers 

CyHy oCyHy* +4H 1.43 and 1.48 2 isomers 

CyHe oCyHyet +4 15 Propane 

CysHe eo CsH yet +4H 13 and 14 2 isomers 

CsHy 9 CsHett+4H Lite ls 7 isomers, mean 1.2 
CyHy oC yy *+3H 14 2 isomers 

C\Hy oC Hy +3 1.42 to 1.46 4 isomers, meawa 1.44 
CyHeo Cs +3H 1.2 and 1.35 2 isomers 

CyoHy oCgHy 43H LA Benzene, appearance 


potential 39.8 volts 


CyHy oCyH 2H 13 2 isomers 

CysHeo Csi +2H 10and 1.2 2 isomers 

CyHy oCsHye* +2H Liltols } isomers, mean 1.2 

CH, CoH H l Ethylene, appearance, 
potential 36 volts 

Coy oC alls H 1.2 Ethane appearance, poten 
tial 32 volts 

CyHe oCsH +H 1.3 and 1.2 2 isomers 

C Heo CH i L2tol4 4 isomers, mean 1.3 

CsHy oCsH +H 1.1 and 1.25 2 isomers 

CH o CoH H 1.2 Benzene, appearance po 


tential 30.4 volts 


Hy oCsHy* 11 Isoprene 
CyHy oC Het 10 Benzene, appearance po- 
tential 27 volts 


to remove each H atom from a hydrocarbon and 
the energy of recombination 2H-+H, is 4.5 electron 
volts, it makes a large difference whether the pairs 
are removed as H, or 2H. As the different radi- 
cals have different ionization potentials and may 
have excess kinetic energy, the appearance poten- 
tials cannot be computed from these considerations 
alone. However, the results of table 11 would 
indicate that dissociation into H, usually occurs 
in double-ionization processes. That is, it  re- 
quires about the same work to remove 4H, and 
SH,+H, 3H, and 2H,+H, 2H, and H,+H. 
From comparison with the few instances in 
which the appearance potential is known, it appears 
that all double ionization processes involving re- 
moval of H or H, require 30° to 36 volts. Processes 
involving removal of H,+H or 2H, require 
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uighly 40 volts. Processes involving removal of 

‘HL.+H and 3H, must fall in a narrow voltage 
inge above 40 volts, whereas removal of 5H, 
iust require nearly 50 volts. 

Molecules with more than five carbon atoms 
-ive values of 7 (70)/i(50) somewhat less than those 
listed in table 11 for removal of a given number of 
Hi atoms. Thus in seven C;H,, eyelies in the 
process C;H,--C;H,**+8H, the ratio ranges 
from 1.3 to 1.5 as compared with 3.5 given in table 
11. In the same isomers C,H,,-C;H,**+6H 
vives a mean ratio of 1.2, as compared with the 
mean 1.8 for the values in table 11 for loss of 6H. 

In table 7 three of the C,H, isomers give large 
doubly charged peaks of mass 97. This could be 
either CJH** or C7H,,**. The ratio i(70)/i(50) is 
1.15 to 1.17, indicating a low appearance potential, 
and for this reason the ion has been ascribed to 
C-H,,**. Removing 15 H atoms would require 
a very high appearance potential. Similarly 
2.2.4,4-tetramethylpentane in table 8 gives ions of 
masses 96, 97,and 98. The ratios are in the range 
1.13 to 1.18, and the ions have been ascribed to 


C,H" re Cri," * and CH" watt 
V. Summary and Discussion 


The doubly charged ion spectra of hydro- 
carbons with three, four, or five carbon atoms 
follow some simple rules that can be summarized 
as follows: 

1. In the most probable ionization processes 
carbon bonds are not broken, but all except two, 
three, or four hydrogen atoms are removed. 

2. The relative intensity distribution is similar 
in all hydrocarbons with the same number of 
carbon atoms. 

3. The magnitude of the peaks is least for 
saturated molecules, greater for monoolefins and 
alkvleyelies and greatest for diolefins and cyclo- 
olefins. 

4. lon peaks tend to alternate in intensity with 
larger values for ions with an even number of H 
atoms, 

These rules are in contrast to the case of 
singly charged ions. For these the most probable 
ionization processes frequently involve breaking 
carbon bonds, the spectra depend considerably 
on the structure of the molecule and are quite 
different for paraffins, olefins, and _ diolefins. 
lons with an odd number of H atoms tend to be 
more abundant. 
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Doubly charged ion spectra of molecules with 
six or more carbon atoms are quite different from 
the lighter hydrocarbons. The most probable 
ionization processes involve breaking of carbon 
bonds in all these except the alkylbenzenes and 
the C; eyelies. Structural differences and the 
degree of saturation give large differences in the 
spectra. C ;H,** is often the maximum peak in C, 
hydrocarbons and in C; paraffins. In other C, 
hydrocarbons and in Cy and Cy hydrocarbons, 
C;H,** is often the maximum peak. Rules 3 and 
4 remain true, in general, with minor exceptions. 

A qualitative consideration of the ionization 
process gives a partial explanation of some of the 
properties of doubly charged ion spectra. The 
least potential required to doubly ionize a large 
molecule would involve removal of two electrons 
from two atoms near each end of the molecule. 
(The Coulomb law of attraction makes the work 
required to remove the second electron least if it 
is removed from an atom as far away as possible 
from the first ion). The resulting ion with a charge 
near each end is subject to large disruptive forces, 
and probably a frequent result of double ionization 
is that the molecule ion dissociates into two singly 
charged ions. In general this process will be 
indistinguishable from single ionization processes, 
but Stevenson and Hipple [8] point out that the 
appearance potential of CH,* in the two butane 
isomers indicates that it comes from the process 


CHa »C,H,* T CH;,* T 2e. 


If, however, a molecule loses most of its hydrogen 
atoms in the dissociation process, then the bonds 
between carbon atoms become double bonds and 
the bonding force is greatly increased. Evidently 
it is then sufficient to overcome the disruptive 
force of the double charge. 

The intensity rules for doubly charged ions in 
molecules with three to five carbon atoms can be 
explained on the hypotheses that all but two to 
four H atoms must be removed to give stable ions, 
and that the probability of this ionization process 
decreases as the appearance potential increases. 
With unsaturated molecules and cyclics the num- 
ber of H atoms to be removed is less, and the 
intensity increases. As hydrogen atoms tend to 
be removed in pairs to give Hy), ions with an even 
number of H atoms are commonly more abundant 
than ions with an odd number. 


377 


oe - are. 


ee 


are 04 © eee ee 














Terres SEs TSS Te. «se 


From table 11 it is seen that removing 10 H 
atoms from C;H, to give C;H,** requires nearly 
50 volts. One reason why different ionization 
processes are involved in molecules with six or 
more carbon atoms is because the work required 
to remove all but two H atoms exceeds the ioniz- 
ing voltage. In the C,H, isomers, C;H,** is 
again the largest peak. This can be obtained by 
the least energy by removing CH, and four H, 
molecules. The fact that C,;H,**, C;H,**, and 
C;H,** are larger than the peaks with an odd 
number of H atoms, supports the view that loss 
of CH, and H, molecules actually occurs. In 
singly charged ion spectra, peaks with an odd 
number of H atoms are larger, and it is known that 
the more probable ionization processes involve 
dissociation with loss of CH, and other radicals 
with an odd number of H atoms. (lonization 
with dissociation into CH, is known to occur but 
is less probable [8]). 

In Cy hydrocarbons doubly charged ions with 
six carbon atoms are found in C,H,» isomers, and 
these ions are the most abundant ions in cyclo- 
hexene, CyHy», and in benzene, C,H,. Cy, ions 
become more abundant as the number of H 
atoms decreases. This is also true of C; ions in 
the C; hydrocarbons. In Cy, hydrocarbons there 
are no doubly charged C, ions in C.H,, and CH, 
isomers. They are strong in 4-ethenyl-1l-cyclo- 
hexene, CH, and predominant in the CyHy,. alkyl- 
benzenes. 

In hydrocarbons with seven or more carbon 
atoms the ion C;H.’* is commonly the most 
abundant ion. The empirical fact that such an 
ion is stable is another factor that makes the 
doubly charged ions of the heavier hydrocarbons 
differ from the lighter ones. Also in the heavier 
hydrocarbons certain configurations of carbon 
atoms seem to favor doubly charged ions such 
as the 2,2,4-trimethyl arrangement in Cy, hydro- 
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carbons. Other configurations, such as ethyl sice 
chains in the octanes, give no doubly charged 
ions. This introduces striking individual dif- 
ferences among isomers that are not found in the 
lighter hydrocarbons. Individual differences, de- 
pending on structure, are the rule in singly 
charged ion spectra. 


VI. Conclusion 


The study of doubly charged ion spectra is of 
interest in that it gives additional evidence as to 
the ionization process. Doubly charged — ion 
peaks are, however, not very useful in the prac- 
tical application of most spectra to chemical 
analysis, although it is essential to know the 
complete spectrum of every pure compound 
analysed. In the case of propadiene and propyne 
there is a striking difference in the doubly charged 
ions that affords the best basis for distinguishing 
the compounds. In nearly every other case the 
peaks are too small to be used in analysis. 
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Heats, Equilibrium Constants, and Free Energies of 
Formation of Cyclopentene and Cyclohexene 


By Morton B. Epstein,* Kenneth S. Pitzer,’ and Frederick D. Rossini ‘ 


For cyclopentene and cyclohexene, values are presented for the following thermody- 


namic properties to 1,500° K: heat-content function, free-energy function, entropy, heat 


content, heat capacity, heat of formation from the elements, free energy of formation from 


the elements, and logarithm of the equilibrium constant of formation from the elements. 


Equilibrium constants are given in graphical form for some reactions of isomerization, 


dehydrogenation, and disproportionation. 


As part of the work of the American Petroleum 
Institute Research Project 44 at the National 
Bureau of Standards and the University of 
California, values have been compiled for the 
thermodynamic properties in the gaseous state 
to 1,500° K of the heat-content function, free- 
energy function, entropy, heat content, heat 
capacity, heat of formation, free energy of forma- 
tion, and logarithm of the equilibrium constant 
of formation for cyclopentene and cyclohexene. 
Calculations have also been made of the free 
energies and equilibrium constants of a number 
of reactions involving isomerization, dehydro- 
genation, and disproportionation of these com- 
pounds, 

The values of the constants used in the present 
calculations are as follows [1,2]°: The calorie is 
the conventional thermochemical calorie defined 
as 4.1840 absolute joules; the absolute temper- 
ature of the ice point is 273.16 +0.010° K; the 
gas constant is 1.98719 +0.00013 cal/deg mole. 

For cyclopentene and cyclohexene, the values 

This investigation was performed as a part of the work of the American 
Petroleum Institute Research Project 44 at the National Bureau of Standards 
ind the University of California 

Research Associate on the American Petroteum Institute Research Project 
i# at the National Bureau of Standards 

Associate Director of the American Petroleum Institute Research Project 
Hat the University of California, Berkeley, Cal 
* Director of the American Petroleum Institute Research Project 44 at the 
National Bureau of Standards 
! Figures in brackets indicate the literature references at the end of this 
paper 
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of the heat-content function, free-energy function, 
and heat capacity were taken from reference {3}. 
As usual, the value of the entropy is given by the 
value of the heat-content function less the value 
of the free-energy function, and the heat content 
is obtained by multiplying the value of the heat- 
content function by the temperature. 

The resulting values of the thermodynamic 
functions are given in table 1, which includes 
values of the following properties from 0° to 
1,500° K: function, 
function, entropy, heat content, and heat capacity. 

Values of the standard heat of formation at 
25° C of gaseous cyclopentene and cyclohexene, 
(solid, 


heat-content free-energy 


from carbon graphite) and hydrogen 
(gaseous), were calculated from the values for the 
standard heat of formation for the liquid state cal- 
culated by Prosen, Yenchius, and Rossini {4] from 
experimental values for the heats of combustion, 
together with values for the standard heat of 
vaporization calculated by Wagman and Rossini 
[6] from vapor pressure data of Forziati and Ros- 
sini |7]. The values for the standard heat of forma- 
tion from the elements, at 25° C, A//f°, from refer- 
ence |4] are as follows: Cyclopentene (liquid), 
1.16 £0.16 keal/mole; eyvclohexene (liquid), —9.70 
0.19 keal/mole. The values for the standard 
heat of vaporization at 25° C, A//e®, from reference 
[6], are as follows: cyclopentene, 6.71 +0.07 
keal/mole; evyclohexene, 8.00 +£0.08 keal mole. 
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The method of calculating values of the standard 
at of formation, the standard free energy of 
ormation, and the logarithm of the equilibrium 
onstant of formation for the different tempera- 
tures in the range 0° to 1,500° K is the same as that 

escribed in section IV, 1 of reference {5}. 

The resulting values for the formation of the 
viven hydrocarbon in the gaseous state, from the 
elements carbon (solid, graphite) and hydrogen 
gaseous), each in its thermodynamic standard 
reference state, are included in table 1, which gives 
values of the following properties to 1,500° K: heat 
of formation, free energy of formation, and loga- 
rithm of the equilibrium constant of formation. 
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In figure 1 are plotted, as a function of tem- 
perature, values of the logarithm of the equilibrium 
constant for the reactions of dehydrogenation of 
cvclopentane to cyclopentene and of cyclohexane 
to cvclohexene, 


TEMPERATURE IN °C 
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Fictre 2 Dehydrogenation of cycloherene to benzene and 


disproportionation of cyclohexe ne to benzene and cyclo- 
herane. 


rhe curves give the logarithm of the equilibrium constant as a function 


of temperature for the following reactions in the gaseous stats 
CH gas, cvclohexene C.Hs (gas, benzene) +2H, (ga 


iCoH gas, cvclohevene CoH, (gas, benzene) +2C Hy (gas, cyclohexane 


In figure 2 are plot ted, as a function of tempera- 
ture, values of the logarithm of the equilibrium 
constant for the dehydrogenation of evelohexene 
to benzene and for the disproportionation of eyclo- 
hexene to evclohexane and benzene. 

In figure 3 are plotted, as a function of tempera- 
ture, values of the logarithm of the equilibrium 
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Preparation of Calcium Melibionate 
By William W. Walton and Horace S. Isbell 


It has been found possible to prepare the comparatively rare salt, calcium melibionate, 


in any desired quantity by electrolytic oxidation of a crude solution of melibiose. The 


melibiose solution, prepared from commercial raffinose, is subjected to electrolytic oxidation 


in the presence of caleium carbonate and calcium bromide, and the calcium melibionate is 


crystallized directly from the electrolyte. 


The optical rotation of a 2-pereent aqueous 


° sa 6 20 
solution of ealeium melibionate, expressed on the anhydrous basis, was found to be fa]p + 


106.2°, 


melibionate forms a basie caleium salt, Cal C),H»,Oj2)2.4Ca0. 


which disagrees with the value, [a]}7-+88.6°, reported in the literature. 


Caleium 


This is analogous to the 


basic calcium salts of lactobioniec and maltobionic acids. 


I. Introduction and Discussion 


The need for calcium melibionate in connection 
with a study of the glycosyl derivatives of ascorbic 
acid led to an investigation of methods available 
Neuberg, Scott, and Lach- 
mann [1] * prepared calcium melibionate by oxidiz- 


for its preparation. 


ing melibiose with bromine, removing the excess 
bromine with a stream of air, neutralizing with 
calcium carbonate, and separating the calcium 
melibionate from the accompanying calcium bro- 
mide by repeated crystallization from ethanol. 
Levene and Wintersteiner [2] treated a solution 
of raffinose with baker’s yeast and malt sprouts to 
This solu- 
tion was purified by treatment with basic lead 


obtain a crude solution of melibiose. 


acetate, hydrogen sulfide, and a decolorizing car- 
bon, and the melibiose in solution was oxidized 
with iodine in the presence of barium iodide and 
A lengthy purification of the 
inelibionic acid involved the use of lead carbonate, 
sulfurie acid, silver sulfate, hydrogen sulfide, and 
barium hydroxide. The calcium. salt 
tained by neutralization of the melibionie acid 
with caletum carbonate and crystallization from 
water by the addition of methanol. 


barium hydroxide. 


was ob- 


Phe subject matter reported in this paper was abstracted from a thesis 
ibmitted by William W. Walton, of this Bureau, in June 1947 to the Faculty 
the Graduate School of the University of Maryland, represented by Dr. 
vathan L. Drake, in partial fulfillment of the requirements for the degree of 
Hoctor of Philosophy. (See also J. Research NBS 41, 119 (1948) RP1910.) 
Figures in brackets indicate the literature references at the end of this 
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obtained crystalline 


then oxidized the 


Levene and Jorpes [3] 


melibiose from raffinose and 
melibiose with bromine in the presence of barium 
benzoate {4}. They purified the melibionie acid 
by successive treatment with silver carbonate, 
hydrogen sulfide, sulfurie acid, and barium hydrox- 
ide, and finally by extraction of benzoie acid 
with chloroform. The solution of melibionie acid 
was boiled with caleium carbonate; the resulting 
calcium melibionate was precipitated by adding 
methanol and was reerystallized by dissolving in 
water and adding methanol. 

The laborious purification steps required in all 
of these methods of preparation made desirable a 
simpler method for the preparation of caleium 
melibionate. Isbell and Frush [5] have developed 
a convenient method for the preparation of aldonic 
acids from aldoses by electrolytic oxidation. In 
the case of glucose, crude solutions can be oxidized 
and the acid separated directly from the elec- 
trolyte as the caleium salt. It seemed possible 
that caleium melibionate might be prepared 
by this method directly from the liquor that re- 
mained after hydrolysis of raflinose and fermenta- 
tion of the resulting fructose. 

Application of the electrolytic method to the 
preparation of calcium melibionate from melibiose 
and calcium carbonate gave a product that erys- 
tallized directly from the elect rolyte as the reaction 
The crystals were easily separated by 
By returning the mother liquor to the 


proceeded . 


filtration. 
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electrolytic cell with additional melibiose and cal- 
cium carbonate, the salt can be prepared in a con- 
tinuous process similar to that used for the produe- 
tion of caleium gluconate [6]. The relatively 
small quantities of the salt required for the present 
study, however, were obtained in several separate 
preparations. The crude product was readily re- 
crystallized from hot water 

Crystalline calcium melibionate is reported in 
the literature to contain from 11 to 15 percent of 
water. The product prepared by electrolytic 
oxidation varied in water content from 6.6 to 13.9 
percent, and it was found that air-dried material 
lost or gained weight from day to day, apparently 
from changes in the relative humidity. An 
anhydrous sample gained more than 10 percent 
in Weight when stored in a vacuum oven at 50°C 
but was restored to its anhydrous condition when 
dried in vacuum over anhydrous magnesium per- 
chlorate at room temperature. The specific rota- 
tion, [a)7’, of the pure salt in 2-percent aqueous 
solution was found to be + 106.2°, expressed on an 
The specific rotation is recorded 
in the literature as +-88.6° at a temperature of 
17° CC. Unfortunately, the wavelength of the 
light used is not given in the original publication 
[1], but the constant is quoted in various hand- 
books [7, 8] as [a]}. 

Inasmuch as basic calcium salts are useful for 
the separation and purification of a number of 
sugar acids [9], it seemed desirable to ascertain 
whether melibionie acid also forms such a salt. 
The basic caleium salts of the monobasic acids 


anhydrous basis 


derived from monosaccharides are of the type 
Ca(anion),.2CaQO. 
charides, lactobionie acid, and maltobionic acid, 


Two acids derived from disac- 


however, have vielded basic calcium salts having the 
formula Ca(anion),.4CaQO. A basic calcium salt of 
melibionic acid has now been prepared and found 
composition Ca(C),H,,O,),.4Ca0, 
which corresponds to the basic salts prepared from 
The new sub- 


to have the 


lactobionic and maltobionic acids. 
stance is difficultly soluble and useful for the 
preparation and purification of melibionie acid. 


II. Experimental Details 
1. Preparation of a Solution of Melibiose 


Six hundred grams of raffinose (purity 99.08 ©&% 
based on optical rotation), equivalent to 1 mole of 
pure raffinose, was dissolved in 3.6 liters of water 
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containing 60 ml of nutrient solution’ Fo. 
drops of a I:1 aqueous acetic acid solution w.s 
added to bring the pH of the solution to 4.5 
The optical rotation of the solution was found |o 
be 90.3° S (2-dm tube). Twenty milliliters of 
invertase solution * was added, and the mixture 
was allowed to stand at room temperature for 24 
hr. The observed rotation, based on the volume 
before the addition of the invertase solution, was 
47.8° S (2-dm tube). This value indicates that 
the sugar had been converted to melibiose and 
fructose to the extent of 97.4 percent. 

To remove the fructose, two cakes of baker's 
veast were added, and the solution was allowed to 
ferment for 3 days at room temperature. The 
optical rotation at the end of this period was 
70.5° S (2-dm tube). For complete hydrolysis 
and destruction of the fructose, the calculated 
rotation should be 70.7° S. The rotation ob- 
tained indicates that the conversion of the raffinose 
to melibiose and removal of the fructose was sub- 
stantially quantitative. 

The solution was heated to boiling to destroy 
the veast, and clarified by the addition of 30 ml of a 
decolorizing carbon, and filtration. The filtrate 
was concentrated under reduced pressure to 1.S 
liters and used directly for the preparation of 
calcium melibionate. 


2. Preparation of Calcium Melibionate 


The solution of melibiose obtained from 600 ¢ 
of raffinose was transferred to a 3-liter beaker, and 
32 ¢ of calcium bromide and 55 g of calcium ear- 
bonate were added. Two graphite electrodes, 22 
mm in diameter, were immersed to a depth of 12 
em, and the solution was stirred mechanically 
while a current of 1 ampere at 10.8 v was passed 
through. The electrolysis was stopped when 
titration of an aliquot with iodine [10] indicated 
that 97.3 percent of the melibiose had been ox: 
dized, 
Higher current 
lower percentages of oxidation. 

The solution was heated to dissolve the portion 


The current efficiency was 82.9 percent 


efficiencies were obtained for 


of the calcium melibionate that had crystallized 
The filtrate was 
concentrated under reduced pressure to a volumy 


and the hot solution was filtered. 


The nutrient solution was made by dissolving 2.5 g of NHgNOy:, 0 
of KH oP O4, and 2.52 of MgS04.7 HyO in 1 liter of water 
‘ The invertase solution is a commercial product made by Difco Laborat 
ies, Detroit, Mich 
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yout SOO ml and allowed to cool to room tem- 
ture. The calcium melibionate that erystal- 
| from the solution was collected on a Biichner 
el and washed with three 25-ml portions of 
twater. Three hundred and twenty-five grams 
rystalline calcium melibionate, calculated on 
The filtrate 
concentrated under reduced pressure to 100 
An addi- 
al crop of 40 ¢ of calcium melibionate, caleu- 
The 


vield of crude product Was 365 gy or 6.7 


anhydrous basis, was obtained. 
and 75 ml of methanol was added. 
ed on the anhydrous basis, was obtained. 


cent. Larger quantities of raffinose were used 
) number of preparations without any difficulty, 
hon analysis, the crude material was recrystal- 
d several times from hot water, and finally was 
dried. The water content of the material was 
termined by drving samples to constant weight 
room temperature over anhydrous magnesium 
perchlorate, and the results reported below are 
ypressed on the anhydrous basis. 
Caleulated for anhydrous Ca(C).H.,- 


(nalvsis 
‘a, 4 38.2; H, 5.6. 


() > B3hok: 


CC: 33.3: 37 


Found: Ca, 5.33: 


lelp 106.2° (water, ¢=2). 


3. Preparation of Basic Calcium Melibionate, 


Ca(C,.H,,O,,),.4CaO 


Ten grams of calcium melibionate was dissolved 
n 200 ml of water, and the solution was cooled to 
Qe? ( 
200 ml of water at O° ¢ 


A suspension of 3.5 g of calcium oxide in 
‘was added, and the mix- 
ture Was shaken for several minutes. It was then 
filtered into a flask through which a current of 
The clear 


solution was concentrated under reduced pressure 


carbon-dioxide-free air) was passed. 
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to 150 ml and finally heated for a few minutes on 
the steam bath to complete the precipitation of 
the basie salt. The salt was separated by filtra- 
tion in the absence of carbon dioxide and dried at 
60° Cin vaeuum. Analysis: Caleulated — for 
Ca(C).H2,O;2)2.4CaO: total Ca, 20.47; CaO, 22.91 


Found: total Ca, 20.45: CaO, 22.56 


The authors express their appreciation to Naney 
B. Holt for the measurements of optical rotations 
and to Rolf A. Paulson for the determinations of 
carbon and hydrogen reported in this paper 
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Polarographic Limiting Currents 


By John Keenan Taylor, Roberta E. Smith, and Irvin L. Cooter 


A cathode-ray oscillograph and specia 


I. Introduction 


theory of currents limited by the rate of 


son of an electrolyzable constituent of a 
on to the electrode has been developed by 
s lon linear diffusion to a plane electrode 
theoretical relation between the instantaneous 


; 


nt and tim is given hy / 


iz 
A nkCt 
V; 


lis the area of the electrode, 77 is the diffu- 
constant, Cis the concentration of the diffus- 
material in the bulk of the solution, » is the 
iber of electrons taking part in the electrode 
on, and F is the 


esponding equation for 


Faraday constant The 
v spherical electrode 


ed radius, Ry {1, 3) is 


tel One teh FI 
: nkC't rRnk DC. 2 
V; 


similar equation for the dropping-mereury 


trode tS] is 


[7] 2] - nPDeCrZt”’, ; 


e dis the density of the mereury and m is 
mass of mere ury emitted by the electrode in 


t time 


larographic Limiting Currents 


Integration over the lf span ob a drop of the 


values of viven 1yy eq » TIVES the expression 


hy I) Com 


where is the average current, fis the drop time 


and / expresses all of the numerical constants 


Equation 4 is used widel 


Vv in polarographic meas 
urements and referred to as the Ilkovie equation 

Relations between diffusion current and time 
have been investigated experimentally. Laitinen 


and Kolthoff [8] found that eq | expresses the 
eurrent for plane ele trodes provided care Is tuken 
lo eliminate the disturbing effects of convection 


and adverse density eradients However these 


workers were not able to obtain the predicted 


behavior with spherical electrodes Ih aqueous 
solutions and only approached it in gelatin gels 
Where disturbances are minimized by the high 
viscosity of the medium 

No conelusive studies have been made of the 
current-time relationship for the dropping-merecury 
Iikovie [6] made some measurements 


and states that the 


electrode 
data are represented by a 
parabola ol approximately the sixth order Hlow- 
ever, he does not give sufficient details of his ex- 
periments to permit a recalculation of the data 
Most of the tests of the Hkovi equation fall into 
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two groups. The first consists in measurements or 
calculation of all of the parameters of the equation 
and comparing the diffusion current caleulated 
with that observed [7]. The second consists in 


determining the constaney of the diffusion current 


constant, J=74 Cm? or’®, as ©, m, and ¢ are varied 
S|. All of these experiments have been performed 
While using highly damped apparatus for the 


electrical measurements, the assumption being 
made that average currents were measured 

The present paper reports the results of some 
experiments designed to determine the current- 
time relationship for the dropping-mercury elec- 
trode Because of their practical interest, the 
relations between the average currents, and those 
measured by several polarographie circuits also 


have been determined 


II. Experimental Details 


A schematic diagram of the cireuits and appara- 
tus used for measuring the polarographic diffusion 
currents as a function of time is shown in figure 1. 
In this figure, P? 
potentiometer for supplying the voltage for the 


isa Leeds & Northrup portable 


circuit, Which consists of the polarographic cell 
When the 


is on the “measure” position, 


and the calibrated resistance box, R). 
selector switeh, S), 
after amplification, 
A pho- 


tographie record of the deflection of the beam as a 


the voltage drop across A), 


deflects the beam of a cathode-ray tube 


funetion of time is made by means of the drum 
eumera 

As the voltage across R, is only a few millivolts, 
and the voltage required on the deflection plates of 
the oscillograph is approximately 100° volts, an 
amplification factor of several thousand is required. 
Moreover, this amplification factor should remain 
constant over the range of applied voltage and 
during the recording time of several seconds. 
These requirements were met by using a direct- 


David W. 
This amplifier, 


coupled amplifier designed at the 
Tavlor Model Basin [9] 


has only resistance coupling so that it can be used 


which 


with aperiodic phenomena, has an amplification 
factor of 1,000, 2,000, or 5,000, with uniform re- 
sponse in the frequency range from 0 to 50,000 ¢ s. 
The output of this amplifier is connected directly 
to the vertical deflection plates of a Du Mont 
cathode-ray oscillograph, tvpe 208. The photo- 
graphic record of the deflection of the beam of the 
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Le M 
kigtre 1 Ipparatus used for measuring lin ngc 
P,, P, | I I Ri, R2, ‘ t / 
tor: Ra abiliz r: Ss . Ss, ta k 
battery 1, a f 0 aph. ¢ ‘ a: ¢ 
\f, mot 


cathode-ray tube is made by a drum = camera 
Essentially, this camera consists of a lens an 
shutter mounted in a light-tight housing enclosing 
a drum having a circumference of 60 in. and 

rim 4in. wide. The drum is mounted on a vertica 
axle and is belt-driven by a direct-current moto! 


When the 


rotating at 2 rpm, the peripheral speed is approxi- 


s 


through reduction drum 


gears. 
mately 2 in. sec. The housing can be opened to 
permit the mounting of a film on the drum. 

For the electrodes used in these experiments 
the time of each drop was approximately 3.5 s 
so that the diffusion current record for each droj 
extended about 7 in. on the film; this permitte 
records of seven or eight successive drops on eae! 
film. 

In order that the recorded amplitude of | 
diffusion current at any point represent the t: 
amplitude, it is necessary that the reference | 
for zero current remain constant during the 
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for each drop Since there is a small but 
ficant drift of the zero position of the amplifier 

the 30 see of recording, it is necessary to 
da sufficient number of zero positions dur- 
his interval to establish a reliable reference 

A simple means of doing this is to incor- 
te a tapping key, S,, into the input eireuit of 
between R and R 


tion for this switeh is closed, but by tapping 


amplifier The normal 
vey the external cireuit of the amplifier can 
periodically opened for a short time interval. 
| ul interpolation between the various zero 
ss was used to determine the zero level cor- 
wnding to a given deflection 
The correlation of the recorded amplitude with 
voltage producing the deflection Was obtained 
means of the calibrating circuit, consisting of 
battery, Boa current-limiting resistor, Ry, and a 


box, R. With 


“calibrate” position, 


brated resistance selector 


teh, Sj. on the known 
ltave Was applied to the amplifier by means of 
calibrating circuit. By means of the current- 
ting resistor, the voltage drop across Ry, was 
isted so that measurable potentials are applied 
the amplifier input. These potentials were so 
eted that, without changing the gain controls 
the amplifier, the deflection due to the eali- 
ting voltage covered, mn several steps, the 
itimum range of the deflection due to the polaro- 
aphie current. These calibrating voltages were 

usured with the precision potentiometer, P. 
The deflection of 


ith due lo 


the cathode-ray oscillograph 


these calibration voltages was re- 
vded on the same film immediately after the re- 
ling of the polarographie current. From = the 
enitude of the calibrating voltages and the re- 
ne measured deflections, a calibration curve 
V be obtained by plotting sensitivity. (my em 
I respect to the deflection in centimeters of 
oscilloscope From the calibration curve, the 
sured deflection due to the polarographic ecur- 
and the value of R, and Ry, it is possible to 
late the diffusion current. 
\ calibration factor for converting length meas- 
ment along the film to time was obtained by 
usurmng the average drop raie of the electrode 
| dividing this by the average measured length 
the current-time curve for a number of drops 
s quotient, when multiplied by the measured 
time interval 


eth interval, would give the 


re was taken to make the vertical axis of the 
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oscillograph-beam deflections parallel to the aNnis 
of rotation of the drum in order to eliminate timing 
errors, Which might otherwise occur. 

The current-time records were measured on a 
special comparator The comparator consisted of 
a heavy metal bed, 4 ft. long on which rested a 
The 
screw with a graduated head, had a range of about 
OO cm of 


traveling carriage carriage, driven by a 


longitudinal motion. Two long glass 
plates for holding the film, thiminated from below 
by a green fluorescent light, were set about 2! in 
below the carriage and parellel to it. On the 
carriage itself was the observing microscope, sup- 
The 


contained the 


ported on a micrometer. slide observing 


microscope, Of magnification 4 
reference cross hairs. By the use of this compara- 


tor, values of the abseissa (time) and ordinate 


current) could be measured with a precision of a 
few microns 

The solution used for the oscillographic measure- 
ments was 3.018 millimolar in cadmium, 0.1 molat 
in potassium chloride, and contained O.OL percent 
of gelatin to suppress maxima. A similar solution, 
except for the cadmium concentration, which was 
2 724 millimolar, was used in all of the other work 
These solutions were made with known amounts 
of cadmium metal and checked by analysis 

\ silver wire wrapped around the dropping- 
mereury electrode served as the anode in all but 
a few cases where the Tereury pool Was emploved 
cadmium 


Diffusion-current measurements — for 


were made at the potential —-0.9 v with respect to 
the silver silver-chloride anode The droptime t 
of the electrode was determined at this potential 
as Was also the flow, m, of mercury from the eapil- 
larv. It was found that the value of m in’ the 
solution and at this value of the applied potential 
was about 4 percent smaller than when measured 
in air, in agreement with previous investigators 
10] 

Purified nitrogen was used to free all of the solu- 


tions of oxygen. The diffusion-current measure- 


ments were made ina thermostated bath at 25° C 


III. Results 


A number of current-time curves were obtained 
in Which the zero was marked on the film after the 
recording of current, but the slight drift of the 
amplifier was sufficient to preclude their use in 
quantitative caleulations. The qualitative con- 
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clusions drawn from them are, however, in agree- 
ment with and support the more precise observa- 
tions. 

Several experiments were performed in’ which 
the zero was put on the film during the measure- 
ments of the limiting current. However, wherever 
a zero mark occurred, no record of the limiting 
current was obtained, so that some of the records 
were spoiled. Thus the 
particular drop, or the portion of the curve where 


start, or ending of a 


the current was changing most rapidly might be 
Any sufficient 


discarding the 


missing. of these defects was 


reason for measurement. As a 
result, only the complete records of two drops were 
accepted and used for quantitative interpretation. 
A number of the other traces were used for meas- 
urements of the value of the maximum current, 
however, or for the study of selected portions of 
the current-time curve. 

A portion of the photographie record is shown in 
figure 2. The ordinate increases in the usual 
direction, but the time axis increases in the right 
to left Three traces of the residual 


current are shown in the upper part of the figure. 


direction 


The current will be observed to start at a high 
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value and to decrease with time, the zero trace 
being always below the residual current line. “Thy 
polarographic limiting current is shown in the 
lower part of the figure. This current is not zero 
at the start of the drop but goes through a mini- 
mum, and then increases rapidly at first and more 
slowly later until the drop is disengaged from thy 


capillary. Figures 3 and 4 are large scale plots o! 


pa 
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Figure 3. Plot of residual current 


current and limiting current, 


The initial behavior of the limiting 


the residual 
spectively. 

current is due to the residual current being great 
than the limiting current during the very ea 
When the instanta: 
is subtract 


period of drop formation. 
ous value of the residual current 


from the corresponding value for the limiting 
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rent, the resulting diffusion current starts at 
o within the error of the experiment. 

The photographie records were measured in the 
nner deseribed, and the abstracted data are 
enoin tables | and 2.) The first column con- 
us the abseissa readings in turns of the lead screw 


the comparator, and the second column gives 
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I ABLS 2 Va iation of diffusion current with time tor 
drop two 
t= 3.488 sec =? 305 mg/sec: ¢ .O18 millimole/liter 


Limiting Residual Diffusion 
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the corresponding deflection of the oseilloscopie 
spot from its zero position. The third column 
gives the elapsed times from the beginning of the 
drop, which were obtained by multiplying column 
one by the abscissa factor. The values for the 
limiting current, /,, in column four resulted from 
multiplying the deflection by the corresponding 
sensitivity for that deflection found from the eali- 
bration curve. The residual current, computed 
ina like manner, was plotted on a large scale as a 
function of time, and values corresponding to the 
tabulated values of time were read and recorded 
in column five. Values for the diffusion current, 
‘a. IN column six are the differences of columns 
four and five. 

Figure 5 is a plot of the logarithm of the current 
with respect to the logarithm of the time. For 
conformance with the Ikovie equation, the plot 
should be linear with a slope numerically equa! to 
one-sixth. The plot is not linear, and a comparison 
with the slopes drawn in the figure shows that it 
never reached this value at any time during the 
life of the drop. This is emphasized by the mathe- 
he", in 
its logarithmic form, was adjusted to the data by 
When the entire 
range of the measurements was used, the value for 


matical treatment in which the function /, 


the method of least) squares. 
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exponent “—O.82 was obtained, but the fit’ was 
extremely poor. In the range 0.1 sec to maturity, 
the value, »—0.51, was found with an average 
deviation of 1.5 percent. For all times greater 
than 0.5 see there resulted »—0.249, with an 
average deviation of +£0.97 percent, whereas in 
the more limited range of 1 sec or greater, 7 
0.227, and the average deviation was +0.70 per- 
cent. If times greater than 2 see are used for the 
calculation, the exponent became n—=0.186, and 
the average deviation was +0. 30 percent. A 
number of attempts were made to fit other simple 
equations of various types to the data, but none 
was found satisfactory 

The average current is of interest, since this is 
the quantity to be expected to be measured when 
using highly damped = galvanometers or other 
measuring instruments. Accordingly, large-scale 
plots were made of the data for two drops, and the 
average current was found by graphical integra- 
tion. For drops having maximum currents of 
27.37 and 27.70 wa, respectively, the corresponding 
average currents were 22.1; and 22.3,ua The 
ratios of the average to the maximum current 
were O.S1, and 0.80, in the two cases, which are to 
be compared with the value 0.859 (6/7) found from 
integration of the Ilkovie equation. Integration 
of the equation ty kt? gives O.SO,; for the ratio, 
in good agreement with the graphical integration. 
Although this latter equation does not fit the data 
during the early time of drop formation, the bulk 
of the area under the curve oceurs during the time 
when the fit is good, which explains the agreement 
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of the value found by integration with that ob. 
tained by graphical means. 

Calculation of the diffusion-current cons! 
[11], I. defined by the expression I tg/Cm? “t 
is a satisfactory means of reducing experiment. 
done under slightly different conditions 
common basis, since this expression has | 
shown to hold true for widely varying charact, 
istics of the cleetrode When the maxin 
current is inserted in the expression, the consta 
Diggs 4.25 ts obtained whereas the average cur 
vields the value /-— 3.44. 

Less precise but confirmatory evidence for ¢| 
deviation of the experimental current-time cury 
from that of the Ikovie equation was obtain 
from measurements with a Hathaway galy; 
nometer oscillograph, model 58 B.A’ cha 
speed of 28) mm see was used, and the galy 
nometer had a sensitivity of 0.38 ga imm and 
frequency response of 20 es, permitting a fait! 


Tapie 3 Ditfusion-current constants from measure 


with the gqalvanometer osc illograph 


Solution Tas l ie | 
(ad i io s 
Lead is! “ “i 


ful record of the current for all but about t 
first 0.2-sec period of the life of the drop. TI 
apparatus Was calibrated in a manner similar 

that used in the work with the cathode-ra 
oscillograph, and diffusion-current) measurement 
were made on the cadmium solution alread 
deseribed and on a solution 3 millimolar in lea: 
with I-N0 hydrochloric acid as the supporting 
electrolyte. From values of the maximum cur 
rent and of average current obtained from graphic: 
integration of the curves, extrapolated to zer 
time, diffusion-current constants were calculate: 
The results are given in table 3.) (The extra 
rte 


polation introduces no significant error in 
gration, since the area under this part of the cu 
isa small fraction of the total area.) 

The residual current, 7,, may be considered 
be composed of a faradayvie current, 7,, and a « 
denser current, ¢,. The former is the sum of th 
diffusion currents due to reducible impurities 
the solution, whereas the latter arises as the resu 
of charging or discharging the condenser resulting 
from the double laver at the mereury-solutio! 
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rface. The condenser current at any time, 


ving the life of a drop should be expressed by 
kt *, where & represents the product of the 
meity of the double layer, the difference between 
applied potential and that of the electrocapil- 
, maximum, and geometrical constants [12, 13]. 
cordingly a plot of the logarithm of the residual 
rent with respect to the logarithm of the time 
ould be linear with slope equal to minus one- 
rd. Such a plot is given in figure 6 and is seen 
» be linear, but the slope is only one-half of the 
eoretical value. Although pure materials were 
d for the supporting electrolyte and care was 
ercised to remove oxygen from the solutions, 
: possible that there may have been some con- 
bution of ¢, to the residual current. The 
wesent experiments do not permit conclusions as 
to the reason for the discrepancy. 
lt is of interest to determine the magnitude of 
the current measured by various types of polaro- 
vraphic apparatus and to study how this depends 
ipon the constants of the measuring instrument. 
The simplest apparatus is one in which a refleet- 
vy galvanometer is placed in series with the cell 
and the deflection read with the conventional 
lamp and seale. It is well known that galvanom- 
elers may be used to integrate the current. In 
the measurement of asymmetrical periodic phe- 
nomena with highly damped instruments or ones 
whose period during operation is long with respect 
to that of the event observed, the excursions of 
the galvanometer spot around an average position 
are small when compared with the total deflection, 
and the midpoint of the oscillation should approxi- 
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TABLE 4 Dependence of the diffusion current on the 


characteristics of the galvanomete ' 


In column 2, the values in italic type indieate the critical damping resistance 
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mate the average current. With relatively short- 
period equipment, this method is no longer ade- 
quate. 

The results of measurements with several gal- 
vanometers are reported in table 4. Galvanom- 
eters with the periods indicated in column 1 were 
damped with the resistances across the terminals 


shown in column 2. The average deflection is 


shown in column 3, whereas column 4 gives the 
values of the corresponding total oscillations of 
the galvanometer spot. The diffusion-current 
constants of column 5 were calculated by using 
the midpoints of the oscillations as the average 
current. For the short-period galvanometers, 
an attempt was made to determine the variation 
of the galvanometer oscillation with time by 
observing the time when the galvanometer spot 
crossed selected points on the seale. The result- 
ing data were plotted, integrated graphically, 
and the values in column 6 were calculated by 
using the average currents so obtained. 

It is interesting to note that those instruments 
that had a period several times that of the drop- 
ping electrode or were highly damped gave essen- 
tially the same values for J when calculated from 
the midpoints of the galvanometer oscillations. 
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The galvanometer of period 1.85 see was of the 
enclosed spotlight type with a short seale, so that 
only small deflections could be observed with a 
corresponding low precision. The data, however, 
are conclusive in showing that, except for the 
highly damped condition, low values are obtained 
when the midpoints of the oscillations are used to 
calculate the current. If the wave form of the 
oscillation is known, however, the integrated 
current is in satisfactory agreement with that 
found with galvanometers of relatively long period, 
and with the true average current found in the 
studies with the oscillograph. 

Diffusion currents can be determined from meas- 
urement of the potential drop across a known re- 
sistance. This potential varies throughout the life 
of the drop, and hence an average value is usually 
measured. The potentiometric balance, using as 
a criterion the symmetrical oscillations of the gal- 
vanometer, is influenced by the characteristics of 
Data illustrating this 
The spotlight 


the indicating instrument. 
dependence ure given in table 5. 


Dependence of the 


potentiometrically on the characteristics of the qalvranameter 


Tariue 5 d {fusion current measured 


Thal 

la I l — mn I 
ohm “ “ 

; ae) 

1] 1 24 

“ s Is 4 1 

‘ - 4 1 ws 

4 Ws ") 

l y a ‘ 

a1 i 


galvanometer of period 1.85 see was used and 
damped as indicated. The diffusion-current con- 
stant found for the highly damped experiments 
was 3.44, which is in excellent agreement with that 
On the basis of 


these and other experiments not reported, it may 


found for the average current. 


be concluded that the same considerations as found 
for the direet measurement of current apply to 
galvanometers used to indicate potentiometric 
balance. 

An interesting confirmation of the potentio- 
metric experiments is given by some measurements 
with the Speedomax potentiometer to record the 
potential drop across a 200-ohm resistance in series 
with the cell. 
to complete a full-scale deflection of 25 em and 


The apparatus used required 1 see 
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In the particu a 
experiments a maximum deflection of 24 em y as 


had a chart speed of 4 in./min. 


obtained, and the magnitude of the oscillations 
was 9 em. By using the average oscillation ‘6 
determining the average current, the low value o/ 
the diffusion-current constant of 3.32 was obtained 
After graphical integration of the potential-tiin 
curves, a diffusion-current constant of 3.49 was 
found, which is in fair agreement with the tru 
average value, considering the precision of th 
measurement, 

Several models of commercial polarographs wey 
used to measure the diffusion current of the cad 
mium solution used in the experiments with th 
oscillograph. The current sensitivities of thes: 
were determined under the exact conditions of us: 
by measuring the recorded deflection produce 
by a known current, the latter being determine 
by measuring the potential drop across a standar 
resistance. Corrections for residual current we 
made by both graphical and exact methods. Wit! 
the Hevrovsky model VIEL instrument the valu 
for J 


Hevrovsky model XII polarograph a large num 
: \ | | 


3.42 was found, whereas with a Sargent- 


of measurements gave [—3.43. 

A Sargent model XX instrument was also use: 
for determining diffusion currents. This appa- 
ratus is provided with circuit elements to dam 
electrically the oscillations of the recorder cause: 
by the periodic growth of the drop. The degre: 
of damping may be varied from position 1, un- 
damped except for the natural period of the re- 
corder, to position 5, which introduces the maxi- 
mum damping. The cadmium solution already 
described and a solution 3 millimolar in lead in a 
supporting electrolyte of 1-N hydrochloric acid 
were studied, and the average of a number of 
determinations gave the results in table 6. The 


PasLe 6 Diffusion-« urrent constants fram measuremer 


with a Sargent model XN polarograph 


Damping position 


Solutior 


Cadmium 2 4s 


Lead 1s as 01 105 105 


constants for both solutions are somewhat higher 
than those found for the other types of apparatus 
The circuit constants were such that even with th 
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iximum damping, large oscillations of the re- 
der were observed Lingane [14] has reported 
iat this apparatus gives high values for diffusion 
irrent under such a condition, which may explain 
he results found. 

Schulman Battey, and Jelatis [15] have de- 
ribed apparatus that incorporates a current re- 
corder of low damping. An instrument of this 
esign, now marketed under the trade name, 
Polaro-analyzer, was loaned to the authors for 
his work. The Esterline-Angus recorder used 
ves a full-scale deflection inh about 0) ° sec 
\lthough this is too slow to give a faithful record 
f the current during the early period of formation 
fa drop, the response is adequate to follow the 
irrent in the later stages of drop growth, when 
ising drop times of the order of 3 see. With the 
millimolar cadmium solution, the value Jima. 
{11 was found. By using the damping cireuit 
vith which the apparatus is equipped, the ratio 
of the average to the maximum current was found 
to be O.S2, in satisfactory agreement with the 
values found when using the oscillograph. It is 
concluded that the peak current found with this 
ipparatus corresponds to the maximum current 


within a few pereent 


IV. Conclusions 


The results of this investigation are conclu- 
sive in showing that the Ikovie equation does not 
represent the instantaneous current for the drop- 
ping-mereury electrode. The deviation is greatest 
at the start of drop formation but is appreciable 
at all times Part of the deviation aay be due to 
an unsymmetrical diffusion field resulting from 
mechanical screening by the body of the capillary. 
Density gradients at the surface of the electrode 
may exert disturbances similar to those observed 
by Laitinen and Kolthoff for solid spherical elec- 
trodes. Experiments are now being planned in 
Which capillaries with substantially the same drop- 
rate but with different values of m will be studied 
These and other investigations may vield valuable 
nformation concerning the nature of the devia- 
Lions. 

Although the Ikovie equation does not represent 
he instantaneous current, several investigators 
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have shown that it explains, quantitatively, the 
variation of the diffusion current with concentra- 
tion of electrolyzable constituent and with the 
characteristics of the electrode, provided a sufli- 
cient concentration of suppresser of maxima is 
present [S] and electrodes with very small drop 
times are not used [16] The conditions vovern- 
ing the quantitative use of diffusion-current con- 
stants have been reviewed in a previous paper 
[17]. To these should be added the precaution to 
make sure that precise values of either maximum 
or average current are measured by the upparatus., 
A preferable procedure is for each worker to deter- 
mine his own values for the constants to be used 
The individual idiosyncrasies of the equipment and 
personal errors of the observer have the least effect 


in this manner. 
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a 
Stress Distribution Near Reinforced Circular Hole 
Loaded by Pin 
By Samuel Levy and Frank C. Smith 
This paper presents a theoretical analysis for the stress distribution in a flat plate near a 
reinforeed circular hole loaded by a pin, and a comparison of the theory with test results 
for a plate of sandwich construction It is assumed that an auxiliary reinforcement at the 
edge of the hole is so rigid that no distortion of the shape of the hole occurs due to the load 
he plate is reinforced in the vicinity of the hole by cireular “doubler” plates Phe theory 
and experiment are found to be in good agreement 
I. Introduction Closely allied with the problem of the stress 
distribution in a flat plate near a reinforced cir- 
his paper presents a theoretical analysis for cular hole loaded by a pin is that of an unloaded 
the stress distribution in a flat plate near a rein- reinforced circular hole in a plate under stress 
reed circular hole loaded by a pin, as shown in This problem has been solved approximately by 
cure L, anda comparison of the theory with test Timoshenko 2), who treated the reinforcement as 
sults for a plate of sandwich construction, [t is a curved beam of constent section. Other solu 
ssumed that an auniliarv reinforcement at the tions of the unloaded hole problem based on the 
ss of the hole is so rigid that no distortion of plane stress theory were obtained by Sezawa and 
the shape of the hale occurs due to the load. The Kubo [3], by Gurney [4], and by Beskin [5] 
plate is reinforced in the vicinity of the hole by Another related problem is that of the stress 
cireular “doubler” plates, which increase — the distribution ina flat plate containing an elliptical 
thickness by a factor 8. The analysis proceeds region filled with material having a greater stiff- 
om the assumption that the reinforced sheet ts ness than the material in the rest of the plate \ 
a condition of plane Stress This requires the solution to this problem was obtained Iy Donnell 
nforcement to be symmetrical on both sides of 6] by using a method that is a generalization of 
he sheet to eliminate bending stresses, and it Inglis’ solution for the elliptical hole 
quires the thickness to be small compared to In the present paper, use is made of the general 
he other dimensions so that the stress can be plane stress solution in polar coordinates The 
onsidered constant across the thickness. The solution given is approximate in that it does not 
tress cannot be assumed constant for points a give zero stress on the free boundaries of the plate 
listance less than one thickness from the boundary However. there is reason to believe that for plates 
tween doubler plate and sheet with moderately small holes, e. g., plates in which 
Previous work on the problem of a plate with a the hole reduces the section of the plate 15 percent 
reular hole loaded by a pin includes a photo- or Jess, this solution gives a good approximation 
astic investigation by Coker and an empirical of the stress distribution, particularly in’ the 
inalysis by Bickley [1]. These investigations neighborhood of the hole 
vere confined to unreinforced holes. It was decided to check the theory by comparing 
ted } th International ( for Applied Mecha computed strains with values measured in tests 
: , on a specimen of Metalite sandwich construction 
. under concentrated load acting in the plane of the 
Hole Loaded by Pin 397 
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plate. 
indicate that reinforcements in this type of con- 


It was hoped that the comparison would 


struction could be designed rationally on the basis 
of the plane stress theory. 


II. Boundary Conditions 


First, a 
solution is obtained for the stresses and displace- 


The analysis proceeds in two steps. 


ments ina cireular ring of constant thickness and 
any inner and outer diameters. Then the solu- 
tion for the reinforced specimen is obtained by 
considering the reinforced parts as one such ring 
and the rest of the sheet as another. Syvmbols 
having the subseript “2"" correspond to the rein- 
forced sheet, whereas S\ mbols having the subseript 
vs! Vip correspond to the unreinforeed sheet. Svim- 
bols without a subscript are general and apply to 
any ring of constant thickness. 

The origin of coordinates, figure 1, is taken at 


the center of the hole with x, yYoas rectangular 
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coordinates for points in the sheet and 7.@as polar 
coordinates. The positive direction of « will be 
taken in the direction of the applied load at the 


The relations between the coordinates are, 
r rcos@ ‘ 
y=r sin 6 


(1) 
r= 72*+y' 





é=tan"' y = 
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Normal and shearing stresses are denoted by th 
Radial displac: 
ments are denoted by u, and circumferential dis 


conventional symbols @ and +r. 


placements by vy. 

Since the edge of the hole is assumed not to 
deform, the boundary conditions for displacement 
at the inner edge of the reinforcement can be take: 
as, 


%.=0, Ff “) 
te—0, r=a ) 
and for forces, 


. 
az 


Bh | 7 


cos 0— Tr, sin 6],,_,,a dé P 


Bh | [o,. SIN 0+ 7,9 COS O],,_ qa da—0. 


At the junction between the reinforced and un- 
reinforced sheet, it is necessary that the displace- 
ments be the same for both the reinforced and 
unreinforced sheet, or, 

t.=-S.; ad) 
aa.) 


ay Ve, Tr 


It is also necessary that, at the junction between 
sheet, the 
normal and shearing forees be the same, or 


the reinforced and the unreinforced 


Sén=€n3 ad) 


Broe—tm; oa) 


The solution of the actual problem, figure 2, a, 
is the superposition of the solutions for the two 
For the prob- 

- 
agmust 


subproblems, figure 2, b, and 2, ¢. 
lem in figure 2, b, the stresses at 7 
satisfy eq 3, and far from the hole, the stresses 
must remain finite, 1. e.; 


o-1, 7, and 7.9, finite; r--@, (6 


For the problem in figure 2, e, the stresses at 
| 4 


r-a must satisfy the condition that the resultant 
load on the hole is zero 


( [o-. COS O— 7, sin O],,_,,d0—0 
0 
| [o,» sin 6+ 7, Cos 4. da—0 
J0 
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URE 2 Solution for the actual problem (a) by the super- 
on of the solutions tor the subp oblems (b) and (le 


nd that, far from the hole, the stress is a uniform 


nsion in the v-direction of magnitude P 2irh, 


here wis the width of the specimen: 


pP nh cos” 
P 2wh) sin? 6 -7—>o, S 

an P 4wh) sin? @ 
For both the problem in figure 2, b, and the 
problem in figure 2, ¢c, the stresses on the free 
ongitudinal edges of the plate should be Zero, 


\uxiliary stress functions of the type used by 
Howland [7] might be used to satisfy this boundary 
condition. This is not done, however, since for 
moderately small values of the ratio a/w of hole 
radius to plate width, the stress functions to be 
used reduce the unbalanced stresses on the free 


edges to a negligible amount. 


III. Solution for Problem in Figure 2, b 


By using the second line of the general stress 


lunetion eq 77 |S, p. 114], 


l—y ! , 
o Myr (og r) cos 6— , # sin @4 
tr 2a 
T ‘ 
cos 4+ cos #4 (9 
477 4r 
vhere 


v= Poisson’s ratio, 

V7, R arbitrary constants to be determined 
from boundary conditions, 

may be shown from eq 34 [8, p. 53] that the 


quations of equilibrium and of compatibility for 
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the stresses in a 


satisfied by 
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(10) 


By using the relations between strains and stresses 


given on pages 62 and 63 of [S], it may be shown 


that the corresponding displacements are 


» » 
-! ” 


»(1—dyp) « 


=V b 


cost logs 


‘os 0 as 
ink 
sin 4 
sin @ log 7 
, 
. iv 
sin @4 ke 


/-= Young's modulus, 
N= integration constant to be determined from 


boundary conditions. 
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Applying the boundary conditions, eq 2, to eq 1] 


M/ 
: ink’ 
\ R) 
tri 
M/ 
irk 
\V/ 
irk 
r 
trl 
N 
tole 
where 
vIVeS 
AM. 
and, 
M,(1-+-2: 


Applying eq 3 to eq 


R 
M, 
J (1 
a 


log a- N, T 
al l Sv 0) 
3+ 2y—Y 
v)i- Ra? 5+ 


10 gives 


AM. P Bh. 


log a Nz / 


v) 


0 


(13) 
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Applying eq 4 to eq 11 gives 


M,(v?—2v—3) log aa + N,4 r l+v)-+ } 


aa 
R,a?a* (1 —3y = 
T. 15 
Wo 2»y—3) log aa N,4 l+p 
aa 
Raa? (1—3v), J 
and 


(16) 
M1 + 2e+-r*) +10, (3 +-2e— Yr { 


log ad N 





]-+» | R cea’ 5 Vv). 
acu J 
Apply Ing eq 5 to eq 10 gives 
» 
p| 3+v)M rT, 4 2R.céa*| 
ou ‘ 7 
» 
) v) M ’ 2R cea’, (17) 
aca~ 
and, 
» 
3| l 1 V/ Rea 
ad P 
» 
l—wA— PR ara’. IS 
aca 
Apply ing eq 6 to eq 10 gives 
R,—0 19 


Equations 12 to 19 may be solved for the eight 


constants M,, No, T+, R,, M,, N,, fe and R,. 


IV. Solution for Problem in Figure 2, c 


Reference [4] 
equilibrium and compatibility in a ring of constant 


shows that the equations of 


thickness are satisfied by the stresses 


0 Fk Ka?) A—3Ca'/r*—2Da?*/r*) cos 20 
og — F— Ka?/r* + (A+ 6Br*/a?+-3Ca'/r*) cos 20 
ra A 3Br a 3C at rt — Da?/r*) sin 20, 


(20) 


where fF, A, A, B,C, D, are constants to be de- 
termined from the boundary conditions. It also 


shows that the corresponding displacements are 
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u,=(r/ BF)i FA—v)— A114 v)a?/r? + [— AC +») | 
ty)at/rt+-2Da?/r*] cos 20) | 


‘ 


re=( BE\[AQ tp) + B(3+-p)r72/a+Cil | 


! yas r 


Dl—v)ja*/r*] sin 28. 


Applying the boundary conditions of eq 2 to eq 


fives 
l1—v) F.=(1+ vk, 22 
(1+-v).A, 2B, + (1+-¥)C,4-2D, 23 
(1+ yp). (3+ y)B (1+ v)¢ (l1—v)D 24 


Applying eq 4 to eq 21 gives 


(l—v)e’ F,— (14+ v) A, = (1l—v)e® Fi, — (14+) Ay, (2 
(1+ v)atA,—2va®R,+ (14+-»)C,4+ 207D, 
1+ v)atA,—2va®B, + (1+7)C,4+-2eD,, 2t 
1+»v)a‘A,+ (3+ »)a°B.+ (1+»)C,— (1—v)a2D, 
(1+ vata, + (34 vafB, + (1+ 7)C,— Ave D, 
Applying eq 5 to eq 20 gives 7 
Bo’ PF’, +- BA, a8 FF, + hy, 28 
Ba’ A, +3 BC, + 2 Ba? D, = ab A, +-3C,+2e°D,, (29 
Ba‘ A, +3 Bab B, —3 BC, — Ba? D, 
at A, +3a°B, —3C,—eD 30) 
Applying eq 8 to eq 20 gives 
fF, =P Awh, 31 
A PAwh, (32 
B,=0. (33 


V. Solution for Problem in Figure 2, a 


Combining eq 10 and 20, we obtain, for th 
stresses in a ring of constant thickness, 


o,=F-+ Ka?/r’+ > 
[—(3+-v)M/4ar— 7/22 — Rr/22\ cosd 4 
[|—A—3Ca'/r*—2Da* r*| cos 26 


0, F— Ka*/r’ 
[(1—v)M/4ar + T/2ar' +3Rr/2x\cos 64 
[.A+6Br?/a?+3Ca'/r'| cos 28, 


T /2xr 
3Cat/r't 


vy) M 4ar 
[A+3Br/a? 


+-Rr/2x| sin 64 
Da*/r*| sin 26. J 
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he values of the arbitrary constants, A, B, C, 
) KM, R, T for the sheet 1 and the reinforce- 
ent 2 can be determined as follows: M,, R,, F,, 





1, By: From eq 14, 19, 31, 32, and 33 
M,—P Bh 
R,—0 
F,=P/A4wh (35) 
A, P/Awh 
B,—0 


\f,: By combining eq 14, 17 and 18: 
M, Pih. (36 


r.. 2% By combining eq 12, 13, 15, 16, and 
17, and making use of eq 35 and 36, three simul- 
lancous equations are obtained: 


+2(3—v)a'R, 1-+-v)*a*P/Bh 
2(1+~)7,4+ 2(3—v)atatR,—2(1+-r)T, 

1 + v)?a*a*?(B—1)P/Bh 

287, +-2Bata'R,+-2T, = 0. 


K,, Fo, Keo: 


making use of eq 35 

l—v) F,—(1+»)k,=0 

|l—v)e’ F,— (1+ ») A+ (1+) Ay = (1 —v) ef P/Awhk 
Ba® PF, +- BA, — Ky =e P/Awh. 


By combining eq 22, 25, and 28 and 


(3S) 


C;, D,, Ae, Be, Co, Do: From eq 23, 24, 26, 27, 29, 
and 30, 


and making use of eq 35, 

1+») A,+ 2vB,— (1+7)C,—2D.=0 

1+») A.+ (34+ 9) B+ (1+)C,—(1—»)D,=0 | 

1 +v)atAl,—2va® BR, + (14+-¥)C,4+ 2D, 
1+») C, 20D, = (1+»)a'P/4uh 
(1+-vatA,+ (34+ v)a®B,+ (14+-¥)C,- 
1—v)o?D,— (1+»)C,+ (1—v) oD, 
1+viatlP/4Awh 

Bat A, +-3BC,+-2 Ba? D,—3C,— 20°), a P Awh 

Bab A, +-3 Ba® B,—3 BC, — Ba? D,+-3C,-4 
oD, o' P/Awh | 

(39) 
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Combining eq 11 and 21, we obtain for the dis- 
placements in a ring of constant thickness, 


u,=(r/E\|FU—v)—K(1+ va? /r*]4 
(cos 6) 4r FE | M("?—2v—3) log r+ 
T(1+-»)/P?+Rr(1—3r) +N] 

[(r cos 26)/E\[— AC +») —2vBr*/a’ 


C(1 + vjat/r*§+2Da?/r*| 


v—([(sin 6) 4r FE] M1 +»)? + M18 + 2v— r*) log r+ 
T+)? +-Rr* (5+ ¥) —N]- 


[(r sin 20)/ FLAG +) + B(3 + v)r*/a* 4 





C(1 + v)at/r§— DA — v)a?/r?| 
(40) 


The additional arbitrary constants, .V,, Vs, are 
obtained by combining eq 12, 13, 15, and 16 with 


eq 39 giving: 


N, 2(1+ vy) (a?—1)a*R,+[(1+)8- 
28(3—yv) log aa—2(3—y) log al 


lw >I Bh | 
(1+ v)P/2Bh | 4) 


1+-v)@?R,+[1+¥+2(3—y) log a] 


VI. Comparison of Analysis With 
Experiment 


The specimen, figures 3 and 4, was furnished for 
these tests by Chance-Vought Aircraft, Dallas, 
Texas It was of Metalite construction, consist- 
ing of two sheets of 0.013-in. 248-TS1 Alclad 
aluminum alloy bonded to a 0.5-in core of balsa 
wood, which had been cut with the grain normal 
to the sheet. The edges of the specimen were 
reinforced by pairs of *- by *- by !.-in. 245-T 
aluminum alloy angles bonded to the sheet. 

The concentrated load was applied by a steel 
pin, of 1.25-in. diameter, making a snug, but 
turning, fit in a hole reamed in a steel bushing. 
The Metalite was reinforced to a radius of 5 in. 
from the center of the loading pin. <A detail of 
the reinforcement is shown at section A-A, 
figure 3. The reinforcement to a radius of 2.5 in. 
was considered heavy enough to prevent distortion 
inside this circle. 
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The following constants were taken from. the 
specimen and w n the analvsis of thy 
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} O02 in 
1) Cbtety ON) 
t} 
he effective width em ot thre mnel was obtammed 
by dividing the gross cross section of the panel, 
neluding ede ing ry hie sheet) thickness, 


The arbitrarv constants were determined by soly- 


' 


i the simultaneous equations in section V, with 


the result 

A 0.3049P Al 0.1295P 
Bb O.00L977P 
0.06584P 
LD) 0.7572~P D) 0.1275P 
Fy —0.3049P FY, 0.1425P 
A, =0.426P A, — 0.0768P 


V4, —38.46P M,=15.15P 


402 





bheeebeeed 


DOREPPPE DED ROebeeeee 
SPEDEDPRDD DDD R DE DD ED: 


TiTiiiiii 


HRUURehabheee 


hee df 


aaaaae 


sehen LLL TT Tre eee 


CS LOMA RAE RP UP EAA ARR RAMA — 4. 





\ 20S 07 \ 63.01P 


O.00P15P 


R () R 
/ 69.04 P 


The stress in the unreinforeced region 1 
was computed by substituting the constants wit! 
subseript Ll into eq 54; whereas the stress in th 
reinforeed region 2 (2.5< r<_5) was similarly con 
puted using constants with subset pt » 

The radial and circumferential strains wet 
computed from the stresses by the well-know1 


formulas 
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/ 10.6 * 10° 1b/in, 


rure 5 shows the measured strains at a load of 
lb, and figure 6 shows a comparison between 


sured and computed longitudinal strains for 


the longitudinal and transverse centerlines and 


for a transverse section 12 in. below the center of 


the hole. The computed maximum strain, just 
outside the reinforcement, differed only about 2 
percent from the measured strain. Larger differ- 


as good as could be expected in view of the assump 


ences were found elsewhere agreement = Is 
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Picure 6 Average longitudinal strain at 10,000-lb load 
X, Measured. caleulated i) Transverse center line b) longitudinal 


center line ©) transverse line 12 in. below center of hole 


tions of the theory that a state of plane stress 
existed, that the central area was rigid, and that 
the hole and reinforcement were small compared 
to the width of the specimen, 


VII. Conclusions 


It is concluded that: 

1. The plane stress theory may be used to com- 
pute the stress distribution in the neighborhood of 
a reinforced circular hole loaded by a pin in a plate 
of sandwich construction, in which practically all 
the load is carried by sheets of metal bonded to 
each other or to a light-weight core and in which a 
relatively rigid central reinforcement is used, 

2. For large values of the ratio of plate width to 
hole diameter and with a nearly rigid central 
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reinforcement for applying the load, the solution 
can be expected to be accurate in all parts of 1 \e 
plate. 

3. For moderate values of plate width, the 
sults for the stresses are not of very satisfactory 
accuracy near the free edges of the plate, but ave 
reliable in the neighborhood of the hole near the 
places of largest stress. 

4. The presence of the reinforcement is likely \ 


move the point of highest stress from the edge of 
the hole to a point in the plate just outside the 
reinforcement. 
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Volumetric Determination of Columbium 
By Howard B. Knowles and G. E. F. Lundell 


The quantitative determination of columbium based on reduction by amalgamated 


zine and reoxidation with potassium permanganate has, in general, produced results that 


lack both accuracy and precision. 


A study of the various factors involved indicates that 


essential conditions require adequate duration of reduction in the presence of an excess of 


zine and a proper concentration of both acid and columbium 


\ procedure Is described 


whereby columbium, in a sulfurie acid solution, is quantitatively reduced by amalgamated 


zine to the trivalent form and subsequently oxidized to the quinquevalent state 


I. Introduction 


\ method for the quantitative determination of 
columbium, depending on the reduction in acid 
solutions of quinquevalent columbium to the tri- 
valent state by zine, and subsequent titration with 
a standard solution of a suitable oxidant, has been 
the object of numerous investigations extending 
over a period of years. In brief, the procedures 
attempted have involved the treatment of a hy- 
drochloric, a sulfuric, or a hydrofluoric acid 
solution of the columbium with zine, usually 
amalgamated, and subsequent oxidation of the 
carefully protected solution directly, er indirectly 
through use of a ferrie salt, with potassium per- 
With but a single exception [1]! no 
reliable procedure based on the foregoing principle 


Ihanganate. 


has been established that permits the use of 
stoichiometric values. A critical study of the 
procedures advocated by earlier investigators [2], 
in conjunction with the results obtained, is of 
little aid in indicating the underlying reason or 


cause for failure to obtain a quantitative reduction 
of this element. The use of material of unknown 


purity, imperfect procedures, doubtful conclusions, 


ind other factors such as inadequate data, have all 


served separately or collectively to retard the 


solution of a rather difficult: problem. 
The present article concerns only the question 
f consistently obtaining complete and uniform 


Figures in brackets indicate the literature references at the end of this 
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reduction of the columbium, such as would permit 
its use for quantitative determinations, and em- 
bodies results that have been obtained in an in- 
termittent study of various related factors. 


II. Experimental Details 
1. Amalgamation of Zinc 


A study of the reducing action of amalgamated 
zine [3] indicates that for a number of elements it 
is requisite that the zine be but lightly coated with 
mercury. Satisfactory amalgamation is accom- 
plished as follows: Transfer 1,000 g of 20-mesh 
zine metal (low in iron) to a stout, wide-mouthed 
glass bottle, add 500 ml of a 2-percent aqueous 
solution of mereuric chloride, mix continuously for 
45 to 60 sec, discard the solution, and wash five 
times with distilled water. Treat the washed 
amalgamated zine with 500 ml of warm, diluted 
sulfuric acid (1+-99)*, mix, diseard the acid, and 
wash thoroughly with distilled water. The pre- 
pared amalgam is preferably preserved in water to 
which a few drops of hydrochloric acid have been 
added. The period of use of an amalgam thus 
prepared is limited. With continued use it even- 
tually loses its efficiency and can no longer be used 
for this purpose. 


2. Sources of Columbium 


The two samples of columbium oxide used in 
the investigation were laboratory preparations. 
One was an oxide on which three different analysts 


? Unless otherwise specified, concentrated reagents are meant 
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Obtained 96.45, 96.58, and 96.42 percent of Cb.O 


by the tannin method [5 A preferred average 
value of 96.43 pereent for this sample would, for 
the purpose of reduction, be increased to an ap- 
parent 97.10 percent by the presence of 0.32 per- 
cent of TiO., 0.40 percent of red 
cent of Cr.O 
prepared from a specimen of columbium metal, 


and O.0S per- 


\ second sample of oxide was 


spectroscopically free from titanium and tantalum 
This specimen, however, contained traces of iron, 
copper, and tin, Which were removed before prepa- 
ration of the final Cb.O,;. This sample was used 


on the basis of L00.00-percent purity. 


3. Permissible Acidity 


\ study of the acid concentration necessary for 
the quantitative reduction of columbium indicates 
that the solution should) contain preferably 20 
percent by volume but never less than 15 percent 
by volume of sulfurie acid. With the latter con- 
veneral tendeney to slightly 


centration, there is a 


low but still satisfactory values. Under these con- 
ditions the reduced solution is clear and possesses 


color. At 


tions of acid the reduced solution exhibits a dark 


un “amethyst-like” lower concentra 
blackish) color and vields low values on titration 

The presence of small amounts of succinic acid, 
introduced as an aid to prevent hydrolysis, if not 


beneficial at least exerts no harmful influence. 


4. Concentration of Columbium Solution and 
Duration of Reduction 


Experiments indicate that essential factors for 
satisfactory reduction involve the time of contact 
of the solution with the amalgam, the amount of 
and concentration of the columbium to be reduced, 
and the amount of available zine. The last, in 
turn, is a function of both the active area of the 
amalgam and the degree of amalgamation. In 
this study the time for the passage of a// solutions 
through the reduector, unless otherwise noted, was 


not less than 25 min. No study was made of the 


rate of reduction, a factor conditioned by other 


variables: but where only a few milligrams of 


columbium ts) involved, satisfactory reductions 


can undoubtedly be effected in less time. Two 
Jones’ reductors of the conventional type contain- 


ng columns of amalgamated zine slightly less than 


1 inch in diameter were used. One column was 


IS in. in length and the other 33 in 


406 


Definite limitations regarding the maxim 
concentrations for use with these or similar red 
tors are difficult to define because the solutic 
when introduced into the reductor, undergo s 
dilution and, due to nonuniform mixing, la: 


homogeneity Although concentrations of as h 





as 1.7 mg of Cb.O, ml have been reduced quar 
tatively at times with the long reductor, a stu 
of all experiments indicates a concentration } 
greatly exeeeding 1.0 mg ml to be desirah 
Similarly, the concentration for use with the sh 


reductor should not exceed 0.7 to O.S hie mil 


5. Temperature During Reduction 


Although heat is not essential to quantitat 
reduction [4], all reductions were conducted at 
temperature of 65 a 4; This temperat 
facilitates the diffusion of zine through the am 
gam [3], and effects more vigorous action duri 


passage of the solution through the reduetor 
6. Results 


The results of experiments obtained with ty 
reductors under Varving concentrations ol a 


and columbium are given in the following tabli 


Prarie 1 Red j 1 ap ed / 
ts 
Hos 
Hy) 
‘ 
4 } Hy 1 i 
x | Hof) 
| 
. Hy 
12 H 
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? Red lior flected " eductor pith 1 S-tn. 


III. Discussion 


, 


A study of the results recorded in tables 1 and 2 


dicates that solutions containing 20 percent of 


furic acid by volume tend to vield. slightly 
sitive values, whereas those containing 15 per- 
nt are somewhat low. No study of this varia- 
on was undertaken, but in a few experiments 
i presence of traces of hvdrogen sulfide was 
dicated, notably in solutions of higher acid 
meentration, by a slight discoloration of moist- 


ed lead-acetate paper op contact with the eases 
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evolved. A similar case has been noted during 
the reduction of chromium solutions [6] 
titanium dioxide, as 


The presence of ndded 


recommended by some investigators [1], 7] as an 
aid in securing uniform and complete reduction 
does not appear to be a requisite. Experiments 6, 
25, 26, and 27 of table 2 vielded slightly low re- 
sults in the presence of this constituent but gave, 
on repetition in the absence of titania, equally 
unsatisfactory values, thus indicating other factors 
as a source of error 

A reconciliation of the differences between ex- 
periments 6, 7, and S with 20 and 21 of table 2 
and experiments 1] and 21 with 9 and 20 of table 1 
is difficult because of the faetors involved The 
low and unsatisfactory values obtained in these 
experiments are believed to be largely due to the 
amount and concentration of columbium that each 
reductor Is capable of handling, although other 
factors, such as a lower concentration of acid and 
the age of the reductors, may have exerted con- 
tributory effects A repetition of the experiments 
vielding low results, in a somewhat larger volume 
of solution, gave satisfactory results as is shown 
by experiment 31 of table 2 and experiments 22 
and 23 of table 1 
that for a satisfactory reduction it is necessary to 


A study of the results indicates 


maintain the proper relationship between both the 
concentration and amount of columbium to be re- 
duced and the size of the reductor that is to be used 

Experiments 6, 10, and 12 of table 1 indicate 
progressively the failure of an amalgam that had 
proved satisfactory in a series of 11 prior deter- 
minations involving a total of 2.66 ¢ of Cb.O 
The smaller reduetor, used in the experiments re- 
corded in table 2, gave no indication of exhaustion 
at the conelusion of a series of 14 determinations 
involy ing lL.46¢ of oxide 

With the reduetors employed in the present 
investigation, desirable limits appear to be 300 
me of Cb.O, contained in a volume of 200 ml for 
the 33-in. reduector and 175 mg ina similar volume 
for a reduetor with an IS-in. column. The above 
quantities may be materially increased if larger 
volumes of the Cbh.O 
the appropriate acid concentration, are emploved 


solution to be reduced, nt 


Solutions containing as impurities elements that 
are reduced to metal by zine should not be intro- 
duced into the reductor, since their presence tends 
to result in an unnecessarily violent evolution of 


hvdrogen 


407 








ree vy 





IV. Recommended Procedure 


The following procedure describes conditions 
suitable for solutions that may contain as much 
as 300 mg of Cb,O, and requires the use of a 
reductor having a long column of amalgamated 
zine. With solutions containing not more than 
175 mg of the oxide, satisfactory reductions can 
be effected in a shorter reductor, such as one 
with an 18-in. column. 

Fuse a dry weighed sample of not more than 
300 mg of the oxide in a 25-ml covered silica (or 
equivalent) crucible with 3 to 5 g of potassium 
pyrosulfate. Cool, transfer the covered crucible 
with contents to a covered 250-ml beaker, add 
20 ml of sulfuric acid, and heat moderately to 
dissolve the melt. Cool, add 1 to 2 ml of hydrogen 
peroxide (30%), cautiously dilute with 100 ml 
of water, and rinse and remove the crucible and 
cover. Cool, add 20 ml of sulfuric acid, dilute to 
200 ml with water, add 2 g of succinic acid, and 
heat moderately, while stirring, until the succinic 
acid dissolves. 

Reduce the columbium in a reductor having a 
33-in. column in the following manner: First 
clean and preheat the reductor by pouring through 
it hot diluted sulfuric acid (1+19) and then hot 
water (90° C). Discard the washings. Place a 
three-fold excess of ferric iron in the receiver and 
attach to the reductor. Pour into the reductor 
100 ml of sulfuric acid (1+4) heated to 65°+5° C. 
With the aid of suction draw the foregoing solution, 
the solution of the columbium, 150 ml of sulfuric 
3 During reduction, particularly during passage of the columbium solution’ 
the reductor alternately fills and empties due to formation of gas. This 
condition, although resulting in a nonuniform rate of passage of the solution, 


is normal and should not be disturbed. With volumes of columbium solu- 
tion exceeding 200 ml, the time of reduction is materially increased. 


acid (1+-4) containing | percent of succinic acid. 
each at 65°+5° C, and 200 to 250 ml of col! 
water slowly through the reductor without inter- 
ruption at such a rate (unchanged) that the tote! 
time for the passage of all solutions is not less thay 
25 min.’ Cool the solution of reduced columbium 
by immersion of the receiver in ice water during 
the operation. Remove the receiver from the 
reductor, rinse and drain the outlet tube in the 
customary manner. Add 10 ml of phosphoric 
acid (85%) and titrate with a 0.1-N solution of 
potassium permanganate (0.0066455 ¢ of Cb.O, 
ml) [S], using orthophenanthroline as an indicator 
Correct the titration by deducting the amount of 
potassium permanganate consumed by a_ blank 
determination, involving all reagents, obtained 
under identical conditions. To insure absence of 
oxidizable impurities each of the  acid-wash 
solutions is permanently tinted, while hot, by 
the addition of a few drops of a solution of 
potassium permanganate. 
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Effect of Concentration on the Viscosity 
of Dilute Solutions’ 
By Robert Simha 


A theory of the higher coefficients of the viscosity-concentration curve of a suspension 


of nonspherical particles is presented. 


Starting with the simple model of a dumbbell, the 


flow around a single particle and its modification due to interaction with other particles are 


considered, 


It is shown that the coefficients a,; in the equation 


n=mno(1-+-a,e+ az,ce*+-a,c°+-. ...) 
are related to each other, namely, 
a,=k,a,?, a3=k,a;*. ..., 
where the k's are independent of molecular size, in agreement with empirical equations. An 


explicit value for k; is obtained that applies also to models consisting of various arrays of 


spheres. 


The variation of k; with molecular shape is discussed, and various factors affecting 


the numerical values and the validity of the last-mentioned equation are pointed out. 


I. Introduction 


The dependence of the viscosity » of a dilute 


solution on the concentration ¢ (expressed as 
volume fraction or as weight by volume) can be 


represented as a power series ot the type 
9=0( 1-+-a,¢+- G20" a3¢°+-....); (1) 


if the solute is a nonelectrolyte. — no is the viscosity 
of the pure solvent. Considerable experimental 
effort 


determination of the coefficient a, 


and theoretical has been devoted to a 
the so-called 
trinsic viscosity, for suspensions of particles of 
shapes 1]! original 


hvdrodynamie treatment of spheres. In 


various since Einstein’s 
these 
livdrodynamie theories, a; is a measure of the 
disturbance of the flow of the solvent due to the 
presence of the solute particles at infinite dilution. 
The total effect under such conditions is the sum 
of the effeets of the individual particles. For the 
range of higher concentrations, a number of 


Figures in brackets indicate the literature references at the end of this 


warded an A. Cressy Morrison lrize in Natural Science in 104s by the 


York Academy, 


of sciences 
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empirical equations have been proposed. They 
can all be expressed as a power series of type given 
in eq 1. As one can readily see, the theoretical 
calculation of the higher coefficients a, a3, ete., 
becomes increasingly complicated. It has been 
carried out explicitly first for the coefficient a, of 
a spherical suspension by the author [2]. a, a 
number independent of size for spheres, measures 
the disturbance of the flow around a particle by 
a second one at a finite distance from the first. 
The method applied in the particular case is 
analogous to the procedure first developed by 
Smoluchowski [3] to obtain the modification of 
Stokes’ law for the resistance of a sphere caused 
by the presence of a second one. Smoluchowski’s 
results, of which we shall make use, are obtained 
by a perturbation in terms of inverse powers of the 
distance between the two spheres. 

Empirically it has been suggested [4, 5] that the 
coefficients a, and a, are related in polymer solu- 
tions by the equation 


a> ke,*, (2) 


where the factor / is independent of the dimensions 
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of the solute particle but may vary with the 
solute-solvent In 
although not always, this rule seems to apply 
It is the purpose of this paper to 


system. several instances, 
satisfactorily. 
examine the validity of the proportionality (eq 2), 
using first the simplest model for a nonspherical 
particle, namely, that of a dumbbell considered 
originally by Kuhn in his treatment of intrinsic 
Viscosity 16]. This extended 
that of a pearl string currently used by various 
authors in the of hydrodynamic 
properties of chain molecules. 

In section IT the velocity distribution around 


model ean be to 


description 


a single particle and the intrinsic viscosity a, are 
obtained because the results are needed for what 
follows. In IIl and lV, 


effects between solute molecules are considered. 


sections interaction 
In section V attention is devoted to the nature of 
the approximations used and to the effect of vari- 
ous factors on the validity of eq 2 and the magni- 
tude of the parameter &. 


II. Flow around an isolated particle. 
Intrinsic viscosity 


Figure 1 illustrates the notations and reference 
2L isthe length of the dumbbell and 
a the radius. It will that a/L<1. 
Let the flow of the pure solvent be represented by 


the X,.\,-plane, 


axes used. 
be assumed 


a simple shearing motion in 
namely, 


i qNo, Me u.” 0, (.>) 


where q is the \v elocity cradient and uy the velocity 
component parallel to the Y-axis. If the dumbbell 
is held in a given orientation relative to the refer- 
the centers P and ) have the co- 
ordinates 2,(P) and a,(@ r,.(P), respectively. 
Equation 3 then gives for the velocities at P and 


u,(P), 


ence system 7 " 


u,(P)= qr P); u,(Q) 


(3a) 


where the z,(2”) are functions of the length 1 and 
the direction cosines of the particle axis. In the 
absence of rotatory Brownian motion, this flow 
will rotate the particle axis into preferred direc- 
tions. If one superimposes upon the motion (eq 3) 
a second one resulting from the spins of the par- 
ticle and considers the resulting force and torques 
acting on the spheres P and @, one arrives at a 
result’ corresponding in three 
Kuhn’s [6] and identical with Jeffery’s [7] formulae 


dimensions to 
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Dumbbell model and 


geometru 


FIGURE l. reference used ar 


ares 


relations involved, 


for the rate of change of the orientation specializ d 
for a thin rod. 
cerned only with the case of so-called complet: 


Here however we shall be con- 


Brownian motion, where any orientation tendency 
is absent. Following a previous procedure |S] wi 
shall consider only the effects of the velocities 
eq 3) and (eq 3a) on the particle. 

The change in the unperturbed flow (eq 3) du 
to a single dumbbell is obtained as a solution of th 
hydrodynamic equations of motion, which neglect 
ing inertia terms, assume in this case the form 


| 
mee | 


oN, 


Here, as in the following, summations from om 
to three are performed over twice occuring indices 
m. p is the At 
vu, has to assume the form shown in eq 3. I! 


pressure, infinite distance /’ 


the velocity u, is transformed to a reference sys- 
tem .Y,(/’) parallel to the first but with its origin 


at 7, such that Y,—N.(P?)+.7.(P), it is seen that 
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he velocities u, consist of two parts if the mutual 
listurbance between the two spheres is neglected. 
The first satisfies the condition imposed by the 
\,(P?)—term. The required solution is of the 
type given by Einstein [1]. The second part, 
satisfying the constant velocity qz,(/’) is of the 
type given by Stokes [9] for translational flow. 
Now for the calculation of the dissipated energy 
of the system particle plus solvent, it is sufficient 
to consider the velocity and pressure distribution 
at large distances compared with the dimensions 
of the solute molecule, and therefore a fortiori 
large compared with the radii a of the spheres. 
As Einstein’s solution contains terms of the order 
(multiplied by a velocity gradient) and 
higher, while in Stokes’ solution terms of the order 
app (multiplied by a velocity) are found, the 
former can be neglected for our immediate pur- 
poses. Physically this amounts to concentrating 
the effect of each sphere on the solvent flow in its 
center P or Q@ The calculation of the perturbing 
field of flow then bears a certain formal similarity 
to that of the field of a rigid dipole, and including 
to the calculation of 


wpe 


hydrodynamic interaction, 
dipole-dipole interaction. It be remem- 
that at the 
particle the equations given below are invalid. 
For the disturbance due to sphere 7’, we find in 
this manner from Stokes’ solution neglecting terms 


should 


bered, however, distances close to 


of the order a*/pp*: 


3 a 3 aX,(P)X,(P) 
u,* u,(P) —6,,——~ u,(P) ; — 
4 Pp 4 Pp 
(5) 
3 AAP) 
p=const.— 5 mu)(P)a et, | 
- PP” } 
where 
» 3) sel 
“1 o6 1. 
The additional velocity u® around sphere Q 


is given by an analogous equation. From eq 


sa and 5 it will be seen that the total velocity 


obtained by superposition of u,(2?’?) and u,(Q) 
will contain terms 
| l 
Pp Pe 
and 
X,(P)Xi(P)_ NX,(Q).X (Q), 
pp* Po 
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At points in the field of flow for which not only 
pe, pera, but also R>L, it is permissible to 
expand these quantities in terms of the ratio 
L/R (see fig. 1) and to discard the third and higher 
powers of L. Omitting details, one finds 


& 
1 1 
a ~ FS 2 , ante 
NX(P)Xi(P)_. V¥i(Q@).X(Q) 
pr Pa 


r (6) 





6.Y ae .. = o 5 r 
R 2s Nire(P) ~ RS [.X,2,(P) + V,2,(P)] 
NX(P)_ Ni(Q) 6X, Sy. ~—— 2r,(P) 

pp py R* oa Nite t ns Rs } 


The velocity to be added to eq 3 is then given by 


uy = ne XP), 
[.V,2,(P) + A alr )}+ a BY 2.(P (5a) 


,  o, U(P)aP.XySON2(P) 
Pp const.— snp Re Rk? 2(P } 


Equation 5a is, of course, a solution of the equa- 
tions of motion (eq 4). It is useful to write it 
in the form of Lamb’s [9] general solution in terms 
of spherical harmonics. It is easily found that 
eq 5a corresponds, using his notation, to: 


3 a ) 
» 


P-s=P, x-2 [N243(P)— X32(P)]. (5’) 


the 
flow (eq 5a), one considers a large 


In order to obtain 


the 


dissipation of energy 
caused by 


; - Ar ; 
spherical volume | » Re around the particle. 
» 


Let p,y and p,y, denote the components of the 
radial stresses exerted per unit area on the surface 
of this sphere. Then one derives from eq 3 and 


the definition of the stresses: 


, Xs. Mt. 
Prx, Nod R Pex," = 09 R, Prx,’=0. (3b) 


Equation 5’, together with an expression cited in 


{9] vield 


(2.X,? 
( R 


1S gl (P)a —T . - : ) a 
Prx = PRAGA: ) Nin (P) (7) 











The rate of dissipation per unit volume is given by: 
id ; P 
( Pex t Prx.)(Um Un dS (S) 


where the integral is extended over the surface of 
the sphere Ry. The 


the order 


result is, a terms of 
1 fy and using eq. 3, 3a, 3b, 5a, and 7: 


dW SOm ar, st 
= nod 1+ (Sa) 


If the 
solution consists of » particles dissolved in the 
such that NVe=n/V 


average number per unit volume, 


if the particle is held in a fixed orientation. 


volume V, represents their 


one has 
2° a@N=c, 9) 


In the pres- 
and hence the 


where ¢ is the volume concentration. 
ent case all effects are additive, 
dissipation of the whole solution is obtained by 


combining eq. Sa and 9 to 


dw 27 2a Py 
dt Nod (1 t 10 a c) (Sb) 


In order to obtain the v ISCOSILV 7 of the solution 
we follow the method of reference [1] It will be 
assumed that the solute molecules are randomly 

Let particle i have 
&. (7) in the X,-system, 

Its contribution u,(7 


distributed in the volume V 
its center situated at a point 
eel, 2, 3;3221,2,...28 

to the flow disturbance at a point XY, is obtained 
from eq 5a if the coordinates .\, are replaced by 


the relative coordinates NX,- Ey i) Hence we can 


write for the velocity distribution in the whole 
solution 
u é q.X, Sou : 5b) 


The \ ISCOSILN of the pure solvent can be defined by 
the relation 


‘ed ly nt Ou, Ou, 10 
"Dt ee, | ae A 


In the same manner we regard the solution as a 


liquid with effective Viscosity n subject to an 


‘fective gradient q, defined as 
dW Ou Ou, ; 
( ”n i | z 
dt ™* ~"| ox, TON, Jy 
\ ,| OU,(7 Ou» i } 
=n , S ee (10a) 
14 >| oN, 0X, |. \ on 
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where the left-hand side is given by eq 8b. A 
0/0.X, 0/0 (7), the sum in eq 10a may b 


replaced by 
+. ,] Ou p Ota! i) . 
| E, (4) 0g, (i) XY, . 


@) 
Ou; 
N LS Of» toe ar 0 


if an integral is substituted for the sum. Th 
integration is to be carried out over the whol 
volume V. The expression can be transformed 


into an integral over the surface S, namely, 


N[ (a ptt p, as 


Substitution from eq 5a then leads to the result 


9 r,(P)? 


Thus the viscosity » of the solution follows from 


9 Fr, we ue 1 


gives the correct viscosity up to linea 


eq 10a as 


: 27 r,( P)? 
” ” 10 «a 
Equation 12 


terms in the concentration ¢«. Averaging over-al! 


(random) — orientations, 2,(P)*)=(L?/3) wi 


finally obtain 

3 
2 a®' 
The intrinsic Viscosity for the model under con 
sideration is proportional to the ratio L*/a? i 
concordance with the familiar result obtained fo: 
infinitely thin ellipsoids or rods not subject to 
any orientation effects by the flow gradient. 

The quadratic and higher terms in ¢ obtaine: 
from eq 12 have to be discarded as long as hydro 
dynamic interaction effects that give contributions 
of the same character and order of magnitude {2 
This has sometimes 
when Ein 


are not taken into account. 
not been observed in the literature 
stein’s results have been extended to higher con 
centrations. However, it can be seen that th 
higher terms in eq 12 will, after averaging th: 


powers of z:(P), contain successively higher powe! 


? The procedure given in reference [7] can equally well be made use of wt 


nteraction effects are absent 
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ratio L*/a*, thus confirming the empirical 


the 


ition (eq 2) which may now be written as 


a,—k,a,?; a,=k.a;’... ., (2a 


re the k’s are, for the model, 


mbers and independent of the dimensions of 


present pure 
molecule 

Equation 13 could, of course, have been derived 
a much more straight-forward manner by con- 

dering the frictional to the flow 

| offered by each of the two spheres. The rate of 

ssipation arising in this way is given by the 


resistance 





oduct of frictional force and relative velocity. 
If the latter quantity is taken from eq 3a, we have 


vy application of Stokes’ law 


dw 


dt n q t2 Y “| N. 


oan aq [a r 
Klimination of VV by means of eq 9 leads at once 


req 13 
III. Interaction between two particles 


The influence of a second particle on the Viscous 


low around a first has been considered first: by 


Smoluchowski [3] for spheres in translational mo- 


on, leading to a modification of Stokes’ law 


The effect so obtained is a function of the ratio 
between radii and mutual distance of the two 
spheres. More exact calculations are reviewed in 


Oseen’s treatise [3]. Rotating spheres have also 
been treated [2] and calculations for nonspherical 
particles have been presented [3,10]. The principle 
of Smoluchowski’s method consists in seeking to 
fulfill by the hvdrody- 
namic boundary conditions on the surfaces of all 
solid That is, 


distributions are superimposed on the previous 


successive “reflections” 


bodies involved additional flow 
flow, to compensate for the additional disturbance 
at the boundary between particle and fluid in- 
Thus the 
boundary conditions are fulfilled to successively 
The first 


proximation that is sufficient for our present pur- 


duced by the neighbors of the former. 


higher degrees of approximation. ap- 


u, U,OU~™(é, U,% 
3 to — $i) 1a 
, Nol } 
l 2 R; 
where 
o¢@. 3a 
{ é » O07 > 
1 PR, is 
Che author ndebted to E. O. kK heckit ji4 
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pose has been shown [3] to fulfill the boundary 
conditions not at each point of the surface, which 
but the 
The velocity 


would make it an exact solution, on 


the 
distribution is in other words such that, if aver- 


average over whole surface. 
aged over the surface of the sphere, it gives no 
relative motion between the sphere and the liquid 
at their boundary. The solution previously em- 
ploved by the author [2] in the viscosity treatment 
is of the same degree of accuracy. 
Smoluchowski’s original equations refer to 
spheres moving in the same direction with equal 
speeds. His results are modified slightly if this 
equality does not hold. Let the position of the 
centers of two spheres be given by &“? and &@ in 
a coordinate system fixed in space, their mutual 
distance by Ry and their respective velocities in 


the U,™. U,° At a 


fig. 2) the velocity and pressure distribution pro- 


£,-direction by point &, 











g 


FiGuRE 2 Geometric relations in interaction of two sphe re 


eq salions 14 and 15. 


duced by sphere 1 at large distances R, has the 
following form: 


(14 
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1°, (&) represents this function at the point 
z,**, 

An analogous expression results for wu; and p®. 
a each consist of two 


(1) 


It is seen that uP and p 


terms. The first is, of course, Stokes’ expression. 
The second one is a function of the mutual distance 
Equation 14 


and represents the interaction. 


leads to a frictional resistance 


. » (2 » » (2) - (1) r°Q - 
+R, En) gi) U, ] (15) 


Thus the resultant resistance is no longer in the 
direction of motion. 

Let us now consider a second dumbbell (VM, N) 
as shown in figure 3. Relative to its center the 
points MW and N have the coordinates 2,(./) and 

ry( MD), which are functions of the orientation of 
the axis ZN in respect to the X,-system. They 
u,(\M), the re- 
spective velocities of the undisturbed solvent flow 
relative to the particle at points M and N. As in 
section IL we shall be concerned with distances R, 


determine u,(\MJ) and wsj0N) 


which are large in comparison with L. Consider- 
ing furthermore only molecules that are far apart, 
the modification of the Einstein-type of motion 
which gives rise to terms of the order a® [2], will be 


With these 
approximations one can make use of eq 14 and 


omitted and only a’-terms retained. 


write the velocity around 7’? to be superimposed 


P = 
on wy In eqo 


Au,”=u,(M)U,,"°(M,P) . Ui? (X, P)+ 
w(N)U ma? (N,P) Uy? CX, P) 





with 
(lf 
BY D 3 a \ 
U,*(X,P)= 4 pp et 
[V,—2r(P)}[Xi—2(P)])- 
pr ) 


Uy’ (M, P) is obtained from the last equation by 
replacing pp by ppy and the coordinates Y, by 
&,+2,(M). 
corresponding Manner, An expression analogous 
to eq 16 results for Au;®. We define the inter- 
action coefficients Pa and qa by means of the 


The last term in eq 16 is obtained in a 


relations: 


pa=Uy?(M, P)—U,P(N,P): qa=U:%(M,Q 


U,%¢N,Q) (17 


and obtain from eq 16 for the total velocity dis 
tribution around the molecule P, Q 
+ (MM) pay Uw? LY, P) 4 


u, =u, (Stokes) 


uy M Mal Fan NV,Q) (l6a 


ry 
rhe re- 


mainder contains functions of the mutual dis- 


where the first term is given by eq 5a. 


tance and orientation of the two interacting mole- 
cules. As stated previously, one sums over 1 


from one to three. 


The quantities (",,,” and (’,,.% in eq 16a are of the same type as Stokes’ velocity distribution and can 


be transformed in the same manner for large distances RP. 


ee 4 3au(M) \ 
/ StoKes) (Pmi—Jmi) 
i 41 eR Pmi—~Jmi) y ; 


3 au,(M) 


p— p(Stokes)+ , R (Pni— Om) 


In analogy to eq 7 this makes a contribution to the Y,- 


9 u,(.M) 


Prxi = Prxi( Stokes) : Noll R 


Nytu(P)bim- 


> Xite(P) tn(P)\- 


(4X, m 
(Pmi— mi) ) 


One finds in this manner in analogy to eq 5a 


Re "D> X,2(P)—[X,.21(P)-4 Nita(P)]\- (18) 


<> ee 
, 


LY,.2(P)-4 R? 


Xizn(P)|+ St Xi2(P) 
k 


component of the radial stress 


Application of eq 8 now yields for the total rate of dissipation per unit volume of the system consisting 
of solvent plus one particle (2, Q@) acted upon by a second one (M,N): 


dw 
dt 


36 ra 6 ra 
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J » 
ul a 1+ 5 V x,(P)?— - \ £2(.M) [82.(P) (puna) - 2n(P)(Pu—da)Iy° (19) 
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Creametric relations in interaction between two 
molec ules, equation lt 

The interaction coefficients py, and gy can be 

evaluated for pairs sufficiently far apart in terms 


of the mutual distance p (see fig. 3). A’ series 
expansion in terms of 1/p gives, using eq 17 and 


disearding higher powers in L/p: 


oa \ 15 ‘ ‘ 
Pa—Qa)* — 6 HE Doge (P)- Doge (M 
plop 
3 Sor.( Pye M)—"y (és (P)+ &0(P)). 
p" — p : 
Sen (M) =") (2M) +8, MDa P) 


- tm, 


r,(P)a2(.M)+2,(P)2,(.M) +6 E > bell). 
p 


4 


[ 


‘t,7,(AZ) S(P)e(M ]}. (20) 


Equations 16b and 20 constitute the complete 


solution of the interaction problem for two 


particles, 


IV. Viscosity of Solution 


To obtain the quadratic concentration term in 
the Viscosity one computes first from eq 19 the 
dissipation of the mixture. If we visualize again 
the centers of the solute molecules randomly dis- 
tributed over the volume V, then for the ith 


particle (P,Q, situated at a point &(/) and acted 
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upon by a fixed particle (WM, N), pal)—gqal) 
becomes a function of the relative coordinates 
&—&(i). It is to be summed over all (P,Q, 
species interacting with one (V7, N) and then mul- 
tiplied by the total number of particles. In this 


manner eq 19 vields 


dW dW 
Ci) 


dt @ Fou 
9 ro(.M) a > . . , 
ayy moe? (34 PY [Puli gu (i) 
27,(P)S5 [palt)— al], 19a) 


where the first term is given by relation (eq Sb) 
The summations may be replaced by integrations 
over the volume V. The use of the expression 
eq 20 in the whole region amounts to neglecting 
the influence of particles close to a given one. 
The error so committed increases, of course, at 
higher concentrations where the average distances 
between solute molecules become small. The two 
integrals can easily be transformed into surface 
Omitting 
all terms that vanish on integration because of 


integrals by means of Gauss’ theorem. 


reasons of symmetry, one finds for instance, that 


the remainder in (p.,;—q) can be written as 


(ou—%: 90 Oy, (P)te(M)-+24(P)( MD} 


O/& S67), 0 ff 3b *k 
| fi: p )+ ae ( p” ; )| 


with a similar expression for p,,;—q,. The final 
result for the dissipation due to particles oriented 
ina given direction «;(/) interacting with particles 


in another direction x,(.V/) is 


dW o Sa 27 af PY 
dad ™  'ITi9 @ | 


0a! [4 (Pay Pyar Mary M) (19b) 


622( P)ao(.M) (22( P)a( M) 4 


r4(P).ty(.M)) +, ( PYM] 2! 


) 

The next step consists in calculating the effec- 
tive velocity gradient g, defined in eq 10a. The 
velocity distribution around molecule ¢ in’ the 
presence of a perturbing molecule at an arbitrary 
point & is derived from eq 16b and 20 by re- 
placing .Y, by XY, 


£.(/) and & by &—&(4). 
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The corresponding result for the whole mixture 


follows again from eq 5b. 
section IT. 


q, is then derived as in 
The expression 


to be integrated over the sphere 2, now involves 
& ¢ t 4 

terms of the order >,, and 3). The integration 
R R,° 


requires the values of the following integrals over 
the surface of the spherical enclosure: 


oe 


RY; Sett"dS= 


~! } AS S 2r7t"dS _ R,*. 


¢ ¢ 
si $2 63 - 
105 


The integrand written above depends upon the 
mutual orientation in each molecule pair. On 
r(P) and r.(M) 
terms, there remains after a lengthy but straight 
which 


averaging independently the 


forward caleulation, omitting all terms 


vanish for random orientation 


wNS ( u S: +u 1 )ds 


R - a 
. 162 a- > 
nA 5 rq R z,( P)*z, M)?. 


Instead of eq 12, one finds now, using eq 19b, that 


27 z.A{P 27 7 x.(P)*x.(.M)* 
10 a 10 5 a’ 
7 ™P,_ 9 2(P) 9 27 2(P iM)? oP 
10 a- , 1 4 a’ . : 
(12a) 


where the denominator represents q,. Expand- 
ing up to quadratic terms in the concentration 
and averaging over the direction cosines, we arrive 
at the following result for the viscosity of the 
solution: 


3 L? 87 Lt , a 
n m (1+5 2 OT en ai ¢* ): (loa) 


On the other hand, eq 12 would have given for 
the coefficient of the quadratic term the smaller 
value 81/100. 
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V. Discussion 


Equation I3a yields for the coefficient /, define 
in 2a the value 
kk, =0.77. (2) 


\< 


Had interaction effects not been taken int 


consideration, eq 12 would have given 
&,=0.36. (2la 


These results refer so far to the model of a dumb 
bell. Let the molecule be represented next by a 
linear array of length 2, consisting of a numbe: 
of spheres with radii a, v per unit length, arranged 
evenly. The simple considerations at the con 
clusion of section II, which lead to eq 13, indicat 
that now the factor L * 


value (Li< L’. 


contribution to the dissipated energy by a spher 


is replaced by an average 


This is seen as follows: The 


located at a distance y from the center of thi 
molecule is proportional to y*. This is to be 
summed over all groups in the molecule and 
multiplied by .V, the number of molecules in unit 
volume. Nis proportional to 2vL 


Thus the intrinsic viscosity a, 


inversely 
becomes propor- 
tional to 

2p L , 

QoL yfdy=(L 


0 
« 


In the case of a randomly coiled chain one aver- 
ages over all internal configurations of the chain in 
addition to averaging over all orientations of the 
whole molecule. As long as coiling is not extensive 
the determining quantity has been shown re- 
peatedly to be the mean square separation ol 
Reverting to the linear array, it can 
be seen from eq 19 and 20 that the interaction of 


chain ends. 


parts of two molecules at distances y; and y, from 
their respective centers will contribute terms pro- 
portional to yy." to the dissipation rate. Con 
sequently, the coefficient (ly will be proportional to 
the square of (L*) as one could have surmised 
A corresponding result obtains for a solution con- 
Thus the 
numerical values (eq 21 and 21a) should be valid 


sisting of moderately coiled chains. 


also for a long extended molecule and a coil. 

The value (eq 21) is larger than those derived 
experimentally from suitable extrapolations of 
viscosity-concentration curves, the latter being 
usually of the order of 0.3-0.7 and more nearly 
comparable with the value (eq 21a). On the 
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her hand, the theory in its present form can 
produce only approximately the contribution of 
‘hydrodynamic interaction to the &, coefficient 
rather to its excess over 0.36. It has been 
essed previously that the evaluation of the 
teraction coefficients py, dp (eq 17) as a series 
4 20) and therefore eq 19b are valid only for 
y distances. One could im- 
ove the situation by dividing the volume V 
round the center of a molecule into three parts. 
The first comprises all particles with spherical 


rge intermolecular 


vroups so close to those of the center particle, that 
an expansion of the mutual disturbance in terms 
of the ratio @a/Ry (fig. 2) becomes altogether 
In other words, solute molecules are 
their 


meaningless. 
touching each other whatever 
From a hydrodynamic point 
these complexes would be treated as 
Such conditions, frequent at 


higher concentrations for large molecules will be 


practically 
relative orientation. 
of view, 
single molecules. 
disregarded here (see also below). In the second 
region intermolecular distances are large in com- 
parison with the size of the individual group in the 
molecule, but are of the order of magnitude of the 
dimensions of the molecule. This would permit 
the application of Smoluchowski’s equation (eq 
14), but not of the expansion (eq 20). All 
further summations are to be performed in terms 
of individual group to group distances without 
using expansions in terms of molecular center to 
In the third region the current 
Of course, 
division is somewhat arbitrary. 


center distances. 
approximations are valid. the sub- 
The conditions 
in region two affect both the numerator and 
denominator in eq 12a. They should be more 
important in the caleulation of the dissipated 
than in that of the effective 


velocity gradient g, (compare eq 10a), because of 


energy (see eq 19a) 


the direct summations over the (py—q,) oceur- 
ring in the former. Also in this region a disturb- 
ance of the random mutual orientation of particles 
merely because of geometric reasons is possible. 
The treatment in section II took no account of 
the hydrodynamic 
in the same molecule. 


interaction between 
For the dumbbell a simple 


ealeulation based on the result (eq 15) leads in a 


groups 
lirst approximation to a correction factor in the 
expression (eq 13) for a, equal to 

9 a 
1 
t20 L 


Viscosity of Dilute Solutions 


fora/L<1. Alsoasafirst approximation, a, remains 
unaltered. Thus only a decrease of k, by a few 
percent results from this effect. In an array of 
spheres the interaction of neighboring spheres 
with velocities of the same sign relative to the 
surrounding liquid will actually give a correction 
factor that is smaller than unity. The inter- 
action between different groups within a strongly 
coiled molecule changes the flow pattern pro- 
nouncedly, as has been shown recently by several 
authors, and this influences also the total mutual 
disturbance between two different molecules. 

It is perhaps noteworthy that the calculation in 
section II already gives as large a value of hk, as 
0.36. Effects other than purely hydrodynamic 
ones may contribute to the slope of the viscosity- 
concentration curve even at relatively moderate 
concentrations. Experimentally the reported ini- 
tial slopes have sometimes been obtained by extra- 
polation at concentrations at which the conditions 
inherent in our derivation are no longer fulfilled; 
the solute molecules of large molecular weights 
are no longer sufficiently far apart on the average 
There are indications in the literature that at 
concentrations roughly below 0.1 g/100 ml solu- 
tions of chain polymers exhibit steeper slopes than 
the extrapolated straight lines. Qualitatively 
similar effects have been observed in tobacco 
mosaic virus solutions [11]. In the latter case one 
is dealing with long rod-like particles and the 
possibility of aggregation is strongly increased at 
higher concentrations. Orientation effects are then 
relatively more important and this reduces the 
viscosity at higher Systematic 
measurements of the dilute 
polymer solutions as function of velocity gradient 


concentrations. 
Viscosity of \ ery 


are desirable to decide the 
curvature, 


and soivent character 
nature of the apparent downward 
found on approaching the region of high dilution. 

It is obvious that the present considerations cdo 
not apply directly if the nature or the identity of 
the individual molecule changes with concentra- 
tion. For if formation of aggregates 
occurs Which can be expressed, say, by an equi- 
librium constant A between double and single 
makes an 


instance 


molecules it is seen that this 
additional contribution to the coefficient a,, namely, 


LK (a,° a; ) 


where a,“—a, is the difference between the 


intrinsic viscosity of double and single molecules. 
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Thus the higher coefficients in the expansion 
(eq 1) cannot be found solely by hydrodynamic 
considerations and therefore eq 2a is then not 
valid. This is particularly the case when the 
effect of velocity gradient in concentrated solutions 
examined. even in dilute 
solution the interaction is affected by 
orientation than the intrinsic viscosity. 

Finally the 
particles may be compared with those previously 
obtained for compact The author's 
method of calculation [2] leads to a value a,=14.1 
[10]. Einstein’s results give a value of 4 for ay. 
Thus, for spherical particles, the corresponding 
values are 4, =2.26 and 0.64, respectively. It 
would be of interest to caleulate the gradual 
increase in the #,-values with increasing sym- 
metry by considering ellipsoidal particles of finite 
thickness. Such results would be of value in the 
case of protein solutions and also of chain poly- 
mers. One should expect a parallelism between 
ky and other parameters characterizing shape 
such as the exponent in the modified Staudinger 
Also cross-linking 
should result in an increase of kh). Such changes 
have been reported. It appears from the fore- 
going considerations that, leaving aside other fac- 
tors, the parameter &, will remain independent of 
molecular weight as long as the degree of coiling 


However, 
more 


is to be 


present. results for asymmetric 


spheres. 


equation. branching and 


dloes not change markedly. 
VI. Conclusions 


In treating the concentration dependence of 
Np, two effects may be considered at low concen- 
trations. Aggregation, or entanglements, and 
hydrodynamic interaction of distinet molecular 
The first makes a contribution to the 


viscosity of the solution, whieh depends on the 


entities. 


frequency of the aggregates and can be described 
in terms of the intrinsic viscosity of an aggregate. 
The second one leads to results that are in agree- 
ment with empirical relations proposed for dilute 
Withour explicit caleulation of the 
cubic and higher terms in ¢ it cannot be seen 


solutions. 
whether »,, is represented by the binomial for- 


mulas given by various authors. The present work 
indicates that they probably are not strictly valid. 
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The numerical value derived for k, (eq 2a and |) 
is of the correct order of magnitude but somewhat 
larger than those encountered for solutions of chein 
molecules in good solvents, to which the mocel 
used is applicable under certain conditions. — [t 
seems to be more nearly in agreement with some 
numerical results for globular proteins. Possible 
experimental and theoretical reasons for this 
divergence have been outlined. 

The parameter /, depends in a characteristic 
manner on the shape, as determined for coiling 
molecules by their structure, the solvent environ- 
ment, and also on the molecular weight. This 
conclusion can be reached on the basis of the 
results obtained for two extreme cases treated here 
and in previous work, respectively. In this con- 
nection it is noteworthy that globular proteins 
appear to give larger values for k, than chain 
polymers. The region between the two extremes 
of an “open” coil, or dumbbell, and a sphere re- 
mains to be treated. On the basis of such a 
treatment it may become possible to derive 
information about the properties of the solute 
similar to that presently obtained from intrinsic 
viscosities. 

It is of interest to examine the influence of the 
mutual disturbance on the viscosity at finite 
velocity gradients, when partial orientation in the 
field of flow occurs. Effects present at infinite 
dilution should be more pronounced. 
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